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delineated.  A simple  equation  is  given  for  the  threshold  of  thermoelastic 
instability  when  both  faces  are  initially  flat. 

For  the  case  of  an  initially  wavy  conducting  face,  the  model  for 
viscous  heating  envisages  a full,  coherent  film  on  the  one  hand  and,  on 
the  other,  a broken  film  with  zero  heating  in  the  cavitation  zones.  For 
both  models,  the  continuous  thermoelastic  deformation  of  the  conducting 
face,  with  changing  sliding  speed,  is  discussed  and  criteria  for  the  onset 
of  instability  are  presented. 

txperiments  are  reported  to  corroborate  the  findings  of  the  theoretical 
analyses. 
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Unless  otherwise  mentioned  in  tlie  text,  the  symbols  have  the  following 
meanings : 


a 


C 

H 

f 

F 

h 

h' 


h^.hz 


i 

^il 


a wave  number,  describes  the  time  lag  of  propagation  of  a wave 
in  the  y-direction 

'a'  for  the  steadily  moving,  constant  amplitude,  wave 
spatial  exponent  of  growth  of  steadily  moving  temperature  wave 
wave  velocity 
specific  heat 

Young's  modulus  of  elasticity 
face  load 
non-dimensional  f 
mean  thickness  of  lubricant  film 
small  surface  perturbation 

amplitudes  of  sine  and  cosine  components  of  h' 

operating  surface  waviness 

film  shape 

small  change  in  h 

initial  surface  waviness 

non-dimensional  h 

non-dimensional  h 

non-dimensional  h^ 

thermal  conductivity  of  metal 

thermal  co  iductlvity  of  lubricant 
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k,,  = tlienii.-il  cliff usivity  ol  metal  — 

M ^ pC/ 

k = thermal  cli  f fusivi  ty  of  lubricant 

oil  ^ 

X = circumference  of  seal 

L = width  of  seal  face 

p = pressure 

p ^ = zero-order  pressure 

p"  = pressure  perturbation 

p'  = zero-average  component  of  p" 

p = constant  component  of  p" 

P = non-dimensional  p 

o o 

P"  = non-dimensional  p" 

P'  = non-dimensional  p' 

^ A 

P = non-dimensional  p 

q = heating 

Q = volume  rate  of  flow 

t = t Ime 

T = temperature  perturbation 

1 = amplitude  of  T 

o 

u = velocity  in  x-direc t ion ; also  non-dimensional  sliding  speed 

u'  =»  perturbation  velocity  in  x-direction 

U = sliding  speed 

V = surface  displacement 

w * velocity  in  z-direction 


width  factor  for  seal  face  ' 
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» distance  along  the  edge  of  contact 
* non-dimensional  x 
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= distance  normal  to  x 
= non-dimons  iona  1 li ' 

= distance  across  tlie  width  of  the  face 

= coefficient  of  thermal  expansion  of  metal 

= t ime -exponent  of  )^rowth  of  T 

= surface  deflection 

= thermal  deformation 

= elastic  deformation 

= non-dimensional  6 , 

th 

= non-dimensional  6,, 

r, 

= wave  number  (2TT/wave  length) 

= non-dimensional  h 

o 

= (ifx  + ay  - >Cct ) 

= (<x  - <ct) 

= viscosity 
= Poisson's  ratio 

= thermoelastic  displacement  potential 
= density 
= normal  stress 
= shear  stress 

= spatial  exponent  of  growth  of  temperature  wave  in  y-direction 
* non-dimensional  location  of  a load  element 
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CHAPTKK  I 


INTRDDI'CTION 


In  recent  years,  it  has  been  rc'ccsnized  that  liigli  speed  sliding 
contact  is  associated  with  macroscopic  instabilities  of  interfacial 
pressvire  and  temperature.  Generally  referred  to  as  "thermoelastic 
instability",  the  phenomenon  occurs  through  the  interaction  of  frictional 
heating,  lieat  transfer,  interfacial  pressure  generation  and  wear.  At 
certain  operating  conditions  and  for  certain  combinations  of  material 
properties,  thermoelastic  instability  can  cause  nominally  flat  surfaces 
to  deform  locally  so  that  the  load  acting  on  the  interface  is  carried 
mainly  by  a number  of  patches  or  asperities.  Tliis  implies  high  contact 
pressure  at  certain  points,  with  surrounding  regions  of  lower  pressure 
and  relative  surface  parting  which,  in  turn,  concentrates  the  frictional 
heating  at  these  points.  The  resulting  high  temperature  at  these  points 
causes  them  to  expand  rapidly,  thereby  aggravating  the  situation  through 
a progressively  worsening  cycle  of  reduced  "contact  area",  more  concen- 
trated heating,  greater  thermal  expansion.  The  surface  configuration, 
once  disturbed,  can,  at  least  in  theory,  undergo  runaway  and  catastrophic 
deformation,  and  hence  the  term  "instability".  In  practice,  of  course, 
a damper  is  expected  to  be  provided  by  the  mechanism  of  wear  and/or 
elastic  suppression  of  the  thermally  induced  growth.  Wear,  which  in 
turn  di'pends  on  interfacial  pressure,  can  wear  down  a prospective  asperity 
faster  than  it  grows;  normal  pre'’sure  on  the  surface  can  elastically 
t latten  or  squash  down  an  asperity  at  the  mt>ment  of  its  conception. 
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lliiluT  of  thfsi’  c i rcums  t .Tiici’s  rcprt'seiii  s a s Lali  i 1 i z i ii>»  loc'dl)ack  arranni’- 
im’nt  wluTobv  any  prospect  of  asperity  formation  is  nipped  in  tlie  bud. 
However,  since  botti  of  tliese  stabili:'ing  mechanisms  depend  on  tile  mag- 
nitude of  interfacial  pressure,  it  is  conceivable  tliat  asperity  forma- 
tion may  proceed  to  a certain  point  tie  fore  the  pressure  is  lar^e  enouKli 
to  ni‘f;ate  further  thermal  urowtli.  In  sucli  a case,  tliermoe  last  ic  insta- 
ll i 1 i t y will  still  be  said  to  occur,  although  ,'i  stable  patcli-like  contact 
configuration  will  eventually  evolve.  Such  patch- like  contact,  in 
fact,  mav  hi  incompatible  with  many  applications. 

Thermoelastic  instability  can  conceivably  affect  numerous  practical 
applications  involving  differing  geometrical  configurations.  It  was 
experimentally  observed  by  Sibley  and  Allen  lIJ  while  studying  the 
friction  and  wear  characteristics  of  unlubricated  ceramics.  Three 
stationary  pucks  were  pressed  against  the  face  of  a rotating  washer,  as 
shown  in  fig.  (1).  This  configuration  approximates  that  of  disc  brakes. 
For  certain  materials  and  sliding  speeds,  bright  hot  spots  were  observed 
to  form  and  move  in  apparently  periodic  patterns.  It  was  concluded, 
further,  that  higher  friction  coefficients  (producing  high  frictional 
heating)  and  lower  wear  rates  were  conducive  to  more  dramatic  displays 
of  hot  spots. 

Much  the  same  conclusions  were  drawn  by  Barber  l2j  who  analyzed  the 
case  of  three  pins  pressed  against  a moving  surface,  as  illustrated  in 
tig.  (2).  Arguing  that  all  contact  problems  consist  of  surfaces  touching 
only  at  a number  of  points  (as  opposed  to  geometrically  true  surface 
contact).  Barber  went  on  to  show  that  even  slight  disparities  in  the 
loads  carried  by  the  pins  (as  is  almost  Inevitable  in  practice)  leads  to 
the  highest  frictional  force  at  the  end  of  the  pin  bearing  the  greatest 
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I lo.id.  riirtt  pin,  tlien,  is  lifati'd  tin*  most,  expands  the  most  and  carries 

an  incieasin^ly  greater  sliare  of  the  total  load.  Kventuallv,  this  process 
can  lead  to  all  the  load  l>ein)^  carried  by  tlie  sinj'le  pin,  with  the  other 
two  liftinji  oft  tile  surface.  If,  however,  the  first  pin  acquires  a 
share  of  tlu-  load  lor  which  its  wear  rate  exceeds  the  rate  of  thermal 
expansion,  it  will  begin  to  be  worn  down.  A second  pin  will  begin  to 
carry  more  and  more  load  while  the  first  will  cool  off  and  lift  from 
the  moving  surface.  Kventually  it  will  be  the  third  pin's  turn,  and 
thus  the  load  will  continue  to  shift  from  one  pin  to  another.  If,  on 
the  other  tiand,  tlie  rate  of  wear  Is  greater  than  that  of  thermal  expansion 
from  tfie  very  outset,  all  three  pins  wil'  establish  simultaneous  contact 
and  a stable  situation  will  prevail.  It  is  also  possible  to  have  a 

[ quasi-stable  contact  configuration  consisting  of  a single  pin  in  con- 
tact, with  a rate  of  frictional  heat  generation  identical  to  that  of 

the  heat  lost  through  conduction,  convection  and  radiation. 

The  above  findings  were  reinforced  by  Dow  L^J,  who  analyzed  a 
similar  model  with  two  pins.  The  pins  were  assumed  to  he  so  short  that 
their  temperatures  were  functions  of  time  alone,  not  of  space.  For 
given  material  properties,  the  criterion  for  the  onset  of  thermoelastic 
instability  was  determined  to  be  a critical  sliding  speed.  It  was  shown, 
further,  that  such  a critical  sliding  speed  is  obtainable  irrespective 
of  whetlier  wear  is  present  or  absent,  and  that,  in  both  cases,  sub- 
critical  speeds  correspond  to  thermoe lastical ly  stable  operation. 

Although  the  pin  model  is  realistic  in  that  so-called  surface  con- 
tact actually  consists  of  multi-point  contact,  it  falls  to  account  for 
the  conduction  of  heat  along  the  surface,  from  one  point  of  contact  to 
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ailotluT.  Itiis  effect  was  considered  by  [)ow  and  burton  who  analyzed 

the  two  dimensional  model  of  a semi -inf initc  knife-edge  or  scraper  in 
sliding  contact  with  the  surface  of  a semi-infinite  half-space  as  illus- 
trated in  fig.  (3).  Wear  was  neglected,  restricting  the  applicability 
of  the  analysis  to  boundary- lubr icated  systems  s\icli  as  the  apex  seals  in 
a Wankel  engine.  Once  again  a critical  sliding  speed  was  determined, 
above  which  instability  would  set  in.  In  a subsequent  paper  L5j,  the 
authors  studied  the  role  of  wear  on  the  same  contact  geometry  and  found 
that  thermoelastic  instability  would  occur  beyond  the  critical  sliding 
speed  and  that  the  resulting  patches  of  contact  would  then  either  move 
along  the  surface  or  oscillate  in  stress  intensity  about  fixed  nodal 
points . 

Mechanical  seals  constitute  a class  of  widespread  applications  in 
wtiich  thermoelastic  instability  can  be  of  crucial  importance  in  that 
surface  distortions  can  jeopardize  the  sealing  effect.  A substantial 
body  of  literature  exists  for  thermoelastic  instability  in  face  seals, 
wherein  the  seal  is  represented  as  a pair  of  thin-walled,  co-axial  tubes 
meeting  end  to  end,  with  tangential  sliding  at  the  interface  (fig.  4). 

The  role  of  wear  has  been  reconfirmed  by  Heckmann  and  Burton  i.6j  who 
have  shown  that,  other  variables  being  identical,  the  system  with  the 
lower  rate  of  wear  will  have  the  lower  critical  sliding  speed  and  that 
if  tlie  wear  rate  exceeds  a certain  critical  value,  instability  will  not 
occur  at  any  finite  speed.  Burton  et  al  [?]  have  established  that  wfien 
a good  thermal  conductor  slides  on  its  own  kind,  unrealistically  large 
values  of  friction  coefficient  are  required  for  tlie  occurrence  of  thermo- 
elastic instability;  they  have  concluded  that,  for  a pair  of  identical 
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conductiiiK  niiUfrials,  thornioi- last  if  ii>stal)i  1 i ty  will  not  occur.  A 
thermal  conductor  sliding  on  a Kood  insulator  will  experience  such  an 
instability  the  most  readily.  burton  ts]  lias  elaborated  on  this  aspect 
and  pointed  out  that  even  wlien  a pair  of  thermal  conductors  are  used, 
thin  deposits  of  oxides  can  render  one  surface  a poor  conductor  relative 
to  the  other.  A disturbance  may  then  be  relatively  stationary  on  one 
surface,  which  will  behave  like  a Rood  conductor,  wli  i le  moving  rapidly 
on  tile  otfier,  which  will  behave  like  a poor  conductor.  .Such  arguments 
have  led  to  greater  emphasis  on  conductor- insulator  pairs  in  face 
seal  studies. 

In  engineering  practice,  flat  surfaces  are  never  geometrically  flat. 
Nevertheless,  it  has  been  shown  that  even  perfectly  flat  surfaces  will 
develop  asperities  l4j.  In  general,  the  roughness  of  a surface  is  com- 
posed of  roughness  components  of  numerous  different  wavelengths,  and 
the  largest  wavelength  is  known  to  have  tlie  lowest  critical  speed  .5j. 

In  the  case  of  the  radial  face  seal,  the  largest  possible  wavelength  is 
the  circumference  of  the  seal.  It  has  also  been  shown  that  initially 
flat  surfaces  are  deformed  by  thermoelastic  instability  so  that  contact 
occurs  through  a number  of  patches  or  hot  spots  which  ultimately  acquire 
a limiting  configuration  [.9J  in  which  the  size  of  a patch  is  independent 
of  load,  but  dependent  on  sliding  speed.  burton  and  Nerlikar  LlOj  have 
shown  that  the  effect  of  increased  load  is  heightened  thermal  effects 
leading  to  greater  surface  deformation  and  increased  leakage.  They 
have  also  demonstrated  LiiJ  that  if  initial  surface  curvature  is  taken 
into  account  and  the  load  at  the  outset  is  assumed  to  be  borne  by  a 


number  of  Hertzian  contacts,  increase  of  sliding  speed  monotonlcally 
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dfCit‘;i  SI' s till'  width  oi  t lu'  contact  patches.  This  effect  lias  been  shown 
to  he  quite  dramatic  it  the  IliTtzian  patches  are  assumed  to  overlap, 
implvinn  full  contact  hut  non-uniform  pressures.  Kilaparti  i.12]  has 
exper imenta  1 Iv  demonstrated  the rmoe i as t i c instability  bv  running  a metal 
cup  witli  a flat  lip  on  a slab  of  glass.  At  high  speeds,  bright  flashes 
or  hot  spots  have  been  seen  to  form  and  slowly  traverse  the  face.  Photo- 
graphs ot  the  metal  face  taken  after  such  runs  (fig.  5)  show  dark  tracks 
indicative  of  tlie  path  traversed  by  the  hot  spots. 

Iheri’  are  numerous  other  applications  wtiere  thermoc  lastic  instability 
may  be  of  consequence.  In  gas  turbines  '.13j,  where  very  high  speeds 
and  temperatures  prevail,  severe  problems  can  conceivably  arise  from 
occasional  contact  between  blade  tips  and  shrouds.  In  electrical  machinery 
involving  high  current  flows,  frictional  and  electrical  heating  can 
significantly  alter  the  contact  at  the  brushes,  thereby  generating  high 
resistance  in  the  current  flow  paths  and  further  exacerbating  electrical 
heat  generation.  Such  problems  have  in  fact  been  reported  and  now  await 
ana  lys i s . 

All  of  the  above  work  was  concerned  with  dry  friction.  With  wear 
assumed  to  be  negligible,  the  results  could  be  extended  to  include 
boundary  lubricated  cases.  Left  out  was  the  large  body  of  applications 
in  which  fluid  films,  liquid  or  gaseous,  are  maintained  between  the 
sliding  surfaces.  The  present  work  seeks  to  fill  that  gap  by  asking 
whether  thermoelastic  instability  can  occur  in  the  presence  of  a viscous 
lubricating  film. 

In  view  of  the  complexities  involved,  the  study  has  been  conducted 
in  three  phases.  The  first  phase,  simplified  by  an  idealized  model, 
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dial''  with  '^uitaci's  tlial  arc  piTtcctly  flat.  No  <)iic  s t i oils  arc  asked  as 
to  whether  such  flat  surfaces  can  exist  in  practice;  instead,  the  ques- 
tion is  posed  as  follows;  "If,  somehow,  tlie  surfaces  are  flat  and 
separated  by  a film  of  liquid  lubricant,  is  there  a critical  sliding 
speed  abov'e  whicl)  I hi  rmoe  las  t ic  instability  will  occur?"  Such  a 
critical  speed  is,  in  fact,  obtained  and  seen  to  depend  on  the  thickness 
of  tile  film  and  material  properties.  Furthermore,  it  is  seen  that 
instability  can  occur  only  for  a range  of  film  thickness.  The  simplified 
model  permits  a fairly  rigorous  theoretical  analysis. 

1 he  second  phase  consists  of  experiments  run  with  the  purpose  of 
determining  whether  thermoelastic  instability  can  be  produced  in  the 
laboratory.  Not  only  have  they  been  observed  and  recorded,  but  also, 
much  to  the  author's  surprise,  the  measured  critical  speeds  have  been 
found  to  be  embarrassingly  close  to  those  predicted  by  the  idealized 
analysis  of  the  first  phase. 

Finally,  it  is  recognized  that  real  surfaces  are  never  geometrically 
flat.  Kven  surfaces  that  are  almost  flat  at  the  outset  will,  under 
thermoelastic  effects,  evolve  through  a succession  of  stable  states  into 
wavy  surfaces  as  reported,  for  instance,  by  Taniguchi  and  Ettles 
Nonlinear  effects  make  the  analysis  somewhat  difficult.  It  is  found, 
however,  that  in  theory  a number  of  operating  conditions  can  be  defined 
to  represent  neutrally  stable  thermoelastic  states.  In  practice,  such 
states  will  not  be  reached  because  they  lie  beyond  the  critical  speed 
di’fined  In  the  first  phase,  and  this  critical  speed  cannot  be  exceeded. 
This  explains  the  close  conformity  betwi'en  theoretically  predicted  and 
expiT  imenta 1 ly  measured  results,  but  not  the  exact  mode  by  which  the 
critical  speed  is  actually  attained. 


CHAI’H-R  II 


INSTAIill.ITY  FDK  INITIALLY  FIJVT  SlIKFACKS 


2.1  Cleometry  iind  Procedure 

Two  circular,  co-axial  cylinders  are  pressed  end  to  end  and  one  is 
made  to  turn  against  the  other  (fig.  4>.  To  simplify  the  analysis, 
the  cylinders  are  replaced  by  two  flat,  co-planar  blades  that  slide 
relative  to  each  other  along  a straight,  common  edge  or  interface. 

This  is  not  a great  departure  from  the  original  configuration  as  long 
as  the  wall-thickness  of  the  cylinders  is  small  compared  to  the  radius. 
This  approximation  will  permit  the  use  of  the  simpler  cartesian  coordi- 
nate system  in  applying  such  equations  as  the  two-dimensional  Fourier 
heat  conduction  equation. 

During  sliding,  if  the  pressure  is  uniform  along  the  interface, 
the  temperature  will  gradually  rise  to  a steady-state  value  determined 
by  distant  boundary  conditions.  If,  however,  a small  disturbance 
perturbs  the  uniform  pressure,  prevailing  conditions  may  cause  this 
disturbance  to  grow,  to  be  damped  out,  or  to  remain  unchanged.  If  it 
remains  unchanged,  the  system  is  in  neutral  stability.  In  that  case, 
the  effects  triggered  by  the  perturbation  will  together  result  in 
regenerating  the  original  perturbation  via  some  sort  of  a "c losed- loop" 
phenomenon.  Since  the  disturbance  Itself  may  be  expressed  as  a Fourier 
series  of  waves  along  the  common  interface,  the  stability  of  the  pres- 
sure (or  temperature)  distribution  may  be  Investigated  in  terms  of  tlie 
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■•tahilitv  I't  till'  ciimporu'nt  wavfs.  In  pracL  icf , siich  perturbations  are 
often  eausecl  bv  small  particles  of  debris  in  the  interfacial  spaci'. 

It  lias  been  |)ointed  out  that  a like-on-like  mati'rial  pairing  pro- 
vides tiu'  best  stability  of  operation,  whereas  a good  thermal  conductor 
sliding  on  a non-conductor  has  the  greatest  potential  lor  instability. 

Kor  this  reason,  the  latter  kind  of  material  pairing  is  being  studied. 
Ibis  alsi>  leads  to  considerable  simplification  of  the  analysis  and  a 
better  chance  of  interpreting  the  relative  influence  of  various  material 
prope  ft ie  s . 

The  problem  is  formulated  as  a linear  one  with  linear  lieat  transfer, 
thermal  expansion  and  elastic  displacement.  The  materials  are  assumed 
to  obey  Hooke's  law  and  to  have  constant  properties.  The  properties  of 
tlu'  lubricant  are  also  assumed  constant  to  begin  with  and  it  is  hoped 
that  the  results  of  the  simplified  analysis  will  provide  valuable  insigh.t 
into  the  effects  of  variable  properties.  The  problem  is  then  dealt  with 
as  follows  (see  fig.  b): 

1.  It  is  assumed  that  the  two  surfaces  are  nominally  flat  at  the 
outset  and  arc  separated  by  a lubricant  film  of  a certain  thickness 

h. 

2.  A Small  perturbation  is  subsec|uent  ly  introduced  in  the  form 
of  a t<'mperature  wave  which  causes  a small  surface  wave  h'  to 
develop  on  the  conducting  body.  Since  in  the  original  cylinder 
configuration  such  a wave  would  have  been  continuous  around  the 
circumference,  in  the  ideal  iiJcd  flat-blade  model  the  component 
harmonics  of  the  wave  must  complete  one  or  more  full  cycles  between 


the  "ends"  of  the  plate. 


Figure  6 


Block  Diagram  for  Analysis  of  Initially  Flat  Surfaces 
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Tho  sxirfatt'  wave  lx'  altx-rs  the  vx'locity  ilistriixxition,  the 


If) 


1.  Tho  sxirfaoo  wave  h'  altx-rs  tin-  vx-locity  ilistriixxition,  tho 
prossiiro  <li  s tr  ihxit  ion  and  tin-  hoat  tr.ansfcr  in  tho  oil  film.  Tho 
onersy  equation  is  xisod  to  I iinctional  ly  rolato  tho  tomperaturo 
por  txirhat  ion  T with  its  ottoct  h'  (or  vice  versa,  as  far  as  the 
practical  mechanism  is  concerned). 

•'».  Tho  rosultinj;  prossxiro  disturbance  is  determined  from  the 
linearized  Reynold's  equation.  Tliis  then  determines  the  elastic 
(or  pressxxre-Renerated)  sxirface  displacement  Sj,. 

5.  The  perturbed  heat  balance  causes  a perturbation  in  the  thermo- 
elastic displacement.  This  is  obtained  by  solving  the  thermoelastic 
equation  and  is  called  the  thermal  deflection  ® • 

fx.  Finally,  it  is  postulated  that  a neutrally  stable  situation  docs 
exist,  for  which,  as  explained  earlier, 

"’='th^S: 


In  this  equation,  the  unknown  is  the  sliding  speed  U.  If  a solution 
exists,  it  represents,  for  the  given  conditions,  the  critical  sliding 
Speed  which  defines  the  threshold  of  instability. 

Kvidently,  no  attempt  is  being  made  to  define  or  determine  the  entire 
time  dependent  behavior  of  any  perturbation.  The  time-dependent  problem 
is  being  confined  to  the  special  instance  where  the  perturbation  neither 
grows,  nor  decays  with  time.  In  other  words,  this  analysis  does 
not  concern  itself  with  such  questions  as  to  how  the  perturbation  grows, 
if  it  grows,  or  how  it  decays  if  it  decays.  Rather,  it  intends  to  seek 
out  the  threshold  at  which  such  a system  will  become  xinstable.  It  is 
hoped  that  in  this  way  a satisfactory  design  or  operational  criterion 


can  be  arrived  at. 
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2.2  Hquations  of  Motion 

The  upper  plate  consists  of  a good  thermal  conductor,  the  lower 
plate  of  a good  insulator.  At  a certain  instant,  a surface  displacement 
h'  occurs  on  the  top  plate.  This  displacement  represents  a small  pcr- 
turhation,  and  has  a wavelength  that  is  largo  compared  to  h,  the  film 
thickness,  and  |h'|  the  amplitude  of  the  perturbation. 

The  upper  plate  is  moving  steadily  with  a velocity  U as  shown  in 
fig.  (7),  the  lower  plate  being  at  rest. 

For  near-parallel,  viscous,  stream-line  flow,  the  Navier -Stokes 
equations  reduce  to  the  two  dimensional  equation  for  creeping  motion 


^ u dp 

^ “T  “ d7 

assuming  that  the  face  is  wide  enough  to  permit  the  ignoring  of  "end- 
leakage"  effects. 

For  cons  tant  u, , 


I 


i^L. 

11  dx  2 


+ Ay  + B 


(2-1) 


The  velocity  profile  will  change  as  shown  in  fig.  (7).  Moreover, 
the  net  velocity  at  any  point  may  be  thought  of  as  the  superimpos ition 
of  a pure  Couette  flow  on  a flow  between  two  plates  in  the  presence  of 
a pressure  gradient.  For  the  Couette  flow,  fig.  (8), 

2y  -F  h 


"2  “ 


tJ 


(2-2) 


2(h  + h') 

For  the  flow  with  pressure  gradient,  fig,  (9),  the  following 
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Figure  7.  Coordinates  System  and  Composite  Velocity  Profiles 


Figure  8.  Couette  Flow,  Definition  Sketch 
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UsinK  ttic  first  ol  tliosc,  i-q.  (2-1)  ^ivcs 
- B 

2 ^ ax  8 

Bounilarv  comlitian  (ii)  nivfs 

/ li  . , ' . „ I dp  h ^ hh ' li ' 

Wk 

Addinii  tom'tiu'r  oqs.  (2-J)  and  (2-4l,  one  obtains 

^ 2^  dx 

Kq.  (2-3)  then  gives 

„ 1 dp  - (h'  h'' 

wtience , 


1 ^ I ^ U 

U.  = r; V y - h 


hh ' h^\ 


(2-3) 


(2-4) 


1 2, 


Then,  since  the  net  velocity  is  given  by 


U = Uj+  U2 


one  obtains 


fV 
u = — ^ 


2v  + h'  ...  ^ I 2 , , hh ' h ' 


U + 


2 (h  + h ' ) 


2p,  dx 


y - h'y  - 


(2-5) 


in  which  dp/dx  has  yet  to  be  determined. 


2.3  Continuity  Equation 

Initially,  prior  to  the  occurrence  of  the  surface  displacement, 
Couette  flow  was  in  effect,  for  which  the  net  flow  is 

,,  h 

Ql  = PU  7 

After  the  disturbance  h'  has  occurred, 

Q2  “ pU  ^ + pu^  "I  (^  + (refer  to  fig.  (9)) 
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For  contiimi  Cy  of  flow, 


Ql=  Q2 


whence 


u « 
c 


2 h' 
h + h' 


U 


(2-6) 


2.4  Pressure  Gradient 

In  determining  the  pressure  gradient  dp/dx,  the  flow  will  again  be 
broken  up  into  two  parts,  one  with  a linear  velocity  profile,  the  other 
with  a parabolic  velocity  profile,  fig.  (10).  The  shear  stresses 
balance  in  the  Couette  flow  part  (subscript  L denotes  linear)  of  the 
composite  flow.  For  the  pressure  flow,  the  condition  for  equilibrium  is 

2t  -dx  = - ^ 1 h + h'  dx  (2-7) 

p >x  \ 


where 


V" 


\hy  / 


or, 


T ' 

P 


2u 

c 

Vi  Z 

h + h' 
2 


, from  eq.  (2-7), 

dx 


- JZ  ' - V 2 

1.  I U I ** 


[h  + h ' 

Noting  that  h'/h  « 1 and  dropping  terms  containing  squares  of 
h'/h,  one  obtains 


(2-8) 
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Subs  t i tut  in;>,  for  dp/dx  from  oq . (2-8)  into  t-q.  (2-5)  gives  u in 
terms  of  ii  . Use  of  eq.  (2-6)  tlnm  yields  the  velocity  distribution 
in  tile  flow.  Dropping  terms  containing  higher  powers  of  h'/h,  one 
t inds 

u = ■ 3 14  ^ ^ - 


^ > 5 h', , 

^ ^ ^ ^ + _ . _ _ u 

h h h h 


The  perturbation  velocity,  then,  is 

, '3h'  2 h’  5 h'^  „ 

u ' = \ y'z?y"4Z"'  ^ 

h h h ' 


(2-9) 


(2-10) 


2.5  Non-steady,  Moving.  Temperature  Perturbation 

In  the  general  case  of  a moving  temperature  wave  whose  amplitude 
is  time-dependent,  the  temperature  distribution  in  the  metal  may  be  of 
the  form 

T(x,y,c)  = Atl(t)Y(y)sin(<x  + ay  - iCct) 
and  must  satisfy  the  two-dimensional  heat  conduction  equation, 

1 

dx^  dy^  Si 

Further,  T will  be  required  to  satisfy  the  boundary  condition  that 
at  the  metal -oil  interface 

3 1 

T s T = T e^  sinC'.x  - Kct) 
s o 

where  T is  a constant  amplitude  of  the  wave.  It  is  easily  seen  that  a 
o 

temperature  distribution  T(x,y,t)  that  satisfies  tlie  above  conditions  is 
T = T e^*'e  ^^sin(''',x  + ay  - ’fet)  (2-11) 


where 


<c 

2V 


(2-12) 


T 


2-3 


(2-r3) 


Kquations  (2-12)  and  (2-13),  wlicn  solved  for  ? and  a in  terras  of 
C , give 


and 


2' 

+ V 


i ' /'f, 


2 2 
/ c 


1/2 


a = ± 


2'' 

+ K 


+ 


2 2-1/2 
c 


(2-14) 


(2-15) 


A few  worils  are  in  order  here  on  the  nature  of  the  parameters 
a and  c.  It  is  evident  fromeq.  (2-11)  that  B represents  the  time-exponent 
of  growth  of  the  temperature  perturbation  T.  The  amplitiidc  of  this  wave 
grows  or  diminishes  accordingly  as  6 is  positive  or  negative  and  remains 
constant  if  B = 0. 

The  quantity  | is  the  spatial  exponent  of  growth  of  T in  the  "y" 
direction,  i.e.  , into  the  body.  It  may  be  noted  that  in  writing 
eq.  (2-11),  no  explicit  boundary  condition  was  imposed,  nor  was  a specific 
set  of  axes  of  co-ordinates  defined.  An  implicit  boundary  condition  is 
imposed  by  oneq,  (2-11).  Since  the  perturbation  T is  dissipated  as  it 
travels  into  the  body  from  its  point  of  application,  it  is  obvious  that 
as  y -♦  *,  T -*  0.  This  requires  that,  ineq.  (2-11),  F must  always  be 
positive.  Ineq.  (2-14),  then,  the  outer  minus  sign  must  be  ignored,  as 
must  the  inner  minus  sign  in  order  to  ensure  real  values  of  f. 


The  quantity  c is  the  traversal  velocity  of  the  temperature  wave  T. 
Equation  (2-11)  is  written  for  a right-running  wave;  c may  be  either 
positive  or  negative  to  represent,  respectively,  a right-running  or  a 


left -running  wave. 
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As  till'  piTturbat  ion  T penetrates  tlu'  i>o<ly  in  the  "y"  direction,  its 
intensity  is  d imin i siicd , At  the  same  time,  however,  at  its  surface  of 
application  (y=0,  say),  the  original  wave  has  mcived  a distance  of  'c.t'. 

Thus  at  any  plane-  y > 0,  the  corresponding  perturbation  q 

only  smaller  than  T in  amplitude,  but  also  lags  behind  T in  the  direction 

of  motion  ol  the  wave  (fig.  11).  Tlte  quantity  a is  a measure  of  this 

lag,  Kquation  (2-12)  indicates  that  c/a  must  have  the  same  sign  as  ? and 

since  imis  t always  be  positive,  c/a  is  always  positive.  Thus  c and  a 

must  alwavs  have  the  same  sign.  Figure  (11)  is  drawn  directly  from  eq.  (2-11  i, 

It  is  to  be  noted  that  since  the  lag  is  a result  of  the  propagation  of 

T,  a non-moving  wave  (c=0)  implies  a=0. 


l.imiting  Cases 

One  limiting  case  is  that  of  c=0,  viz.  the  non-moving  wave.  In 
this  case,  as  is  evident  from  eqs,  (2-14)  and  (2-15), 


a = 0 


so  that 

/p  2 

-y  — + K 

Pit  M 

T = T c e sin  vex 

o 

Another  limiting  case  is  that  of  p=0,  viz,  the  case  of  the  steadily 
moving,  constant  amplitude  wave.  In  this  case 


, r 1 2.  VC  / 2 C l 

I 2 2 ' . 2_j 


(2-lb) 
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Ttu-  tomporaliiri'  ti  is  Lr  ihut  icni  is  f'ivfii  by 

T = T c ''^sinCt^x  + a.y  - <ct) 
o 1 

so  that  on  tiu-  motal  surface  (y=())  , 

T = T sin/i^x  - /ct)  T sin  •, 
o ^ o 1 


(2-17) 


(2-18) 


(2-19) 


2.6  Temperature  Distribution  in  Oil  Film 

The  temperature  profile  in  ttie  oil  film  is  approximated  by  a poly- 
iKwial  of  tile  second  degree  (see  fig.  (12)  and  Appendix  A): 

T = A'y^  + B'y  + C (2-20) 

in  which  A',  B' , and  C are  functions  of  x and  which  conforms  to  the 
boundary  conditions: 

at  y = + h 

h , 

at  y = •2-  + h 

h 

at  y = - 7 

The  first  two  boundary  conditions  derive  from  the  fact  that  at  the 
oil -metal  interface  the  temperature  is  the  same  for  the  metal  and  the 
oil  and  tlie  heat  conducted  in  each  medium  is  the  same.  Boundary  con- 
dition (2-23)  refers  to  the  fact  that  the  second  body  is  a thermal 
insulator. 

From  eqs.  (2-20)  and  (2-23),  one  obtains 
B - Ah 


ot  . 

T = T c sin  X, 
o 1 


T'T  Si  9t._  . , ^ . 

— = - re  (-  sin  a cos  X,) 

>y  K , , o - 1 1 

Oil 


0 


(2-21) 

(2-22) 

(2-23) 


Figure  12.  Temperature  Distribution  in  Oil  Film 
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Kr('m  k-qs.  (2-20)  and  (2-22), 
K 


A = - 


M 


2K  o 

o 1 1 


pt/ 

T e ( ^ s in  > , - a cos  > , — > 


1 ^ h’' 


h 


Kromfqs.  (2-20)  and  (2-21), 


° t D — 2 - . 

C = T C sin  ) ,-  B 7-  h + 2hh 
o 14 

so  that,  iinally,  corrected  to  the  first  order  of  the  pi-rturhation , 


r = T e s in  ^ - Trr T e ( F s in  > - a cos  > — i hy  + y - -r-  h (2-24 ) 

o 1 2K  . , o \ ' 1 1 / — 4 

oil  h 


As  before,  two  limiting  cases  may  exist.  For  c=0,  the  non-moving  wave. 


Pt 


Si  „ ot  le  2 

T e k — + V 


T = Te'  sinKx  - T e ' — + v sin  Fx 

" 2K  .,h  " 'Si 

Oil 


2 3 

n,y  + y - - h . 


For  the  steadily  moving,  constant  amplitude  wave,  °=0,  eq.  (2-24) 


gives 


T = T sin  ^ 9 _ (hy  + y^-  | h^j  (b  sin  a^cos  (2-25) 


2K  .,h 
oil 


2.7  The  Knergy  Kquation 

Under  the  stated  thin-film  assumptions,  the  energy  equation  takes 


the  form 


l"  6y 

y\ 


dy 


[ — I dy 
^ Piy/ 


(2-26) 


Here  C is  the  specific  heat  of  the  oil.  The  first  and  second  term 
on  the  left  hand  side  are,  respectively,  the  convective  and  local  con- 
tributions to  the  increase  of  internal  energy  of  the  oil.  The  third  term 
on  the  left  is  the  negative  of  the  heat  conducted  into  the  oil  while  the 
right  hand  term  is  the  heat  generated  through  friction. 
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Dof ine 


'’c/  j P ^ If 
^1 

^1 


K 


2 

^ T , 

oil 

^y 


f2  f^u'  ^ 

■1 

Using  cqs,  (2-2A)  and  (2-9),  one  finds 


U'^T  e 


?t 


'’cc  2k 


• K . ,li  cos  >.,+  7 


1 r2  f 


oil 


"oil"  """  "r  4 V"  I®  'H'*'  “ ^ 


iJj 


For  the  non-propagating  wave,  (c=0) , eq.  (2-27)  reduces  to 
ot 


Ui'T  e 
o 


cc  2k 


oil 


^ — 1—2  Jh  ? 

•^oilh  4 V ■X'"''' 


(2- 


(2 


Returning  toeq.  (2-27),  the  steadily-moving,  constant  amplitude 
wave,  (P=0) , gives 


Ui'T 


q = rr. ' K . ,h  cos 

^cc  2k  oil 

oi  1 


" ^1^  Z X,+  h cos 


(2 


Equation  (2-24)  yields 
ot 


T c 
o 


^c?  k 


oil 


ot 


''oil-"  I ■" 

oil 


27) 


28) 


-29) 


(2-30) 


JO 


■’or  tiu'  non-propaRatiriK  wave  (c=0)  , this  Oecomos, 
Rt 


T Re 
o 

IT— 

oi  1 


f-r  \ r2  I?  2'-  , 

I li  K , .+  T 7—  + '•'  ,sin 


oil  J 


Si 


/■ 


(2-il) 


The  steaiiily  moviriR,  constant  amplitude  wave  (?=0)  kivcs,  from 

efj . ( 2-  10  ) , 

T 


or 


01 1 


/-  1 — 2 ^ 

1 li  K . ,+  — K^h  b , cos  > ^ 


oil  3 


(2-32) 


The  frictional  clissipation  is  obtained  from  eq.(2-9)  for  the 
velocity  distribution: 


I 


fA 


2 

U_ 

_2 

h 


h' 


Nothing  has  been  said  yet  about  h'  other  than  that  it  is  a small 
surface  deformation  that  corresponds  with  the  temperature  perturbation 
T.  Obviously,  h'  and  T cannot  both  be  arbitrarily  chosen,  as  one  is  a 
function  of  the  other.  In  this  case  an  arbitrary  T has  been  chosen  and 
h'  will  follow  from  it.  At  this  stage,  however,  it  is  appropriate  to 
define  the  form  h'  will  take,  it  being  known  that  it  must  take  some 
wave  form. 

Let  h'  = hj^sin  >.^+  h2Cos  where  h^  and  h2  are  unknowns  to  be  sub- 
sequently determined. 

Then , 


n 

g 


- u. 


li^ 

h 


h.sin 

V 1 


> j+  h2Cos 


(2-33) 


The  above,  derived  from  the  velocity  distribution  alone,  Is  dependent 
on  the  nature  of  the  temperature  perturbation  only  through  h^  and  h2. 
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Tims  tiir  basic  form,  Riven  bycq.  (2-33),  is  true  for  all  cases.  The 
heat  conducted  into  the  oil,  is  derived  fromeq.  (2-24)  and  found  to 


be 


■at 


V ■ Vo‘' 


in  > a cos  >. 


1 j 


(2-34) 


For  c=(),  the  atiove  becomes 


at  /a  2 

’k-  - Vo'  ■ r + ^ 

M 

Kor  '2=0,  oq,  (2-34)  givc.s 

V ■ Vo  ’ ^ r ''1. 

In  all  of  the  above  equations,  constant  fluid  properties  have  been 

assumed.  If  q , q . . q,  > and  q are  substituted  for  in  eq.  (2-26),  tiie 
cc  ^cH  'k  'r 

sine  and  cosine  terms  on  the  left-hand  side  must  be  equal,  respectively, 
to  th(kse  on  the  right-hand  side.  Thus,  using  eqs.  (2-26),  (2-27),  (2-30), 


("2-33)  and  (2-34),  one  obtains; 

[i  i V 0 - 1 V 


h^  e^’'^ 


- 


— r~ 

uirk 


oil 


h^ 


at 


(2-36) 


and 


'r  't„h+  i V 0 - (^lT i V 0 

“■J  bn  I 


—2 

- j K^h  pa  , + ^ 

p,U 


(2-37) 
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Kor  c=0,  i.e.,  for  a non-movin>>  wave,  the  above  two  equations  become: 

-2 


h 't  e"^  , 
h|=  - 


( 1 _ to  2'  ■'  h^  Qt  /a 

^ — Vo^' 


oi  1 


Si 


(2-38) 


and 


r*,. 

f i-’ 

o Ui^ 


uU  k 


oil 


— 7"  ' ^ii+  { V' 


(2-i9) 


Tlie  corrc‘spond  in>’,  relations  for  the  steadily  moviny,,  constant 
.impli  tilde  wave  (==())  are 

h’^To  - I 1 _ - •};2 

'‘r  ■ “X — 1 8 V '"r  3 V - -H  Vo'’ 

uU  k^.j 


o 


’’l'  ■ 2 k 

;lU  01  1 


1 n „ 1 \ 


i-  U - 3 cy  + b_ 


(2-40) 


and 


—9 

h‘‘T  — 

V,  r f tiu-n 


- -r  [■  ( "oi,"  7 V V - (^11+  7 V 


oil 


oil 


>J 


•Si^ 


1j 


(2-41) 


2.8  Thermal  Deflection 

Under  the  action  of  the  temperature  wave  (eq.  (2-11))  generated  In 
the  metal,  its  surface  undergoes  a thermal  deflection  which  may  be 
determined  by  solving  the  thermoelastic  equation: 

_ 2 2 

^ = (1  + v)Qrn^T  (2-42) 

>iy 

where  T is  given  byeq,  (2-11)  and  y is  measured  into  the  metal  from  its 


surface. 


n 


Kquatlon  (2-42)  has  tlie  homogeneous  solution 

c'  *"0  cos  ■>  j+  I)  sin  >,  j'y 

and  a particular  solution 

• = P e e cos  > + Q e'  c •'sin  > 

P 

where  C,  D,  P and  Q are  constants.  In  order  for  to  satisfy  eo  (2-42) 

p - y 

one  must  have 


and 


2 2 2 
K c + ? 


(2-43) 


Q = 


22  2 

/ c + o 


(2-44) 


Then 


■;;  = y + V 
H ‘P 


(2-45) 


whence 

V = cos  X^+  D sin 

^ y=0  ^ 

8t/ 

- i^e  (^P  cos  >.  Q sin  ^ 

+ ae°*'^Q  cos  > j-  P sin 

Here  "v"  is  the  normal  displacement  of  the  surface.  If  this  displacement 
is  suppressed,  a thermoelastic  normal  stress  and  zero  shear  stress 
will  appear  at  the  boundary.  Putting  v^O  in  the  above  equation,  the  sums 
of  the  coefficients  of  the  sine  and  cosine  terms  must  each  be  zero. 


This  leads  to 


}4 


and 


Thus 


C = - (a  Q - 'P) 


= - 7 (f  Q + aP) 


Y = e°''o  C cos  > D sin  ) + e *"0  ^^(P  cos  > + Q sin  X ) 


(?-46j 


(2-47; 


(2-48) 


in  which  C,  D,  P and  Q arc  given  by  eqs.  (2-46),  (2-47),  (2-43)  and 
(2-44)  respectively. 

The  thermal  stress  corresponding  to  the  displacement  potential  V 

in  on.  (2-48)  can  now  be  derived  from  Hooke's  Law; 

2 2 

K ^ V V 

a = '—y  j + V — j - (1  + v)a  'r 

The  surface  stress  is  given  by 


"^s  '^y'y=0 

and  a corresponding  imaginary  pressure  p^^  can  be  defined  such  that 


Pf  " - ‘^s 


which  holds  the  surface  flat,  but  stressed.  In  order  to  determine  the 

thermally  induced  deflection,  under  zero  a , p must  be  relieved  by  the 

s 1. 

super Impos i t ion  of  another  imaginary  pressure  p^  such  that 

p + p =0 

The  surface  can  then  be  thought  of  as  being  stress-free  and  svibjected  to 
the  fullest  deflection.  Dow  [3]  has  shown  that  P2  gives  a deflection; 

2P2  2p^ 

^th“  ~ “ ■ ‘eIT 


(see  Appendix  B,  eq.  (B-5)) 
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('arrvinR  out  the  above  steps,  one  obtains 
2^  aTc- 


/ c + 9 


^ / / 


+ • 1 - T + — : — ■ cos 


(2-49) 


This  gives,  for  the  case  of  tiic  non-propagating  wave  (c=0), 


otj-K-T  e 
M M o 


et 


= 

' th  9 


f. 


Ip  2\ 


sin 


(2-50) 


If  the  temperature  wave  is  moving  steadily  with  a constant  anTjilitude 
(9=0),  cq,  (2-49)  gives 


--a,sin'>,+  (/-b)cos>,- 
th  Xc  1 1 ^ L 


(2-51 ) 


2.9  Elastic  Deflection 


Reynold's  equation  of  hydrodynamic  lubrication  gives 


in  which 


h = h + h' 


(2-52) 


and,  since 


— « 1 , 


,3  ^ 4.  a I'*' 

h * II  I 1 + 3 — , 

h ' 

in  which  h is  constant.  Then,  integrating  Reynold's  equation  twice  and 
dropping  tlie  constants  of  integration,  one  has 

bu.U 


\ 

p = i h2sin  Xj-  h^cos  \^j 


(2-53) 


5f) 


Ihf  const  aiiLs  (>t  i nt  I'l^rat  ion  are  dropped  because  o(  the  need  for  p to 
be  single-valued  around  the  c i rcumlerence , and  because  a pure  constant, 
added  to  p,  will  merely  represent  a uniform  pressure  that  produces  no 
surlace  curvature.  It  has  been  shown  that  (eq.  (H-5)) 


l.quations  (2-5  5)  and  (2-54  > then  give  the  elastic  displacement 
%=  (iSSin  > hjcos 


f2-54) 


(2-55) 


2.10  The  Narrow-Lip  Seal 

In  the  preceding  pastes,  no  account  has  been  taken  of  the  width  of  the 
seal  face  (lip),  although  pressure  and  velocity  have  been  assumed  constant 
across  the  lip.  The  preceding  analysis,  therefore,  applies  to  seals 
with  lips  that  are  wide  enough  to  permit  the  ignoring  of  "end  leakage." 

As  with  "long"  and  "short"  bearings,  we  can  here  treat  the  "narrow-lip 
seal"  as  a separate  class  of  face-type  seals.  Indeed,  in  the  class  of 
seals  studied  here,  narrow  lips  are  of  the  greater  likelihood. 

The  simplified  Navier-Stokes  equations  are 

2 


2 

^ w 
a ^ 


and  the  continuity  equation  is 


(2-56) 


^u  ^w 


0 


(2-57) 


J7 


With  the  present  co-ordinate  system,  fig. 
(lit  ions  are 


(15), the  boundary  con- 


w = 0 

w = 0 


P 


P 


max . 


The  above  equations  may  be  solved  first  for  u and  w and  then  for  p,  to 
obtain 


,,  = i_  iE  ^ y2  _ + 


111'' 

h ' 


w 


p 


3^,U  dh'  ( 2 

-3  dx  ' ^ ■ 4 ’ 
n 


(2-58) 


where  it  has  been  assumed  that 

l£  lE 

and  where  second  and  higher  order  terms  have  been  dropped. 
Along  the  center  line  (z=0) , the  pressure  is  maximum, 


. 3uUL  dh' 


max. 


dx 


(2-59) 


4h 


while  the  average  pressure  is  (fig.  14) 


- 2 

P 


T P 

3 '^max. 


uUb  dh' 

37- 
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111  till'  1 (' 1 1 owi Tif; , p is  treated  as  tlie  working  pressure.  The  elastic 
deflection  is  tlion  given  by 


whence 


= 12 
“K  K< 


» = _ 

"K 


.tn.  dh ' 

3 dx 


If  h'  - hjSin  '‘^+  h^cos  ^ , 


^j..= = 3 i hjco.s  /j-  h2Sin  >,^ 


(2- 60) 


The  temperature  distribution  is  found  to  remain  unchanged  and  is 
still  given  by  eq . (2-24),  it  is  seen  that  q,  , q , q , and  q all 

-iR’  ^g’  ^cc’  ^CJt 

remain  as  before  and  are  given  by  eqs.  (2-33),  (2-32),  (2-27)  and  (2-29) 

respectively.  This  means  that,  in  the  energy  equation,  all  terms 

remain  the  same  and  h^  and  h2  are  still  given  by  eqs.  (2-35)  and  (2-36). 

Thus  the  only  effective  difference  that  the  narrow- lip  consideration 

2 2 

makes  is  in  the  elastic  deflection  term  which  is  now  L“<‘'/12  times  the 
value  given  by  eq.  (2-55).  This  is  a non-dimensional  multiplier  whicii 
may  be  introduced  in  eq.  (2-55)  to  denote  the  departure  of  the  configura- 
tion from  the  wide-lip  one.  Let  us  call  this  factor  W.  Then 


(),.  = — l!(h  sin  X.  - h cos  \ ) 
K ^^2j^3  2 111 


(2-61) 


Here,  for  a narrow-face  seal,  W wi 1 1 have  the  appropriate  value.  The 
wide-lip  seal  is  a limiting  case  of  the  narrow-lip  one  and  is  represented 


by  W = 1 . 
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.'.11  Jii.siahi  Lity  riirosliold  t)t  tlio.  Steadily  Moving  Perturbation 

In  tlif  absence  of  adequate  vinderstanding  of  the  phenotnona  that  occur 
in  face-type  seals,  no  prediction  can  be  made  as  to  whether  a perturba- 
tion will  necessarily  bo  moving  or  stationary.  This  is  why  the  tempera- 
ture pe  rturbat  ion  (eq . ( 2- 1 1 ) ) has  been  defined  to  include,  both  possibilities. 
However,  since  the  intention  of  this  study  is  to  determine  the  conditions 
under  which  the  perturbation  will  neither  grow  nor  decay,  we  shall  require 
that  ?=0.  Therefore  we  shall  look  at  two  distinct  cases,  one  in  which 
the  perturbation  is  non-propagating  (c=0)  and  another  in  which  it  is 
moving  (c-'O).  Although  in  the  preceding  pages  these  two  cases  have  been 
referred  to  as  limiting  cases,  they  arc,  in  fact,  the  only  cases  one  needs 
to  investigate. 

As  has  been  explained  in  section  (2.1),  the  threshold  of  instabil- 
ity is  determined  by  that  combination  of  conditions  for  which 


Remembering  that 


h'  hjSin  >, h2COs 


and  using  eqs.  (2-51)  and  (2-bl),  one  has 

^mVo 


h-S in  ) ,+  h^cos  X 
I 1 Z 1 Kc 


-n.sin  1',+  (K  - b)cos 
11  U 


12  I W ( \ 

+ U I h^s  in  X - h^cos  X ^ J 

I '.w:  h^ 

Til  is  equation  can  be  broken  up  into  two  equations  in  the  two 
orthogonal  sine  and  cosine  components  as  follows: 


■1 


wc 


Uh, 


VX  h 


(2-62) 
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k T 1 ■>  1 1 


and  h^  art'  Riven  hy  eqs.  (2-40)  and  (2-41)  respectively, 
for  hj  and  , one  obtains,  irom  eq.  (2-b2), 


Siil»s  t i tut  inR 


2/“'*'m'Si  - b , 3 *Si' 


3 V'’  ^ ^ { i - ' 

+ o '1,  '■  ^ ; K . ,+  7 K hb 

2 h i'"  k , , 1 / 2m  k , , o I 1 4 m 

oil  oil 


4k  W h 
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12k.  W -3 


7: — ; a , c + — s — b ' 

Ki--  k 1 ,^27 

oil  h h 


(k  . -T  k tib  c 

jj,k.,\oil  3 >!  j 


Similarly,  oq.  (2-62)  gives 


(2-64) 


2'’^^mSi  ‘*1  6W  ' 1 — " r 12W  ! 1 — 

^ ^ ~ ly  k 4 V’^y/  ^ ■ E<  k 3 

oil  oil 


M 2r  Sii  k 3p,  k . ■ p ^ 


(2-65) 


Equations  (2-64)  and  (2-65)  represent  two  algebraic  equations  in  U and  c. 
Owing  to  tbe  complicated  nature  of  these  equations,  it  is  not  feasible 
to  solve  them  simultaneously.  However,  each  can  be  solved  for  a scries 
of  values  of  c,  both  positive  and  negative,  to  provide  corresponding 
values  of  U,  A simultaneous  solution  is  represented  by  that  value  of  c 
for  which  both  equations  give  the  same  U.  These  solutions  are  carried 

out  in  the  computer  and  graphed  in  fig. (15)  in  which  the  chain-dotted 

-4 

lines  indicate  the  area  of  practical  interest  and  W = 10  , which,  for 
a seal  diameter  of  2 inches .corresponds  to  a seal  lip  width  of  0.014 


inch.  The  materials  are  aluminum  sliding  on  glass,  with  K = 1.0  in. 


for  one  complete  wave.  Each  of  the  solutions  is  discussed  separately. 
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Figure  15.  Llne.s  of  Neutral  Stability,  Numerical  Solution  (W  = 10  ) 
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Ttu'  roj;i(in  rcprosontod  by  fiy,.  (15)  may  bo  broadly  subdivided  into 

■/ones  where  one  or  more  of  certain  effects  may  be  predominant.  There 

will,  for  instance,  bo  one  such  /.one  where  convection  effects  are 

ne^ilij’ibly  small.  It  is  apparent  fromaq.  (2-2hi  that  in  this  zone  one 

must  either  ha'’  the  condition 

a)  f q -t  q is  very  smal  I 
cf  ^cc 

or 


b)  0 . = q and  each  is  very  small 
'cl  ^cc 

From  eqs.  (2-24)  and  (2-32)  it  is  seen  that  (b)  will  be  true  when  c and 
h are  both  very  small,  namely  in  the  bottom  left-hand  corner  of  fig.  (15). 
Condition  (a)  requires 


c g U 


and 


J 


which  come  the  closest  to  being  satisfied  when  h and  U have  the  highest 

values.  In  fig.  115),  this  will  be  the  area  in  the  upper  right  hand  corner. 

Pressure  effects  will  be  minimal  when,  as  seen  from  eq.  (2-61),  either 

W is  very  small,  or  h is  very  large  or  U is  very  small.  Tlie  bottom  right- 

hand  corner  of  fig.  (15)  represents  a region  where  the  last  two  conditions 

are  satisfied  and  here,  for  any  reasonable  value  of  W,  pressure  effects 

(l.e.,  % ) will  bo  expected  to  be  insignificant.  Pressure  effects  will 
h 

also  be  insignificant  for  the  trivial  case  of  vanishingly  small  h^  and 
h^ , i.e.,  for  h ^ 0 and  U » 0.  This  would  be  represented  by  the  ".pper 
left-hand  corner  of  fig,  (15). 

It  is  found  that  the  solutions  graphed  in  fig. (15)  do  fairly  fulfill 


the  above  expectations  and  it  is  possible  to  approximately  represent 
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each  family  of  H<’hiti(ins  by  simple  alRcbraic  relations  involving  the 
critical  sliding  speed  Ucri  t . > corresponding  wave  speed  and 

h,  as  follcAi/s 


Family  I (jig.  16) 


U . 
c r 1 r . 
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(Pressure  effects  are  negligible) 


Family  II  (fig.  16) 
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Family  III  (subdivided  into  parts,  fig.  17) 
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Figure-  16.  Lines  of  Neutral  Stability,  Families  1 and  II 
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Figure  17.  Lines  of  Neutral  Stability,  Families 
and  IV  (W  = lo"^) 
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(Kotli  cdnvfction  and  pressure  effects  are  significant) 
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(Pressure  effects  arc  unimportant) 


(2-73) 


Family  IV  (subdivided  into  parts) 
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(Pressure  affects  the  critical  wave  speed  but  has  no  significant  effect 
on  the  critical  sliding  speed) 
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(Convection  affects  the  critical  wave  speed  but  has  no  significant  effect 
on  the  critical  sliding  speed) 
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Only  familifs  III  and  IV  tall  within  the  ri-ginn  of  practical  interest 
and  enclo'ie  an  area  wliere  instability  prevails.  Interestingly,  they 
also  enclose  a region  wliere  I',  as  determined  by  eq.  (2-b4),  is  complex. 
Inside  this  complex  region  the  actual  system  cannot  support  the  values 
ot  c supplied  in  the  computer  calculations  to  determine  11.  I’igure  (18) 
shows  the  I'gg-shaped  region  containing  the  values  of  c and  'a  for  wliich  U 
becomes  complex.  The  peripheral  points  are  the  solutions  representing, 
on  this  plane,  families  111  and  IV. 

figure  (19)  shows  all  the  families  of  solutions  in  terms  of  — and 

c 

(h<).  The  box-shaped  region  is  the  one  referred  to  in  the  preceding 
paragraph.  Using  the  algebraic  expressions  given  by  eqs.  (2-bb)  through 
(2-7h),  it  is  possible  Co  broadly  define  this  region.  fairiilies  I and  11 
coalesci“  into  the  horizontal  line  below  the  box. 


2.12  Instability  of  the  Non-moving  Perturbaicion 

In  the  conclusion  to  the  last  section,  it  was  noted  that  for  very 
small  wave  velocities  and  negligible  convection  and  clastic  effects,  a 
threshold  of  instability  does  exist  (Family  IV,  fig.  15  ) . To  determine 
on  which  side  of  this  line  stable  conditions  prevail,  one  may  check  for 
values  of  the  exponent  P close  to  and  on  either  side  of  this  threshold. 
In  the  energy  equation  (eq.  2-26)  , q^  and  q^^  remain  unchanged, 
while  q^^  and  q^^ , being  small  by  comparison,  drop  out.  Kquat  ion  (2-39) 
becomes 


h2=  0 

while  eq,  (2-38)  becomes 
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Figure  18.  Domain  of  c and  h in  Wiich  U is  Complex 
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Fit?ure  19.  Linos  of  Neutral  Stability  (W 
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Tlio  surface  wave-  h'  now  bocotnos 
h ' = h I -s  i n /X 

and,  since  elastic  effects  are  small  compared  to  the  thermal  effects, 
the  threshold  condition  is  Riven  by 

h-  = (2-781 


Using  eqs.  (2-50),  (2-77)  and  (2-78),  one  obtains 

. ^ i . JSj s' 

P 2 ° ^ ^ 

? r—  + K 

' 'Si 


(2-791 


Upon  transposing  the  second  term  in  eq,  (2-79)  to  the  right 
hand  side  and  then  squaring  both  sides,  one  can  obtain  a quadratic 
equation  for  non-zero  3 which  can  be  solved  to  yield 

.2,  7. .2.  “mSi 

..  °'mV  ''  bi  “ "h 

' 5,  ' 5,  ^ 


(2-80) 


The  + was  introduced  by  the  process  of  squaring  mentioned  above 
and  it  will  be  shown  subsequently  that  the  minus  sign  is  the  correct 
one . 

As  the  threshold  of  instability  is  of  principal  interest  here, 
one  may  seek  the  condition  that  gives  R=0.  In  eq . (2-79),  the  limit 
of  the  left-hand  side  for  p -+  0 is 
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Introclncinj;  this  in  cq.  (2-79),  one  obtains 


'f,'  ^ 

where  the  superscript  * refers  to  the  point  whore  P=()^  i.e,,  the  thres- 
hold of  instability.  It  shmild  be  noted  that  eq.  (2-81>  is  identical 
with  eq.  (2-74)  lor  U 

crit. 

A look  at  eq.  (2-80)  reveals  that  o will  be  non-zero  if  the  plus 
sign  is  chosen.  Moreover,  if  (U/h)  from  cq.  (2-81)  is  substituted  into 
eq.  (2-8()i,  becomes  zero  only  when  the  negative  sign  is  chosen.  There- 
fore in  eq . (2-80)  the  minus  sign  in  the  + is  the  correct  one. 

° may  be  expressed  in  terms  of  U/h.  Let 


(2-81) 


and 


/*-  fiT* 
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ThcMi  eq.  (2-80)  can  be*  shown  to  give 
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The  expression  under  the  radical  sign  in  eq.  (2-82) 
the  desired  degree  of  approximation  depending  on  the 


(2-82) 

may  be  expanded 
size  of  V/V  ; 


V , of  course,  is  constant  for  a given  geometry  and  given  materiel 
properties.  Knowing  the  sliding  speed  would  then  enable  one  to  deter- 
mine P readily  from  eq.  (2-82).  The  critical  or  threshold  condition  is 
represented  by 


y_ 

* 

V 
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which  also  corrospomls  to  Q=0.  ('or  V/V  « I,  cq.  (2-H2)  Kives 

2, 

5 = - (2-83) 

In  tig. (15),  I'amily  IV  represents  V/V  = 1,  On  either  side  of 
this  line,  on('  may  write 


V__ 
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+ e 


where  € 
h quat ion 
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V 

mav  bo  positive  (above  the 
(2-82)  can  then  be  written 

Si  r 2 f 

- 1 - 4(1  + €)  + ■ 1 


1 ine)  or  negat ive 


+ 8(1  + e) 
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(below  the  line). 
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If  c > 0,  3 > 0.  If  e.  < 0,  Q < 0.  Below  the  lino  representing  Family  IV 
in  fig. (151,  Che  system  is  stable.  Above  the  line,  instability  prevails. 


2.13  Remarks 

In  the  lubricated  face-type  seal,  a perturbation  represented  by  a 
steadily  moving  wave  may  become  unstable  within  a reasonable  operating 
range.  The  sliding  speed  corresponding  to  neutral  stability  (Family  IV, 
fig.  15  ) is  given  by 


U . 
cri  t. 


(2-84) 


Tile  wave  velocity  is  very  small  and  so  are  pressure  and  convection 

effects.  Instability  occurs  when  U > U . and  is  caused  by  viscous 

crit. 

heating  and  heat  conduction  into  the  metal.  The  disturbances  witli  the 
largest  wave-lengths  (smallest  K)  become  unstable  first. 
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Tln'  t-liri'slKild  condition  ^ivcn  by  eq.  (2-84)  exists  only  for  a 
range  of  t i Ini- th  i cknoss  h,  which  can  tu‘  dptermined  hy  using  eqs.  (2-70), 
(2-75),  (2-7  3)  and  (2-7b).  Tlie  lower  lioundary  of  this  range  is  given  hy 


_ Oh  ^iW  .,1/4 


l'.  Or 


M 


and  the  upper  boundary  by 


h = 


K 

oi  1 


Blok  L 15J  has  demonstrated  a different  mechanism  for  thermal 
instability  of  liquid  films  in  which  viscosity  varies  with  temperature. 
The  onset  of  such  an  instability  is  determined  by  a parameter  r^  which 
expresses  the  heat  flux  through  the  slower  (or  stationary)  wall  as  a 
fraction  of  the  total  heat  flux  across  the  film,  and  a dimensionless 
parameter  N defined  as 


N 


ik' 

dT, 


oil 


If  the  stationary  surface  is  a thermal  insulator,  r^=  0.  It  is  shown 
that,  for  r^=  0,  the  temperature  dependence  of  viscosity  docs  not  affect 
the  shear  stress  or  heat  flux  in  the  film  for  N <"  0.4.  I’ur  thermore , the 
critical  value  of  N corresponding  to  r^=  0 is  2. 

In  the  present  study,  has  been  assumed  constant  across  the  film, 
for  which  a corresponding  (-d^i/dT)  can  be  obtained  from  a chart  [lb]. 

For  the  threshold  condition  represented  by  Family  IV  (fig.  15  ),  it  is 
seen  that  N varies  between  0,016  and  0.26.  This  seems  to  justify  the 
constant-viscosity  formulation  of  the  present  study.  It  is  also  worthy 
of  note  that  the  present  values  of  N do  not  approach  the  critical. 


CHAPTKK  m 


KXPKRIMKNTAI,  OHSKRVATION  OF  THKRMOKLASTIC  INSTAlilLlTY 


i . 1 Hie  oh  joe  t i Vf 

K xptr imt'nts  wore  porlorniod  in  the  laboratory  t(j  verify  the  occur- 
rence of  thermoe  laHtic  instability  in  hydrodynamica 1 ly  lubricated  face 
seals.  The  prediction  of  Cliapter  II  was  based  on  an  analysis  of 
initially  flat  surfaces,  and  while  idealized  theories  often  provide 
enlightening  insiglit  to  complex  practical  problems,  the  proof  of  the 
validity  of  such  idealizations  lies  in  practical  demonstration.  No 
plans  were  made,  nor  any  hopes  entertained,  for  the  precise  quantitative 
corroboration  of  the  predicted  critical  speeds.  The  intention  was  to 
confirm,  through  repeatable  experiments,  that  the  phemonenon  predicted 
in  theory  was  existent  in  practice. 

Moving  hot  spots,  believed  to  be  caused  by  thermoelastic  instability, 
had  been  observed  and  filmed  by  Sibley  and  Allen  Ll.i  and  Kilaparti 
How  . IJ  had  briefly  reported  experimentally  measured  temperature  non- 
uni lormit  ies , clearly  originating  from  thermoelastic  instability,  in 
the  tw(5-pin  model.  In  a very  recent  paper,  Dow  and  Stockwell  have 

described  experiments  with  the  two-dimensional  scraper  model.  The  speci- 
mens used  included  commercial  apex  seals  used  in  Wanke I engines.  The 
edge  of  a scraper  was  pressed  aga inst  a rotating  drum  (fig.  20),  Inter- 
facial temperatures  were  mi'asured  by  means  of  embedded  thermocouples 
as  well  as  by  an  infrared  heat  sensor,  pointed  at  thi‘  interface  and 
sweeping  back  and  forth  along  it.  The  intcrfacial  clearance  was 
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Fluure  20.  Experimental  Set-Up  of  Uow  and  Stockwell  (Ref.  17) 
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ini’, isur  I'd  by  n I'oiipli'd  1 u.m’ r -phot  ocf  1 1 (.iimbi  niiL  i on  th.nL  also  swi’pL  along 
t-lio  interface  and  monitored  the  amount  of  liglit  passing  through.  The 
experiments  clearly  demonstrated  that,  above  a certain  speed,  tempera- 
ture and  clearance  nc^n-uni i ormit ies  of  large  magnitudes  occur  and  move 
along  the  interface.  Some  ol  these  observations  were  made  with  the 
interface  lubricated  by  .S/\b  10  oil  fed  by  a wick  ^presumably  boundary 
rather  than  hydrodynamic  lubrication). 

ihe  present  experiments  i.l8j  provide  evidence  of  thermoelastic 
instability  in  the  regime  of  hydrodynamic  lubrication. 

3.2  The  Experimental  Set-up 

The  rotating  face  consists  of  a cylindrical  metal  cup  with  a flat, 
thin  lip.  It  is  rigidly  held  on  a vertical  spindle  (see  fig.  21)  having 
freedom  of  motion  in  the  axial  direction.  The  stationary  face  consists 
of  a flat  pyrex  disc  held  in  gimbals.  Thus,  the  mean  film  thickness 
can  change  and  the  stationary  face  can  tilt  in  response  to  any  initial 
misalignment  or  to  changes  in  the  circumferential  pressure  profile. 

The  vertical  spindle  is  driven  through  a pulley  by  a d.c.  motor 
whose  speed  can  be  controlled  to  vary  smoothly  from  0 to  above  4000  r.p.m. 
The  spindle,  in  turn,  is  connected  by  a rubber  band  and  pulley  to  a 
small  generator  wtiose  voltage  output  is  calibrated  to  measure  tlie  spindle 
turning  speed. 

'Ihe  oil  is  held  in  a shallow  pool  on  tlie  glass.  Placed  in  two 
holes  in  tlie  glass,  and  rigidly  affixed  to  it,  are  two  eddy-current 
proximity  probes.  Tlie  first  of  these,  probe  A,  reads  the  distance 
between  its  tip  and  the  metal  face;  as  the  cup  turns  through  one  revo- 
lution, the  entire  circumferential  profile  of  the  film  thickness  is 


Figure  21.  Experimental  Set-Up  Showing  Probe  Uocations 
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. rtu'  si'concl  pn)l)i-  ( proln'  HI  is  loc'.itocl  on  tlu'  axis  of  Liu-  cup 
and  mcasiirt-s  vihrations  or  tiouncinp  ol  flic  cup  as  a wliolc.  fioth  probes 
are  connected  to  a dual-l>eani  oscilloscope  to  display  cyclic  patterns; 
each  is  also  connected  to  a hi  ^;li- impedance  d.c.  voltmeter  for  the 
recording  (it  mean  values. 

3 . 1 JJ le  Mea.siirc’mt  nt  of  Oil-lilm  Thickness 

It  has  belli  shown  in  (lliapter  11  that  the  critical  sliding  speed  is 
proportional  to  the  mean  f i Im-thickness , h (eq.  2-K4).  Previous  attempts 
at  nu'asuring  such  thin  films  have  had  mixed  successes  and  in  a few 
casi's  have  involved  complicated  ar r angemc'nts  . 

The  most  popular  means  of  measuring  f i Im-thickness  are  the  elec- 
trical resistance  method  _19,20j,  the  voltage  drop  or  discharge  volt- 
age method  .21j  and  the  capacitance  method  l22,23J.  These  techniques 
have  the  serious  drawback  that  cavitation  in  the  film  upsets  the  most 
careiul  calibration  of  the  instrumentation.  More  recently,  interfero- 
metric techniques  have  been  employed  wliich,  besides  calling  for 

very  sophisticated  set-ups,  can  do  no  better  than  measure  relative 
changes  in  the  f i Im- t h ickne ss . Winney  l24J  has  reported  the  develop- 
ment and  use  of  a rather  elaborate,  though  seemingly  effective,  probe 
that  de[)ends  on  the  variation  of  the  reluctance  of  a magnetic  circuit. 

Ihe  film  comprises  a small  gap  in  this  circuit  and  cavitation  induces 
little  error  because  oil  and  air  are  non- ferromagne t i c and  have  com- 
parable permeabilities. 

In  the  present  instance,  the  use  of  the  eddy-current  probes  pre- 
sents a major  difficulty.  Viscous  heating  in  the  oil  causes  thermal 
expansion  of  the  probes  themselves  (particularly  probe  A)  and  thus  not 
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I'ulv  alters  till'  pi'f-si’t  wap,  l)ui  also  shills  llio  ca  I i hrat  i on  curve  for 
tlu'  probe.  However,  the  calilnation  ecjnstant,  and  lienee  the  dynamic 
response  of  the  probe  to  tlu-  passing  metal  surface,  is  not  affected. 

A technique  has  bi-en  devised  to  circumvent  the  problem  of  zero-drift. 

An  initial  runnin>'-in  period,  at  low  speed,  is  allowed  until  the  d.c. 
voltmi-ter  reading  for  probe  A stabilizes.  Thereafter,  the  speed  is 
gradual Iv  increased  in  steps  and  lime  is  allowed  at  each  step  for 
stabilization  of  the  d.c.  voltmeter  reading  (indicating  mean  film- 
thickness)  for  probe  A.  Simultaneously,  at  a few  of  these  steps, 
after  the  d.c.  voltage  stabilize.s,  the  rotation  of  the  cup  i .s  suddenly 
arrested,  causing  it  to  drop  and  rest  on  the  glass.  Voltages  (d.c.) 
before  and  after  this  stoppage  are  noted  and  simultaneous  photographs 
of  the  oscilloscope  traces  are  taken.  To  this  voltage  drop  is  added 
the  mean  of  the  fluctuating  (a.c.)  component  as  indicated  by  the 
oscillogram  for  probe  A before  stoppage,  to  obtain  h.  The  calibration 
constant  for  probe  A is  0.0024  in. /volt  and  the  digital  voltmeter 
reads  millivolts,  thus  providing  a resolution  of  approximately  2,4  p.in. 

j.4  Procedure 

In  order  to  approximate  the  idealized  conditions  of  Chapter  1 1 as 
tar  as  practicable,  the  glass  and  the  metal  are  ground  and  polished  to 
a mirror  finish. 

Kach  probe,  connected  to  its  own  18-volt  power  supply  and  high- 
impedance  voltmeter,  is  calibrated  with  respect  to  the  metal  surface 
it  measures.  Figure  (22)  shows  such  a calibration  cvirve  for  a probe 
with  a tip  of  1/8  in.  diameter;  it  also  shows  a probable  Initial  setting 
when  mounted  on  the  test  apparatus.  The  two  probes,  A and  B,  are  found 


Figure  22.  Calibration 
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to  liavi'  sliKl'tly  ditleront  c a i i l>r  a t i on  constants  lor  a given  metal  cup. 

Tlif  cup  is  mounted  on  the  spindle  and  the  glass  slab  on  the  gimbals. 
Ihe  two  prolies  are  attached  and  set  to  chosen  initial  gaps.  The 
osc i 1 loscopi'  is  then  adjusted  so  that  the  signals  from  the  two  probes 
are  displayed  to  the  same  scale  (volts  per  inch  of  gap  deviation), 
i'he  oscilloscope  is  set  to  display  the  d.c.  signal  from  probe  H but  tlu- 
a.c.  from  probe  A.  Thus  one  observes  the  clianges  in  the  reading  of 
probe  B throughout  the  duration  of  tlie  run,  i.e.,  any  slow  rise  and  fall 
as  wi' 1 1 as  any  bounce  or  vibration  of  the  cup  as  a whole  are  displayed. 

■flu“  trace  for  probe  A,  by  contrast,  shows  only  the  surface  profile 
relative  to  the  mean. 

The  spindle  is  started  from  rest  and  the  speed  gradually 
increased,  as  explained  in  jJ.'i.  Voltage  readings  from  the  voltmeters 
are  recorded  and  simultaneous  photographs  of  the  oscilloscope  traces 
are  taken. 

No  attempt  is  made  to  control  the  film-thickness,  it  being  deter- 
mined by  the  operating  conditions.  Wlten  probe  B indicates  no  bouncing 
of  the  cup,  the  time -dependent  output  of  probe  A is  taken  as  the  measure- 
ment of  the  contour  of  the  surface  passing  over  it.  Temperature-dependent 
variations  in  such  properties  as  the  viscosity  of  thi-  oil  and  coef- 
ficients of  thermal  conductivity  and  expansion  of  the  metal  are  not 
considered.  Therefore,  no  attempt  is  made  to  measure  temperature 
variations.  Instead,  the  bulk  temperature  of  the  oil  is  measured  after 
shut-down  by  a simple  thermometer  so  as  to  obtain  an  estimate  of  the 
overall  viscosity  of  the  oil  (SAT  10)  from  viscosity  charts.  Pressures 
are  not  measured.  The  only  evidence  of  thermoelastic  instability  being 
sought  is  macroscopic  deformation  of  the  surface. 
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).5  Results  an  cl  Observation s 

Scidcifn  .incl  very  dramatic  de  t ormat  i ons  of  tlie  surface  are  fecund  to 
occur  at  a certain  speed,  interpreted  as  tlie  critical.  I’rior  to  the 
onset  of  the  instability,  the  surtace  wave  grows  steadily  in  magnitude, 
w!i  i le  the  oscillogram  lor  prohe  11  shows  no  hounce  or  vibration.  1 he 
film-thickness  prior  to  the  occurrence  of  the  instability  is  measured 
a',  outlined  in  jl.3  and  is  plotted  against  the  critical  sliding  speed 
as  shown  in  fig.  (2Jj.  In  this  figure',  the  data  points  are  experimentally 
obtained,  wliile  the  closed  regime  of  instability  is  drawn  on  tlie  basis 
of  the  deri'.ations  of  Chapter  11,  hut  for  the  geomcitry  used  in  the 
experiment.  i'he  physical  parameters  and  material  properties  used  arc- 
listed  in  Appendix  (1. 

t in  tfie  course'  of  the  expe  r imc'nt  s , the  well  known  phenomenon  of 

inward  pumping  (sec‘,  for  instance',  ret.  L25J)  has  been  repeatedly 
observed.  1 wo  some-what  unexpe-cted  observat  ic.)ns  are  that,  for  a given 
load,  a wider  seal  face  requires  a higher  starting  torque  and  an 
abundance  of  oil  in  the  poc'l  usually  leads  to  higher  values  of  film- 
thickness  than  when  there  is  less  oil.  These  effects  perhaps  merit 
further  studies  and  some  explanation,  but  are  being  ignored  here  because 
the  critical  condition  alone  is  of  interest. 

In  all  the  oscillograms  that  follow,  the  upper  trace  represents 
probe  B,  the  lower  trace  probe  A.  An  upward  trend  represents,  in  tlie 
upper  trace,  a rising  of  the  cup  from  the  glass;  in  the  lower  trace  it 
represents  a prot ruberance  on  the  face  of  the  cup. 

iigure  (24)  shows  the  evolution  of  a surface  profile  from  a low, 
subcritical  speed  (a)  to  a supercritical  speed  (e).  Large  spikes  or 


Figure  24.  Oscillograms  Showing  Surface  Profile  at  Different  Sliding  Speeds.  i = One  Revolution 


r i I i (■ --  ;iro  sei'ii  at  tlu’  supc  r<  i' i t it;  .i  I .speed.  As  an  asperity  passes 
liver  prolie  A and  dips  into  the  sliv.ht  depression  above  the  probe,  it 
sets  oil  a vibration  that  is  clearly  indicated  by  tlie  upper  trace  in 
tig.  t24e).  This  is  damped  out  until  the  next  asperity  passes  over 
probe  A and  triggers  another  vibration.  The  lower  trace  records  a 
vibration  of  the  same  phase  and  magnitude.  i'he  photographs  of  fig.  (25) 
duplicate  these  aspects.  That  tin.'  asperities  of  fig.  (24i  ) are  of 
thermal  origin,  and  not  the  man i 1 es ta t i on  of  permanent  surface  damage, 
is  proved  by  fig.  (24f).  This  picture,  taken  after  slowing  down  fol- 
lowing lirief  operation  in  the  deformed  state,  shows  that  the  asperities 
have  all  but  disappeared.  More-over,  the  surface  profile  closely 
resembles  that  in  fig.  (24n),  taken  at  the  same  sliding  speed  but  during 
--tart  -up . 

Kxaminatioti  of  the  metal  face  at  the  end  of  the  experiment  indicate 
tliat  the  asperities  break  through  the  oil  film  and  rub  directly  against 
the  glass.  If  only  one  asperity  grows  and  touches  tin-  glass,  static 
equililirium  requires  that  the  glass  plate  must  tilt  and  touch  the  metal 
lace  at  a di  ;une  t r ica  1 ly  opposite  point.  The  metal  cup  tlien  rests  on 
tfie  glass  on  two  "legs”,  the  minimum  required  for  dynamic  equilibrium. 
Sometimes  three  "legs"  are  formed.  l-igure  (2b)  shows  an  aluminum  face 
with  three  black  spots  on  it,  the  locations  of  three  burnt  asperities. 

A very  interesting  feature  of  the  operation  at  speeds  just  below 
the  critical  is  the  growth  and  movement  of  small  surface  disturbances. 
Witli  the  cup  running  at  constant  speed,  small  growths  are  seen  to  arise 
and  traverse  the  face.  The  effect  is  very  graphic  when  the  oscilloscope 
trigger  is  set  to  commence  tlie  sweep  at  a particular  point  along  the 


rams  Showing  Surface  Profile  at  Different  Sliding  Speeds.  Arrow  Shows 
Blemish,  Present  at  the  Outset,  Relative  to  Wliich  Small  Surface  Irregu 
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c i rcumteronco  o(  ihe  cup.  It  lias  lu'oti  repeatedly  obscTved  and  is  clearly 
a surface  phencmienon . It  is  illustrated  by  fig.  (24b, c,d).  in  fig.  (25), 
the  arrow  marks  a surface  blemish  that  was  present  at  the  very  outset 
and  provides  a convenient  reference  point  relative  to  which  moving 
features  can  be  observed.  The  fact  that  tlie  upper  trace  indicates  no 
bounce  prov<'s  that  this  is,  indeed,  a surface  effect. 

The  close  conformity  betwi'en  predicted  and  measured  results  (fig.  23) 
is  startling.  There  is  some  uncertainty  as  to  whether  the  circumference 
of  the  cup  contains  one  surface  wave  or  two  (see,  for  instance,  figs. 

24c  and  ei.  In  tig.  (23),  the  solid  line  is  for  a single  wave  (<  = .66  In 
and  the  dashed  line*  for  two  waves  (K  = 1.33  in  ^)  around  the  circumference 
The  pertinent  dimensions  and  properties  are  listed  in  Apptuidix  C. 

On  the  basis  of  fig.  (23),  it  is  believed  that  the  analysis  pre- 
sented in  Chapter  II,  although  idealized,  does  take  proper  account  of 
the  principal  effects  at  work. 


CHAITKR  IV 


INSTAHll.lTY  /\NAI.YSIS  FOR  INITIAI.FY  WAVY  SURFACKS 


■4 . 1 riif  Ho  elf  1 

It  was  ln'lioved,  even  a few  years  ago,  tfiat  physical  contact  betweea 
the  laces  of  a seal  was  responsible  for  tile  prevention  of  leakage. 

Almost  invariably,  one  face  of  a seal  used  to  be  made  til  graphite,  wliose 
se 1 t - lubr icat ing  property  was  possibly  responsible  for  the  satisfactory 
performance  of  face  seals.  Witli  the  advent  of  higher  and  higher  duty 
requirements,  liowever,  the  soft  graphite  faces  began  to  perform  witli 
greater  and  griater  unpredictability  and  harder  alternative  materials 
began  to  be  sought.  Research  on  seals  was  stepped  up  and  gradually  the 
realization  dawned  that  some  form  of  face  lubrication  was  required  for 
satisfactory  operation.  There  was  strong  evidence  linking  failure  of 
lubrication  to  seal  failure. 

Numerous  theories  have  since  been  advanced  in  efforts  to  explain 
how  a face  seal  supports  load  and  effects  sealing.  dne  such  theory, 
inadequatt' 1 y understood,  is  boundary  lubrication,  which  is  characterized 
by  a degree  of  direct  contact  between  the  sliding  surfaces  and  can  be 
envisaged  under  conditions  of  very  high  load  or  when  the  sealed  fluid  is 
highly  volatile.  Nau  t26j  has  reviewed  the  theories  of  lubricated  face 
seals  and  concludes  that,  in  general,  face  seals  operate  in  the  hydro- 
dynamic  reglm»',  characterized  by  a coherent  lubricating  film  separating 
the  two  surfaces. 
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Hamilton  et  al  [27j  have  presented  a theory  of  face  lubrication 
depending  on  micro-asperities  on  the  surface;  the  fluid  flow  is  stated 
to  generate  positive  pressures  on  the  upstream  side  of  these  micro- 
asperities and  cavitation  pockets  on  the  downstream  side,  thereby 
generating  a net  positive  pressure.  Reiner  [28]  has  tried  to  explain 
the  existence  of  a hydrodynamic  film  through  the  "inward-pumping"  effect. 
This  phenomenon,  frequently  observed,  consists  of  a radial  inflow  across 
the  faces  of  a seal;  it  implies  some  leakage  and  Reiner  has  attempted 
to  explain  it  in  terms  of  non-Newtonian  effects  in  the  film.  Non- 
Newtonian  models  are  difficult  to  formulate,  but  Tanner  [29]  has  shown 
that  they  cannot  adequately  explain  the  ways  in  which  a face-seal  per- 
forms. Another  theory  postulates  the  generation  of  hydrodynamic 
pressure  by  faces  that  either  have  an  angular  misalignment  or  vibrate 
or  oscillate  with  respect  to  each  other  [30, 3l]. 

Other  theijrii's  consist  of  the  "thermal  wedge"  in  which  circum- 
ferential temperature  variations  cause  density  and  thence  pressure 
variations,  "pressure-viscosity"  effects,  radial  eccentricity,  and 
simply  a radial  hydrostatic  pressure  gradient  between  the  sealed  fluid 
and  its  "surroundings".  More  recently,  attention  has  been  focused  on 
the  strong  likelihood  of  hydrodynamic  pressure  being  generated  by  a 
mechanism  that  fits  the  classical  theory  of  hydrodynamic  lubrication 
presented  by  Reynolds  in  1886.  Stanghan-llatch  and  others  [25,32]  have 
•i  'wn  that  hydrodynamic  pressures  can  be  generated  by  wavy  surfaces 
' th'  ifHTation  of  seals  with  wavy  surfaces  can  be  explained  by 

< :■  "r  V , 

• ’ . . It  lias  been  assumed  that  the  two  faces  are  perfectly 


.1 


I'i.it  ••naliles  the  use  of  analytical  deductions. 


it  is  a highly  idealized  configuration  that  does  not  occur  in  practice. 
Because  t)ie  nominal  film  thickness  h is  very  small  in  real  face  seals, 
even  a very  slight  waviness  is  likely  to  be  significant  when  compared 
with  h.  Such  waviness  may  comprise,  in  some  instances,  an  angular 
misalignment  of  the  two  faces.  As  mentioned  above,  there  are  several 
possible  modes  of  operation  for  a face  seal,  ranging  from  boundary 
lubrication  to  hydrodynamic  action.  A large  number  of  applications  fall 
in  the  latter  category,  and  the  generation  of  hydrodynamic  pressures 
depends  on  surface  waviness.  While  the  flat  face  idealization  may 
hold  in  a few  real  cases  involving  very  small  waviness  and  light  face 
loading  or  a hydrostatic  pressure  gradient,  an  analysis  of  the  wavy 
surface  is  essential  to  the  understanding  of  thermoelastic  effects  in 
most  real  seals. 

The  upper  surface  in  fig.  (27)  is  a good  thermal  conductor  with  a 
waviness  defined  by  h^  at  a given  operating  condition.  This  waviness 
is  not  the  initial,  machined,  surface  profile,  but  rather  the  modified 
profile  comprising  the  initial  shape  as  well  as  subsequent  thermal  and 
elastic  distortions  compatible  with  the  operating  conditions.  The 
stationary  bottom  surface  is  a flat  thermal  insulator.  This  surface 
is  chosen  as  flat  because  it  does  not  undergo  tliermal  distortions  and 
any  initial  waviness  on  it  will  remain  the  same  because  elastic  effects 
are  usually  negligibly  small.  The  lace-load  f,  the  sliding  speed  U , 
h^  and  h together  represent  a quasi-static  state  or  operating  condition. 
Since  f and  U are  nominally  fixed  in  operation,  departures  from  the 
quasi-static  state  will  involve  changes  in  either  h^  or  h,  or  in  both. 
Thus,  in  order  to  determine  whether  thermoelastic  instabilities  can 
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occur,  one  may  impose, as  before,  a small  surface  perturbation  h'  on 
and  examine  whether  the  resulting  changes  in  heat  and  pressure  generation 
tend  to  increase  the  amplitude  of  h'.  If  so,  the  operating  condition 
will  be  said  to  define  an  unstable  state.  If  not,  the  seal  will  be 
said  to  be  operating  in  a thermoe lastical ly  stable  regime. 

Most  real  face-type  seals  have  low  aspect  ratios,  i.e.,  the  face- 
width  is  small  compared  to  the  inner  radius.  Therefore,  the  narrow- 
face  geometry  will  be  considered  here.  It  is  then  possible  to  replace 
the  cylinders  (fig.  4)  with  flat  blades  as  in  Chapter  II.  Furthermore, 
because  the  pressure  gradients  in  the  x-direction  are  considerably 
smaller  than  those  in  the  z-direction,  a fair  approximation  of  the  flow 
is  obtained  by  considering  it  to  be  purely  Couette.  This  may  not  seem 
very  obvious  for  sizable  values  of  jh^j/h,  but  is  true  at  least  in 
those  cases  where  the  surface  waviness  has  a large  wavelength  and  the 
face  width  is  small  by  comparison  [33]. 

The  velocity  distribution  for  the  x-direction  is  then  given 
simply  by 


where 

h(x)  = h + fi  + h (x)  + h'(x)  (4-2) 

o 

(i'(x)  being  the  small  perturbation  whose  decay  or  growth  will  indicate 
the  stability  or  otherwise  of  the  system.  The  introduction  of  h'  will 
alter  the  hydrodynamic  pressure  distribution  in  tiie  film,  tlms  clianging 
the  net  pressure  or  "load  support"  generated  in  the  film.  Since  the 
face- loading  f is  constant  for  a given  system  and  must  be  balanced  by 
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the  net  pressure  in  the  film,  the  two  faces  will  tend  to  part  or  come 
closer  to  negate  the  increase  or  decrease  of  the  net  pressure.  This 
relative  axial  motion  between  the  two  faces  is  represented,  in  eq.  (4-2), 
by  h,  which  signifies  a small  change  in  h,  occasioned  by  the  small 
cliange  h'  in  h^.  It  is  assumed  here  that  both  faces  are  held  in  rigid 
mountings  that  permit  axial  motion  without  tilt.  If  one  face  is  flexibly 
mounted,  as  in  the  experiments  described  in  Chapter  111,  it  will  tilt 
in  response  to  any  asymmetry  in  the  pressure  wave  resulting  from  the 
introduction  of  h'.  The  tilting  case  will  be  dealt  with  later. 


(4-4) 
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If  ont'  defines  q = q + q'  where  q denotes  the  zero  order  heating 
g o o 

obtained  in  the  unperturbed  operating  condition  and  q'  the  perturbation 
resulting  from  the  introduction  of  h',  one  can  write,  from  eq,  (4-4), 


q = — 
^o  — 


(4-5) 


h + ti 


::  (h'  + fi) 


(4-6) 


(h  + h )■ 
o 


4.3  Heat  Conduction  and  Convection 

The  heat  generated  in  the  film  through  the  process  of  viscous  dis- 
sipation is  removed  through  convection  and  conduction.  It  was  shown  in 
Chapter  11  (see  also  Appendix  A)  that,  for  a narrow-face  seal,  con- 
vection has  little  effect  on  the  temperature  distribution  in  the  metal. 
Georgopoulos  l34j  has  proved  this  through  a rigorous  solution  of  the 
energy  equation  incorporating  the  convection  term.  It  has  been  found 
that,  for  the  narrow-face  seal,  tlie  convected  heat  is  only  a very  small 
percentage  of  tliat  conducted.  In  view  of  the  above,  it  will  be  assumed 
here  that  all  the  heat  generated  is  conducted  into  the  metal  face.  Thus, 
eqs.  (4-4),  (4-5)  and  (4-6)  all  represent  the  heat  input  into  the  metal. 


4.4  Thermal  Deformation 

The  relationship  between  the  heat  input  to  a surface  and  the  corre- 
sponding surface  curvature  is  given  by 


d^v 

dx^ 


•Si 


(4-7) 


where  v denotes  the  thermal  deformation  of  the  surface.  This  applies  to 


plane-stress  problems  and  was  shown  by  Burton  et  al  ; an  erroneous 
multiplier  of  2 was  inadvertently  introduced  and  later  corrected  [35j. 
In  order  that  eq.  (4-7)  may  reflect  the  fact  that  uniform  heating  will 
lead  to  zero  curvature,  an  extra  term  should  be  introduced.  In  terms 
of  tlie  present  nomenclature,  eq.  (4-7)  may  be  rewritten  as 


^ ('th>  v 


dx 


•Si 


- c 


(4-8) 


whe  re 


o ‘Si 


dx 


(4-9) 


Evidently,  eqs.  (4-8)  and  (4-9)  are  written  in  terms  of  perturbation 
quantities.  The  integration  from  x = 0 to  x = 2 implies  summing  around 
the  circumference  of  the  real  seal.  Using  eqs.  (4-6)  and  (4-9),  in 
eq.  (4-8),  one  obtains 

S ^\h‘  = • T — 2 

dx^  ^ Si  (h  + h )^ 

o 


. M .,2  1 h'  + h 
+ ^ P.U  j — : 

^ o (h  + h )■ 

o 


dx 


(4-10) 


which,  if  integrated  twice,  will  give 


4.5  Hydrodynamic  Pressure 

The  tangential  pressure  distribution  in  the  film  is  obtained  as 
outlined  in  (Tiapter  II.  As  mentioned  earlier  Lsee  ref.  33J,  greater 
care  is  called  for,  in  the  case  of  the  wavy  surface,  in  applying  the 
assumption  of  ^p/Sx  « ^p/dz..  It  should  be  emphasized,  once  again,  that 
such  an  assumption  is  quite  valid  whenever  the  face-width  is  very  small 


ccmipared  to  the  wavelength  of  the  surface  wave.  Similar  considerations 
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the  basis  i>f  the  well-known  "short  bearing"  theory.  (Cameron  ilbj 
s shown  that  for  \/L  8,  the  sliort  bearing  equations  compare  very 

well  with  results  obtained  from  a numerical  solution  of  the  full  Reynold's 
equal  ion.  The  total  pressure  is  given  by 


P 


3ui;  dh  ! 2 

3 dx 
h 


(4-11) 


Also 


*^max  ^ 


and  the  average  pressure 


IS  given 


P = 3 


by  (fig. 

P 

max 


14) 


so  that  one  obtains 


- = . liLL  i2  ^ 
2h'  ' 


(4-12) 


in  wliich  h is  as  defined  by  eq.  (4-2). 

It  should  now  be  pointed  out  that,  for  a symmetric  surface  shape 


such  as  a sinusoidal  wave,  the  pressure  wave  given  by  eq.  (4-12)  will 
also  be  symmetric  because  both  h(x)  and  dh/dx  are  symmetric.  Con- 
sequently, one  will  obtain  positive  pressures  in  the  convergent  portion^ 
of  the  liquid  film,  and  symmetrically  opposite  (negative)  pressures 
in  the  divergent  portions,  as  indicated  in  fig.  (28).  The  negative 
and  positive  parts  of  the  pressure  wave  then  cancel  each  other  and 
the  net  pressure  is  zero,  implying  tliat  no  load  can  be  supported. 

In  trying  to  explain  how  a face  seal  supports  load,  Stanghan-Batch 
et  al  132]  reported  experiments  which  strongly  suggested  the  existence 
of  a hydrodynamic  film  between  the  two  faces.  The  measured  pressure 


waves  showed  positive  peaks  where  the  film  converged,  and  essentially 


Figure  28.  Circumferential  Pressure  Distribution 

a)  Sinusoidal  Surface  Shape  giving  symmetric  h(x) 

b)  Corresponding  theoretical  pressure  wave 

c)  Pressure  wave  with  negative  regions  suppressed. 
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zero  pressure  where  it  diverged.  Staiighan-Batch  and  J ny  subsequently 
published  a hydrodynamic  theory  for  face  seals  L25].  It  was  shown  that 
if  one  of  tlie  two  faces  was  wavy,  a hydrodynamic  pressure  wave  was  gen- 
erated, but  tliat  cavitation  occurred  in  the  diverging  portions  of  the 
film  and  extended  into  the  converging  regions  as  wi?  1 1 (fig.  29).  While 
it  is  known  that  filaments  of  fluid  continue  to  flow  through  and  around 
the  cavity  pockets,  such  flow  is  difficult  to  model  or  describe  mathe- 
matically. Most  commercial  oils  are  known  to  release  dissolved  air  at 
about  four- fifths  of  atmospheric  pressure.  The  authors  found  that,  by 
analytically  calculating  the  pressures  and  then  simply  suppressing 
the  negative  parts,  a strikingly  good  approximation  to  experimentally 
measured  values  could  be  obtained,  as  shown  in  fig.  (30),  This  infor- 
mation will  be  used  in  the  present  study. 

Kquation  (4-12)  can  be  expanded  as 
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(4-13) 


If  one  defines  p = p + p'  where  p denotes  the  zero  order  (operating) 
pressure  and  p'  the  first  order  (perturbation)  pressure,  one  obtains, 
f rom  eq . (4-13), 


Roto'. -on  of 
KQvy  *ocf 


Figure  29.  Typical  Location  of  Cavitation  Pockets 
(Stanghan-Batch  and  Iny,  ref.  25) 
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(4-15) 


where  all  second  order  terms  have  been  discarded  as  being  small. 

Since  denotes  the  film  pressures  in  the  unperturbed  operating 
condition,  its  negative  values  indicate  the  zones  of  cavitation.  It  will 
be  assumed,  as  stated  earlier,  that  in  the  cavitating  zones  the  film, 
for  analytical  purposes,  does  not  exist.  Everywhere  else,  the  film 
exists  and  has  positive  pressures.  The  perturbation  p'  represents  a 
small  change  in  p^  and  has  meaning  only  in  those  regions  where  the  film 
exists.  In  the  ca  itating  regions,  p'  must  necessarily  be  zero,  while 
in  the  positive  pressure  regions  even  negative  values  of  p'  are  per- 
missible because  they  will  be  small  and  have  little  effect  on  the  larger 

positive  values  of  p . 

o 

One  may  consider  for  a moment  the  case  of  a seal  whose  faces  are 
so  mounted  that  there  can  be  no  relative  axial  motion  between  them,  i.e., 
h 0.  It  will  be  seen  that  the  perturbation  pressure  p"  for  such  a 
set-up  will  be 
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2(h  + h 

o 


•^dh’ 

\dx 


3h' 
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(4-16) 


Using  eq . (4-16i  in  eq.  (4-15),  one  may  write 


p'  = p"  - p 

(4-17) 

whe  re 

2 

- 3^ilIL 

dh 

u 2 

(4-18) 

P - 4 

2(h  + h ) 
o 

dx 

It  may  then  be  argued  that  since  p"  has  a wave  form,  it  can  be  thought 
of  as  composed  of  a zero-average  wave  and  a constant.  By  definition, 
then,  the  zero-average  component  is  incapable  of  supporting  load,  but 
capable  of  producing  face  distortion.  The  constant  (uniform)  component, 
on  tlic  other  hand,  is  incapable  of  causing  distortion,  but  produces  a 
net  force  (positive  or  negative).  Since,  by  definition,  p^  supports 
the  fixed  load  f on  the  face,  any  additional  axial  force  generated  withi 
the  film  must  be  neutralized  by  relative  axial  motion  between  the  faces 
(wfiere  the  mechanical  set-up  allows  such  motion,  as  in  almost  all  cases) 
and  such  motion  is  represented  by  h.  Thus,  if  one  subtracts  from  p" 
its  constant  component  (as  represented  by  the  right-hand  side  of  eq . 
(4-17)),  one  is  left  with  its  zero  average  component  which,  according 
to  eq.  (4-17),  is  p' , That  means 


I 


1 ^ 
I 0 


xs=0 


P 


' dx  = 0 


so  that  eq.  (4-17)  gives 
i 


I 

I - 


x=0 


p dx  = - h 


'’i  ’’  3p.UL 


dh 


^ ' T /Ik  V. 

o 2 (h  + h ) 
o 


o 1.  ^ 

° = I . 

o 


(4-19) 


In  eq.  (4-19),  the  right  hand  side  is  known  because  p"  is  defined,  in 
eq.  (4-lb),  in  terms  of  operating  conditions  and  an  impressed  h'.  The 
only  unknown  on  the  left  hand  side  is  h,  which  can  then  be  calculated. 

In  view  of  the  fact  that  the  film  is  supposedly  non-existent  in  portions 
of  the  Interval  of  integration,  the  definite  Integrals  require  some 
qualification,  but  that  will  be  discussed  later. 

For  the  present,  it  is  sufficient  to  recognize  that  the  pressures 


that  cause  face  distortion  are  given  by  p'  and  that  these  can  be 
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dctormiiicd  hy  first  obtaininR  p"  I roin  ('q.  (4-16),  then  li  from  eq . (4-19), 
and  then  using  eq.  (4-171  for  p'. 

4 . ti  Klastic  Deformation 

As  pointed  out  in  the  priueding  section,  elastic  distortion  of  the 
face  is  caused  hy  p'.  The  re  lat  ionsli  ip  between  pressure  and  the 
deflection  produced  by  it  is  Lib] 

?=b 

9 f 

v(x)  = - — p ( s)  Jtn  I X- ~ 1 d? 

r,  • - 
■ = a 


where  the  interval  a to  b defines  the  extent  of  the  loading  zone,  and  | 
the  location  of  a load  element,  as  in  fig.  (31).  This  relationship  is 
valid  everywhere  except  for  ? = x,  and  so  the  region  (x-e)  < ' (x  + e) 
\ is  excluded  from  it.  Differentiating  under  the  integral  sign  by 

[.eibnitz's  rule,  one  finds 


As  the  region  around 


^ ^ L.  P'Aai..  d' 

dx^  “a  lx-?l^ 

? = X is  of  interest  here. 


this  region  is  treated 


independently  by  exploiting  the  linearity  of  the  load-displacement  re la- 

I 


tionship  and  adding  the  component  of  curvature  produced  by  a uniform  load 
over  this  region.  This  is  done  by  differentiating  the  corresponding  dis- 
placement curve  '„36j  and  found  to  be 


cp(e)  = - 


4jV 

nEe 


It  follows  that,  for  small  e. 


d^v 


d^i). 


2 2 
dx  dx 


TIT,  ijj  _ I , I 2 


h-?r 

, 2_  ) 


X + e X' 


21111 


d? 


. 


(4-20) 
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. 7 Nortnaliy.al  ion  of  Variables 

To  reduco  all  the  preceding  equations  to  non-dimensional  form,  the 
followinK  dimensionless  variables  will  he  used: 


X = I'.x 

Y(X)  = ^ 
h 

H 

li 


H(X)  = 


h ( X ) + h 
o 


and 


X = 


crit 


wliere 


U •-  = hK 
crit  , pD" 


M 


as  in  eq . (2-84 ) , 


Also, 


th 


th 

h 


K 


further,  if  both  the  surface  wave  h^  and  the  perturbation  h'  are  defined 
to  have  a wavelength  equal  to  the  circumference  i of  the  seal, 

I = ^ 


so  that  the  dimensional  interval  i is  replaced  by  the  non-dimensional 
interval 

2TT  = Kl 


as  shown  In  fig.  (32). 
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4.7  Normalization  of  Variables 

To  rcduco  all  the  preceding  equations  to  non-dimensional  form,  the 
following  dimensionless  variables  will  be  used: 


X = <x 


Y(X)  = - 


11  (X)  = 


h (x)  + h 
o 


X = 


and 


wfiere 


cr  1 1 


cr  1 1 po* 


M 


as  in  eq.  (2-84) , 


Also, 


th 


th 

h 


E 


further,  if  both  the  surface  wave  h^  and  the  perturbation  h'  are  defined 
to  have  a wavelength  equal  to  the  circumference  i of  the  seal, 

so  that  the  dimensional  interval  A is  replaced  by  the  non-dimensional 
interval 

2tt  = Ki 


as  shown  in  fig.  (32). 


90 


t follows,  from  the  above,  that 


^ dX 
dx  = — 


and 


9 9 

d^h'  -2  d'-y 
= h<  — 

dx"  dX 


etc . 


Kquation  (4-10)  then  reduces  to 

.2 


dx' 


d^A 


th 


a 

V 

2 , 


IV  0,11 

H \ 1 


^ Zl) 


'Y 


il 


- 2'n  . , 2 ' l) 

K^h  o H W ' 


dX 


2n 


dX 


u ' , r,  , u Y + 11 


dX 


(4- 


In  non-dimensionalixing  the  equations  for  pressure,  it  should  be 
noted  that 

# = 0 
dx 


so  that 


dh  d(h  + h) 
o o 


dx 


dx 


Then,  eq.  (4-14)  may  be  written  as 


■■o  „3  dX 


or. 


^ U_  J|_  ^ 

Po  " ' .,r2  U _ -vpa^  dX 
2h  crit  M H 


iZ*}  I ^*Sl  u dll 
2h  -'"m  h'"" 


and,  since  — = W as  defined  in  Chapter  II, 


21) 
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h [ M _ ^ 

’o  W = ■ ^3  dX 


Dt' fining  tlu'  non-dimensional  zero-order  pressure  as 


r , tT* 
h ^ M 

^o  W . ^tK^ 


(U-12) 


one  finally  go  t s 


f)u  dll 
■ - „3  dX 


(4-23) 


Similarly,  by  non-dimens iona 1 izing  eq.  (4-16),  one  obtains 


^ /dy  Y dll 
^j3  -dX  ■ ^ II  dX; 


(4-24) 


where 


p..  = p..  ii  '_il 
P WV.K^ 


(4-25) 


From  eq.  (4-18),  it  follows  that 


p _ . JJil  Fl 

- ,,4  dX 


(4-26) 


whe  re 


— , o 

D - h IH 

'■  = ” -w 


(4-27) 


and,  from  eq.  (4-17), 


p'  = p"  - p 


(4-28) 


where 


. h lIh 

" " ''“Si 


(4-29) 


For  the  non-dimensional  elastic  deformation,  one  has,  from  eq.  (4-20), 


. 2- 1!  ' ffa  K l*^'  ) 

■ dX^  |X-xl^- 


where  1 has  the  meaning  defined  on  pg.  86. 
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or  , 


dX 


•■>  o > ' ')  T * 2 y 


' » <V 

M 


x-xl 


It  will  he  recalled  from  Chapter  II  that  the  lower  limit  of  f i Im- th ickne ss 
for  which  instability  was  predicted  was  given  by 


t 

'■l  = c 


“2 

K t>  iC 

M 


1/4 


if  one  defines 


— , one  finds 
h 


T 

dX- 


(X) 

2 ^6  n 1x-xl^ 


(4-30) 


In  cliapter  11,  it  has  been  shown  that  hj^  is  representative  of  the  edge 
of  the  domain  in  which  elastic  effects  are  consequential  to  instability. 

In  the  present  case,  therefore,  one  would  expect  elastic  deformation 
to  be  important  for  Ii|^  = 1. 

4,8  Stability  Analysis 

As  in  Chapter  II,  the  question  of  whether  or  not  the  system  is  in 
thermoi- last ic  stability  will  again  be  posed  in  terms  of  a perturbation. 

Will  the  introduction  of  the  small  perturbation  Y(X)  result  in  thermal 
and  elastic  (4|,)  deformations  whose  combined  effect  is  greater 
in  amplitude  than  1y|?  If  so,  the  chosen  operating  condition  is  un- 
stable; if  not,  it  is  stable.  As  before,  three  conditions  can  be 


de  f ined : 


9J 


stable  : 

>\h 

+ 4^.1 

• lYl 

neutrally  stabliT 

'■^th 

1 = 1y! 

14-31) 

unstable : 

l^h 

> hi 

'llu>  second  condition  represents  the  so-called  threshold  of  instability 
that  was  sought  in  Ciiapter  II  bv  solving  the  corresponding  equation  for 
the  case  of  initially  flat  surfaces.  In  the  present  case,  such  a solu- 
tion by  analytical  means  is  unfeasible  in  view  of  the  non-linearities 
present.  Therefore,  Y(X)  will  be  supplied  and  and  A^^fX)  will  be 

separately  computed  from  eqs.  (4-21)  and  (4-30)  respectively.  Then 
Uti^(X)  and  li^(X)  will  be  summed  and  1 I will  be  compared  with 

jv]  to  see  which  one  of  conditions  (4-31)  prevails. 

In  accounting  for  cavitation,  it  will  be  noted  that  some  flow  does 
occur  in  the  cavitation  zones,  but  that  the  heating  effects  of  that  flow 
are  uncertain.  One  may  say,  however,  that  such  effects  will  lie  between 
the  extremes  of  full-film  heat  generation  and  zero  heating  (corresponding 
to  a truly  non-existent  film).  Stability  analyses  will  be  made  for 
both  extremes  and  the  true  situation  will  be  expected  to  fall  in  between. 

The  procedure  is  outlined  in  fig.  (33).  For  the  purpose  of  numerical 
computation,  the  derivatives  are  written  in  terms  of  the  central  dif- 
ference formula,  while  the  integrations  are  done  by  the  trapezoidal 
rule.  Since  the  film  equations  are  meaningless  wliere  the  film  does  not 
exist,  the  integrals  betwi-en  0 and  2rr  are  treated  as  improper  integrals 
and  the  integrations  are  carried  out  only  over  those  regions  of  the 
domain  in  which  the  film  exists  (analogous  to  the  "half  Sommerfeld"  con- 
di t ion) . 


J 


Kigurt-  'J3.  Block  Diagram  of  Instability  Analysis 
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I- qua  I ion  (4-21)  is  expaiuiiul  as 


it-l 


( I)X)^ 


Y.  2 n Y. 

2 _i  u 1 

'■  2 (n-l)IJX  2 

il.  1=0  II. 

1 1 


(4-32) 


whoro  n c.l('not('s  tlio  nunihor  ol  pivoLr.  I points,  l)X  tho  pivotal  spacing  and 
V ri'pri.  sent  s the  correcti'd  (or  new) 


(Y) 

new 


(Y) 


old 


+ II 


If  one  wiites  eq.  (4-12)  for  all  the  i's,  one  gi'ts  n equations,  of 
which  two  have  to  be  replaced  by  boundary  conditions.  The  first  boundary 
condition  derives  from  the  fact  that  must  be  single  valued  around  the 

seal  circumference.  Thus 


(\h\  = (\h^ 

1 n 


(4-33) 


The  other  boundary  condition  is 


(a  ) =0  (4-34 ) 

th  ^ 

whicii  is  not  an  arbitrary  imposition.  Tlie  idea  is  that,  after  Y(X)  is 
correcti'd  to  include  H,  there  will  be  no  rigid-body  displacements,  and 
the  deformations  may  be  measured  relative  to  any  fixed  point,  here 
chosen  as  the  first. 

The  resulting  n equations  are  solved  foe  the  unknown  I’Y  using 

a subroutine  called  MVS,  wbich  is  one  of  a litirary  package  that  goes 
under  the  name  of  MIDAS  (Matrix,  Inverse,  Determinant,  and  Algebraic 
System).  As  stated  earlier,  two  cases  are  solved,  one  assuming  zero 
heating  in  the  cavitating  zones,  the  other  assuming  full-film  heating 


throughout . 
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I'cir 


Llif  I'lasLic  do f omia t i ons  , oq. 

1-f  1 1 1^1 

(DX)^ 


(4-50)  Kivi-s 


■>v 

I 


1”. 

_JL 


"k 

whcro  1 reprosonts  a summation,  with  the  l)oundary  conditions 


(4-35) 


‘ <5:’ 

1 n 


(4-36 ) 


and 


(4-37) 


Tile  same  pivotal  spacings  are  used  for  the  a.'s. 

In  the  actual  event,  tiic  search  for  the  threshold  of  instability 
( 1 + Zij.  1 = 1y1)  has  proved  to  he  time-consuming  and  expensive.  it 

has  also  been  found  that  Aj.  (X)  is  orders  of  magnitude  smaller  than 
Ath(X).  (hinsequent  ly , elastic  deformations  have  been  dropped  from 
consideration  and  stability  analyses  conducted  on  the  basis  of  the 
dominant  A , alone.  In  so  doing,  plots  liave  been  made  of  Zi  , / Y 

versus  u for  a number  of  discreet  values  of  e where 


The  resulting  curves  are  shown  in  fig.  (34).  For  each  value  of  e, 
the  critical  u is  given  by  the  point  where  the  corresponding  curve 
crosses  the  1a^I^|/|y1  = 1 ordinate.  These  critical  values  are  then 
plotted  on  a e versus  u grapli  as  shown  in  fig.  (35),  which  also  illus- 
trates the  difference  between  the  assumptions  of  full-film  heating  and 
zero-heating  in  the  cavitation  zones.  The  true  solutions  lie  in  the 
shaded  region. 

Appendix  I)  shows  the  program  used  in  the  computer  calculations. 


(Actual  Solutions  Lit*  in  Shaded  Regime 
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4.9  Nott's  cm  Various  Paraint'lors 

dni'  1)1  till'  most  important  paramotiTs  is  tlif  wavo  length . In  Chapter  11 
it  has  been  shown  that  the  largest  wavelengths  are  the  most  susceptible 
to  instability,  i.e.,  have  the  lowi’st  critical  speeds.  Tlie  results 
represented  by  fig.  (35)  are  based  on  the  assumption  of  a single  wave 
around  the  circumference,  i.e.,  a wavelength  equal  to  the  circumference. 
Inis  is  akin  to  a misalignment  between  two  perfectly  flat  faces  as  shown 
in  fig.  (36).  If  one  of  the  two  faces  is  flexibly  mounted,  as  is  the 
common  practice,  the  resultant  single-lobed  pressure  wave  will  tilt  that 
face  and  neutralize  the  initial  misalignment.  Thus,  a single-lobed  wave 
can  exist  only  when  both  faces  are  rigidly  mounted,  allowing  axial  but 
not  tilting  motion.  Wliere  a flexible  mounting  is  used,  a surface  wave 
must  have  at  least  two  lobes  in  order  that  hydrodynamic  pressure  may  be 
generated.  t:a  leu  lat ions  show  that  with  wavelength  equal  to  half  the 
c i rcumterence , the  critical  speeds  plotted  in  fig.  (35)  will  be  doubled. 

As  shown,  fig.  (35)  provides  the  critical  speeds  for  a rigidly  mounted 
seal  witli  a single-lobed  wave;  doubling  the  values  of  u will  describe 
the  threshold  of  instability  for  <a  flexibly  mounted  seal  with  a two- 
lobed  wave.  Lach  will  represent  the  lowest  critical  speed  obtainable 
for  the  given  mounting. 

An  analysis  was  made  and  a computer  progriun  used  to  allow  for 
tilting  of  tfie  laces  during  operation,  and  are  given  in  Appendix  K.  Use 
of  this  program,  however,  was  discontinued  subsequent  to  the  realization 
of  the  above  facts.  A symmetrical  two-lobed  wave  will  always  remain 
symmetrical  and  produce  no  tilt. 

To  test  whether  the  equations  for  the  wavy  surface  can  be  reduced 
to  those  derived  in  Chapter  11,  let 
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e» 


y.  £ 


■1 

aroic  Fi  Nutat 

eVative  Vac 

,y,  ^ 


tends  to 
ion  (Stanghan- 
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li  = 0 
() 


sinnifyinj'  a flat  surface.  Then  eq.  (4-10)  becomes 


d^6 


dx 


' 1 

(ir  + r,i  - 4 (h'  + h)dx 

1/  u ^ ^ ^ 


(4- J8) 


Then,  let 


li ' = A sin  •'.x 


whence 


J 


i h ' dx  = 0 


so  that 


a 2 

— 2 ^^th^  “ ’ ~ -2 
dx  ” ' 


A sin  <x 


Thus  , 


a 

'th  = “TT  "2  ^i’^  + S 


whe  re 


= 0 


and 


C^x  + = 0 


in  keeping  witli  tlie  boundary  conditions  described  by  eqs.  (4-33)  and 
(4-34).  Tfien,  in  order  that  6^^  may  be  equal  to  h',  one  must  have 

o2 


-2  2 


1 
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J LlO' 

^ M 


which  is  the  same  as  ttiat  obtained  in  Chapter  II. 

This  linearized  derivation  shows  that,  witli  initially  flat  surfaces, 

vertical  ri^id-hody  motions,  when  allowiui,  do  not  affect  the  critical  speed. 

'Ihis  minht  lead  one  to  believe  that  fi  has  minimal  or  no  effect.  However, 

numerical  solutions  show  that  if  fi  is  assumed  to  lie  zero,  V . Rets 

cr  1 1 

smaller  as  |h  j gets  larger,  while  the  opposite  occurs  if  fi  is  allowed 
to  ht  non-zero.  Tliis  would  seem  to  be  an  anomaly,  because  larger  jh^j's 
should  be  expected  to  lead  to  greater  disparities  in  q'(x)  along  the 
film  and,  consequently,  lower  critical  speeds.  A linearized  simplifi- 
cation of  the  heating  equation  suggests  an  explanation. 

If  h^  h,  eq . (^-b)  reduces  to 


, ,2  2h  h'  2h  fi 

q = - ' h'  + h ^ — 

" h h ^ 


(4-19) 


It  fol lows  that , 


for 


h 0 
o 


q'  = 


^ (h-  + fi) 
h 


(4-40) 


in  wtiich  case,  as  already  shown,  fi  does  not  affect  the  critical  speed. 

On  comparing  eqs.  (4-39)  and  (4-40),  one  might  guess  tliat,  with  non-zero 
h^^  and  fi,  the  last  term  in  the  parenthesis  of  eq.  (4-39)  is  the  one 
that  makes  all  tlie  difference. 

Recalling  eq.  (4-21),  and  defining 

Y = ^ sin  X 
h 


H = 1 + e cos  X , e < 1 
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,ind  St  ipu  1 ;)t  in>>,  further  , tliat  11  is  zeri>,  one  has 


2 2’’ 

^ ^th  u"A  ' sin  X 1 ''  sin  X I 

— = . — — ^ - — j dX  , 

dX  h (1  + t cos  X)  o (1  + e cos  X) 


(4-41) 


I'his  can  he  integrated  to  Kive  ^i7j 


u'a  , 


1 + e cos 


_ . -L  ^ ^ + c: 

X 2"  . 1 + e cos  X j()  1 


£A  i + 

l+€cosX  1 
he 


(4-42) 


Integrating  again,  one  obtains  L37J 


.l-c^ 


-1  ' e + cos  X 

cos  — ; — r + C.X  + t._ 

, 1 + e cos  X_,  1 2 


(4-43) 


5=1  for  sinX^O 


sin  X ^ 0 


To  determine  the  constants  of  integration  and  C^,  one  has 
use  the  boundary  conditions 


X = 0 . A^,^  = 0 


X = 2-^  . = 0 


which  give,  frimi  eq.  (4-43), 


0^  = 0 


he  V l-e 


(4-44) 


The  implication  of  eq.  (4-44)  is  that  the  maximum  and  minimum  values 


of  will  not  correspond  exactly  with  the  peak  and  the  trough  of  Y(X), 

2 

but  will  be  displaced  slightly,  corresponding  to  the  value  of  \ l-e  . 
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Kf  t urn  i iiK  to  i'<)s.  (4-41)  and  (4-42),  one  can  now  write 


th 


2 

n^'A 

he 


-i^  e + cos  X 

C OS 


\1  + e cos 


x; 


V l-e 


(4-45) 


th 


dX 


u^A 


7"  \1  + e cos  X 

he 


(4-46) 


l-e 


2' 


At  the  points  wliere  attains  its  maximum  and  minimum  values, 


one  has 


dii 


th 


dX 


= 0 


which,  from  eq,  (4-46),  requires 


cos  X = 


l-e*"  - 5 
e6 


frcjm  which  it  will  he  seen  readily  that,  of  the  two  values  of  6,  only 
the  positive  is  admissible  in  this  case.  Thus, 


and 


cos  X = 


X = cos 


1-e^  - 1 


■1  1 - V l-e‘ 


(4-47) 

(4-48) 


so  that  one  sets 


e + cos  X _ 1 - vT-e‘ 
1 + e cos  X e 


(4-49) 


UsinK  the  above  relationships  in  eq.  (4-45),  one  obtains 


> u^A  6 -1/1  - 1-e^ 

^ = ; cos  I 

he  --2  ^ ® 

..l-e 


J 


-lf_  1 - V 1- 


.2 


(4-50) 
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in  wliicli  llu-  inaximuin  corresponds  to  6 = -land  tl\f  minimum  to  5 


One  can  write,  tlicn. 


th' 


lu  , 1 = 1a  - a . 1 

'til  max  min 


2u^A  I * 1 / .j  ^ / a ,1 

cos  (y)  - cos  (-0)| 


he  \ 1-e 


where 


e = 


l-e‘ 


2A  I I 

and,  bv  definition,  ^ = 1y1 

h 

Thus , 


th' 


Icos  ^(9)  - cos  ^(-9)| 


e V 1-e 


and,  for  neutral  stability  = IyI) 


e ^ 1-e 


J/2 


‘ jcos  ^(9)  - cos  (-9)|^ 

It  will  be  seen  that,  for  small  e,  i.e.,  small  amplitudes  of  h^. 


e.f 


so  that  eq.  (^-52)  gives 


e .y  1-e 


1/2 


- il-' 


cos  I - cos 


V 2 


yi 


and,  as  illustrated  by  fig.  (37), 


- 1 e - 1 e - 1 , , 

cos  - cos  -y  • - e,  for  small  e 


(4-51) 


(4-52) 


. u 


(4-53) 


Figure  37.  Illustration  of  Trigonometric  Relationship  for  Arc  cos 
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It  i ■;  cli'iir  t r(im  oq.  (4-31)  that 

(i)  u = 1 for  a flat  surface  (t;  = ()) 

and  (ii)  u • 1 for  a wavy  surface  and  larger  values  t)f  e give 

smaller  values  of  u. 

Keturning  now  to  the  case  where  11  4 0,  i.e.,  relative  axial  motion 
prevail'-  between  the  two  rigidly  mounted  faces,  it  will  be  seen  that 
analytical  integration  of  eq.  (4-21)  liecomes  very  complicated.  Furtfier- 
more,  H can  he  evaluated  only  from  the  pressure  distribution  (see,  for 
example,  eq.  (4-19))  and  is  best  done  numerically.  .Since  the  purpose, 
in  tfiis  section,  is  to  illustrate  tfie  nature  of  the  dependence  of  u on 
H,  it  will  be  appropriate  and  sufficient  to  simplify  the  analysis  to 
include  only  vi’ry  small  values  of  e. 

Uef ining 

Y = ^ sin  X 
h 

11=  l + esinX 

and  using  eq.  (4-21),  one  notes  that,  for  € « 1, 

= 1 - 2e  sin  X 
II 

then  becomes 

= - ^ sin  X + ri;(l  - 2e  sin  X) 

'h 

2tt 

+ ^ sin  X + H (1  - 2e  sin  X)dX 

^ TT  t.'  \—  / 


Kquat ion  (4-21) 


2 2 
u dX 


108 


wlienci' 


and 


— — — = (-  - 2e[l  y cos  X 
2 tlX  / 

u h 


^ sin  2X  f <: 
2h 


-^r  j , = — - 2eH  sin  X + — cos  2X  + C.X  + (:_ 

2 til  s-  J .-r  12 

11  li  4h 


L'sin^  the  boundary  conditions 


X = (I  : A , = 0 

th 

X = 2'’  ; A = 0 

th 


one  finds 


eA 

4h 


and 


C 


1 


0 


(4-54) 


(4-55) 


At  the  points  where 


th 


is  maximum  and  minimum 


dA 


th 


dX 


0 


so  that  eq.  (4-54)  gives 

^ - 2eH j cos  X - ^ sin  X cos  X = 0 
'h  h 

wltence  either  cos  X = 0 

or  sin  X = — - 268, 

eA  V / 

h 

The  second  condition  is  inadmissible  because  it  yields  | sin  x|  > 1. 
The  first  condition  implies 
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tin  whicti,  usitij;  oq.  (4-53),  it.  is  sot-n  tiiat 


tl.' 


2,,  2t;iiii| 

— " 1 ' - I 


liy  definition, 


lY  = 


2A 


so  tliat,  for  neutral  stability  and  C 1, 


— 1/2 
^ 2eHh  ' 

^ ‘T~y 


(4-5b) 


Kquation  (4-56)  shows  clearly  tiiat  u is  always  greater  than  unity  wlien 
11  is  present.  Furthermore,  it  lias  been  numerically  determined  that 
llh/A  gets  larger  as  e increases.  It  follows  tiiat  larger  values  of  e 
will  result  in  liigher  values  of  u. 


CliAPIT'.K  V 


?RKOii;ill)N  OF  INSTABILITY  FROM  CIVF.N  INITIAL  CONDITIONS 
FOR  TUL  INITIAl.l.Y  WAVY  SURFACK 


j.l  i.^per.'it  1 n^;  and  Initial  Conditions 

It  lias  btHMi  statod  in  section  (4.1)  tliat  tlie  snrfaci'  wavincss  h (x) 

o 

IS  not  t lu'  initial  or  machinod  profilo,  but  ratlior  tlu-  profilo  observed 
at  any  Riven  moment  of  operation.  Together  with  the  face  load  f,  tlie 
slidiiiR  speed  F and  tlie  mean  f i Im- th  ickne  s s h,  h^^  describes  a given 
operating  condition,  i.e.,  an  assumed  quasi-static  state.  The  question 
aski'd  and  answered  in  the  preceding  chapter  was  wliether  this  state  was 
stable.  The  results  represented  by  fig.  (35)  descrilie  in  effect,  a 
series  of  neutrally  stable  operating  conditions  (though  load  lias  been 
ignored);  they  do  not  tell  the  designer  liow  to  avoid  designing  a seal 
that  will  operate  in  the  unstable  regime. 

In  order  to  be  able  to  predict  whether  a given  seal,  as  manufactured, 
will  be  stable  in  operation  under  a given  load  and  at  a given  speed, 
the  data  of  fig.  (35)  have  to  be  "translated"  so  that  they  show  the 
effect  of  the  initial  surface  profile.  It  will  be  recalled  that  h 

o 

represents  the  cumulative  effects  of  thermal  distortion  on  the  initial 

waviness  h^.  Thus,  if  one  calculates  the  thermally  produced  component 

of  h (x)  atid  subtracts  it  from  h (x),  one  will  have  the  initial  profile 
o o 

hj^(x)  which,  together  with  the  specified  face  load  f and  the  operating 
speed  i; , will  represent  the  initial  conditions. 


no 


Ill 


3.2  Dt-tcnninaLion  of  Initial  Surface  Profile,  l.oad  and  SLicHn>>  Spued 
In  till-  assumed  quasi -static  state,  ttie  lieat  K‘‘n*'ration  is  given 
by  eq.  (4-5) 


= _t 


4' 


h + h 


which  corres[)onds  to  a thermally  generated  surface  curvature 


2 Ji 

d 6 O'  q , , O'  q 

th  M^o  1 K.  c 


dx 


i . 


o 'S, 


dx 


(5-1) 


Written  in  terms  of  non-dimensional  quantities,  eq.  (5-1)  becomes 

2 2rT 

d^A^,  2 , / 2 

th  _ u 1 u 

,y2  H ■ 2n  . H 

dX  o 


(5-2) 


and  can  lie  solved  numerically  as  outlined  in  section  (4.8),  using  the 
boundary  conditions  (4-33)  and  (4-34).  Recalling  that 


H(X)  = 1 + 


h^(X) 


it  is  clear  that 


H.(X)  = H(X)  - \j/X)  - 1 


whe  re 


H.(X)  = 


li^(X) 


To  determine  the  load  f,  it  will  be  noted  tliat 


(5-3) 


df  = (x ) • L- dx 


= p L -— 

*^0  K 


J 
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Ki'calliii)^  eq.  (4-22),  one  lias 


f 


^ I — y P dX 

liK  ' ^'m  ■ o 


2n 

r» 

winner  F =:  , |>  (X)dX 

V o 
o 


wluTi*  the  non-dimensional  load  is  det'ined  as 


fM 

' w.  • 1-1^/ 


f 


(5-4) 


(5-5) 


First,  P (X)  is  determined  numerically  from  eq.  (4-23).  Then  K 
o 

is  numerically  calculated  from  eq.  (5-4),  the  intej'ration  heinx  carried 
out  only  over  the  non-cavitating  regions  of  the  fluid  film.  The 
dimensional  load  f is  then  obtained  fri>m  eq.  (5-5).  The  program  is 
listed  in  Appendix  K. 


5.3  Results 

In  fig.  (38)  are  shown  graphs  of  u versus  log  K for  different  values 
of  e.  Kach  point  on  each  of  the  curves  represents  a possible  operating 
condition.  The  threshold  of  instability  is  the  locus  of  all  operating 
conditions  that  are  neutrally  stable  and  is  drawn  from  fig.  (35).  It 
is  seen  that  while  initially  the  non-dimensional  critical  speed  increases 
with  increasing  c and  F,  it  decreases  sharply  beyond  e 0.7  and 
log  F 1-3  2.1. 

On  fig.  (35),  one  can  now  draw  a number  of  constant  load  lines  as 
shown  in  fig.  (39).  Using  the  computed  values  of  and  figs.  (38)  and 
(39),  a comprehensive  map  can  be  drawn  (fig.  40)  showing  the  effects  of  all 
the  parameters  and  their  relationship  with  the  threshold  of  instability. 


Solutions  Lie  in  Shaded  Regime') 


Filiure  39.  Critical  Conditions  Defined  bv  Non-Dimensional  ('pirating  Snriace  Wavines 
Sliding  Speed  and  Load.  (Actual  Solutions  Lie  in  SVaded  Rigin'e) 


Sinc«'  tin-  piirpdst'  is  to  do  f i nt'  tho  critical  Llircsliold  in  tcrnis  of 
a nivi'ii  It^ad,  oporatinn  speed  and  initial  surface  waviness,  the  plots 
ol  tit"  t40»  have  to  he  renderi'd  in  dimensional  terms.  It  will  then 
!h-  possible  to  predict,  for  a ><iven  load  and  initial  surface  waviness, 
tlie  critical  speed  as  well  as  the  f i Im- th  ickiK’ ss  and  surface  waviness 
corresponding  to  the  critical  state.  This  has  been  done  in  two  ways, 
ifie  first  consists  of  writing  equations  for  the  critical  speed  in  terms 
of  tile  parameters  K,  1 11.  1 , e and  then  solving  them  simultaneously  to 
eliminate  h.  The  other  is  a sequence  of  direct  crossplotting  from  the 
data  of  fig.  (40).  In  the  latter  approach,  a value  of  |h,l  is  assumed 
and  values  of  h are  calculated  for  a number  of  points  on  the  threshold 
of  instability.  For  each  of  tliese  points,  is  then  readily  determined, 

whence  f and  U can  be  found  by  using  eq.  (5-5)  and  the  relationship 


u 


L' 


h^ 


• ptv 


M 


The  results  can  tiien  be  plotted  to  show,  for  a given  initial  amplitude 
IhJ,  the  critical  values  of  L'  corresponding  to  given  valui'S  of  f,  as 
in  fig,  (41).  It  is  also  possible  to  draw  graphs  (figs.  42  and  4.1)  showing  h 
and  [h^l  versus  load  corresponding  to  the  neutral  states,  for  a given 
jh^l.  In  order  to  be  able  to  predict  the  onset  of  instability,  liowever, 
the  first  approach  siiould  be  taken.  The  resulting  empirical  equations, 
of  course,  must  generate  the  graphs  of  figs.  (41),  (42)  and  (43). 

It  is  seen  from  fig.  (40)  tliat  the  tlireshold  of  instability  con- 
sists of  two  legs,  one  below  F j 2.1,  the  other  above.  The  two  legs 
are  represented  by  two  sets  of  equations  in  terms  of  nondimc-nsional 


less  for  a Given  Load  and  Initial  Surface  Waviness 


OJ 


P ^ CO  ro  O)  O'}; 

rO  rO  fO  ^ »r>  itI 
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Figure  43.  Critical  Operating  Surface  Waviness  for  a Given  Load  and  Initial  Surface  Waviness 
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variables.  Kor  tlie  first  or  leftward  Leg,  it  is  found  that,  assuming 
no  lieating  occurs  in  tlie  cavitating  regions. 


(5-6) 


i'  13.33 


u = 1.96(1 11^1  + 1.5) 


.286 


(5-7) 


111  J = 0.  158  e 


5 .66 


(5-8) 


These  relationships  are  empirical,  and  the  properties  and  dimensions 
listed  in  Appendix  C will  be  used  in  reducing  them  to  dimensional  form. 
Equation  (5-6)  gives 


i' 


1 h<  / M 


- 13.33  WL  '..K^ 


whence  , 


-2  - 5 ^ i ^ 

h = 1.  1454  X 10 


(5-9) 


From  Eq.  (5-7),  one  has 


r,. 


- h. 


.286 


I'  = 1.96  h<  — - i — ^ + 1.5 
tjPr..  r ' 


or. 


M h 


.Ih.  1 .286 

II  = 4.691  X 10  li  ( — — + 1.5  ! 

' h 


(5-10) 


From  eqs.  (5-9)  and  (5-10), 


4 / I I 2 1 

II  = 1.588  X 10  / I 2.955  x 10  . jh^lf  + 1.5y 


.286 


(5-11) 


Equation  (5-8)  may  be  written  as 

1,,J  ,^Xni-L 


. 158h 
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.iiul,  suhst  i til  L i UK  fur  li  from  oq.  (5-9),  one  lias 


1 1 = 6.04)5  X 10“^  -7^  ^1.87  x LO’^  . '1  h . 1 f y 

O ' t 1 / 


(5-12) 


Ihus,  knowing  i*'jl  dimensions  and  properties  (sucli 

iis  those  given  in  Appendix  C)  one  can  determine  the  mean  f i Im-thickness , 
the  sliding  speed  and  the  amplitude  of  the  surface  waviness  at  the 
critical  condition  by  using  eqs.  (5-9),  (5-11)  and  (5-12). 

For  tlie  case  where  heat  generation  is  calculated  on  the  basis  of 
a full  lubricating  film  (the  pressure  wave  is  still  assumed  to  be  com- 
prised solely  of  positive  values),  the  above  approach  leads  to  the 
following  relationships  (see  figs.  (38),  (39),  (44)) 


i'  8.51 


(5-13) 


u = 1.246(lnJ  + 2.2)' 


(5-14) 


|.C|  = i-L  <e‘-»h 


(5-15) 


liquation  (5-13)  gives 


h'  = H.51^.''- 

< ' O' 


In' 


or,  using  the  properties  and  dimensions  of  Appendix  C, 


h = 2.7041  X 10 


-3  / '"i 


(5-16) 


U = 1.246  M ,i— 
po- 


K |h. 1 . .378 

^ 2.2) 

»..  r / 


M h 


From  eq . (5-  14  ) , 


12'J 


or , 


V = 2.9H2  X 10  h ^ +2.2 


wlu'tu'o,  usiiiK  oq.  (5-l0),  oiik’  obtains 


r = b.obib  X 10 


i '"i 


i.b9H  X 10"^  '111.  1 f'  + 2.2 


i7H 


(5-17i 


inally,  usinH  oqs.  (5-15)  and  (5-lb),  it  is  found  tliat 


, Ih.l  . 

Ih  I = 4.006  X 10  in  3.291  x 10  , [hji 

o f 1 1 


(5-18) 


Fistiros  (42),  (43),  (45)  shows  comparisons  between  the  above  relationships 
and  corresponding  curves  obtained  through  direct  cross-plotting. 

()btaining  similar  equation.s  for  the  second  leg  of  the  thresliold 
(fig.  40)  proves  to  be  a little  more  difficult.  Wlien  heating  is 

assumed  to  be  zero  in  the  cavitating  portions  of  the  film,  it  is  found 

tfiat,  for  all  practical  purposes,  the  sliding  speed  U,  the  mean  film 
tliickm'ss  h and  the  operating  surface  wave  amplitude  ^te  independent 

of  load  For  all  loads  greater  than  that  at  wliicii  tlie  two  legs  of  the 

tlireshold  meet.  For  this  reason,  and  also  because  of  tlie  complicated 
equations  tliat  are  obtained  for  the  second  leg,  these  equations  are 
not  lieing  presented.  They  are  neither  informative,  nor  useable. 

Howk'ver,  figs.  (41  ),  (42),  (4  3)  sliow  the  second  leg  as  obtained  by  direct 
cross -p lot  t i ng. 


5.4  Remarks 

The  results  presented  in  this  chapter  apply  to  "floating"  but 
non- flexible  mountings,  l.e.,  tlie  mean  film-thickness  can  vary  but 


T 


LOG  U (ft/sec) 


f (lbs.) 

45.  (Comparison  Bc-twi-en  (Calculated  Critical  Speeds  and  Corre- 
sponding empirical  equations 
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in’ithiT  fact.'  can  tilt.  OtluT  npcratinn  parameters  Iteing  identic.al,  .a 
flexit>ly  mounted  face  must  necessarily  have  at  least  a two- lohed  profile 
and  the  critical  sliding  speeds  will  then  he  double  those  represented 
by  fig.  (A  1 ) . 

(ine  might  expect,  intuitively,  tliat  the  critical  speeds  will  he 
greater  for  the  case  when  full-film  lieating  is  assumed  than  for  the 
case  of  zero  heating  in  the  cavitating  regions,  I'ecause  the  disparities 
in  heating  are  mori-  abrupt  in  tin.  latter  case  than  in  the  former,  as 
illustrated  in  fig.  (Ah).  However,  figs.  fJ9)  and  (Al)  belie  such 
expectations.  The  explanation  is  not  apparent.  At  any  rate,  as  has 
been  stated  earlier,  the  true  solutions  lie  in  between  the  two  cases, 
in  the  shaded  regimes  shown  in  the  figures. 

Finally,  it  should  be  noted  that,  in  obtaining  the  present  results, 
a number  of  "operating  conditions"  wc're  selected.  Wliile  such  operating 
conditions  doubtlessly  represent  mathematically  admissible  quasi-static 
solutions,  it  cannot  be  said  yet  whether  they  do  occur  in  practice, 
i.e.,  whetlier  a given  "initial  condition"  can  evolve  into  its  corre- 
sponding neutrally  stable  "operating  condition".  Tlie  following  chapter 
seeks  to  answer  that  question. 


^ — V \ k — ^ — r — "T  r: 

a 


b 

full  film 

CAVITATING  FILM 

C CAVITY 


KiK”'"*-'  Meat  Generation  Based  on  Full  Film  and  Gavitating  Film 
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CMAITKR  VI 


KVOUTION  OF  OPKRAI'ING  CONDITIONS  FROM  GIVKN  INITIAI,  CONDITIONS 


I II  Lho  last  chapter,  it  lias  liccn  shown  tliat  , for  a I'ivun  initial 
wavini'ss  and  load,  a neutrally  stable  operatin>;  cc^ndition  can  he  defined 
in  terms  of  sliding  speed,  mean  f i lin- 1 h ickni' ss  and  surface  waviness.  To 
determine  whether  sucli  an  operating  condition  can  actually  he  attained 
in  practice,  one  needs  to  follow  the  growth  of  the  thermal  deformation 
as  speed  is  increased  from  rest.  This  involves  the  integration  of 
eq.  (4-8)  which,  written  in  terms  of  the  zero-order  quantities,  is 


d^  _ %%  _ 

dx^  S 

where  v denotes  tlie  thermal  deformation, 
malizing  procedure  of  section  (4.7),  eq. 


(6-1) 

Using  eq.  (4-5)  and  the  nor- 
(6-1)  can  be  written  as 


where 


2, 

dX 


2 

u 

II 


C 


C 


dX 


(6-2) 


o 

If  one  knows  H and  u,  one  can,  presumably,  determine  \’  from  eq.  (b-2). 
However,  this  equation  is  difficult  to  integrate  when  11  is  the  independent 
variable.  This  difficulty  can  be  circumvented  by  choosing  V as  the 
independent  variable  and  then  asking  what  the  corresponding  11  must  be. 
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1.0  t 


V = V CD.s  X 
o 


Thon 


V cos  X 
o 


so  tliaL  (“q.  (h-2)  can  l>o  written  as 


whence  , 


2 

C - V cos  X = 

o II 


II 


u"  ''  1 • 

C V 1 

1 - ^ cos  X 


(6--3) 


Now,  by  definition  (sec  fig.  27), 


L 

h = -j  h dx 
o 


or,  in  non-dimensional  terms  (fig.  32), 


1 


1_ 

2v 


2n 

H dX 

c 

() 


(6-4) 


Substitution  of  H from  eq.  (6-3)  into  eq.  (6-4)  gives 

. 2tt 


1 = 


dX 


u __ 

2nc  c V 

o , o V 
I - — cos  X 


and,  following  integration 

1 = 


/ /V  2 

r ^ -2. ' 

^ - V,C  ) 


c = + 

\ o 


1/2 


whence 


(6-5) 
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From  eq.  (b-  5),  it  is  set-n  that 

dX 


u V sin  X 
o 

! ^2 
C - V cos  X I 

r\  ! 


(h-6) 


The  amplitude  e of  the  surface  wave  will  he  the  difference  between  the 
maximum  and  minimum  values  of  11,  for  which 


= 0 
dX 


Kquation  (h-h)  then  shows  that,  for  11  to  have  its  maximum  and  minimum 
values, 

sin  X = 0 

which  implies  cos  X = - 1 

Then,  in  keeping  with  the  definition  of  e (see  fig.  32), 

e = 4 I - » • 1 

2 max  min 

1 2 2 
1 I u U I 


2 'C  - V C + V ' 
o o 


V u 

1_2 1 

I 2 2 ' 

C - V 

o 


whence,  using  eq.  (6-5),  it  follows  that 

V 

€ 


(6-7) 


Rewriting  eq.  (5-3), 


= H - V - 1 


or. 


T-T, ^ - V cos  X 

C - V cos  X o 


(6-8) 


lio 


from  wliich 


clU. 


u V sin  X 
o 


(r  w V 

' C - V cos  X 
o 


- + V sin  X 
2 o 


(6-9) 


I t 


dll . 
1 

dX 


= 0 


equation  ((i-9)  >iivcs  either 


or 


sin  X = 0 


u = - ((;  - V cos  X) 
o 


(6-10) 

(6-11) 


Krom  eq.  (h-11),  it  lollows  that 

|nj  = 2u 

wtiich  cannot  be  correct  because  it  calls  for  an  unique  u for  a given 
jn^j.  Kquation  (6-10)  implies  that 


2 ' C - V 


. 1 - V 


o ^ ' "oi 


which  reduces  to 


^ - ilv. 


C - V 


and,  upon  substituting  for  C from  eq.  (6-5), 

1 - 


1 - 


In  the  last  equation,  ( j is,  by  definition,  the  amplitude  of  the 


initial  waviness  of  the  surface.  The  second  modulus  symbol,  however, 


i 

I 


I 


1 31 


w.iN  int  rociiitfd  Lo  around  tiu-  qut“;tion  ol  wliicti  of  the  tw<j  values 

of  H.,  corresponding  to  zero  slope  (eq.  (6-10)),  was  tfie  maximum  and 
wfiicli  the  minimum.  This  modulus  symbol  sfiould  he  dropped  now  to  allow 
to  have  positive  and  negative  values,  which  imply  opposite  phases. 
Thus 


1 


(6-12  ) 


I ** 

or.  111.  1 = V (6-13) 

u" 

wtiich  shows  that  IL  and  V will  be  in  or  out  of  phase  accordingly  as  u 
is  less  than  or  greater  than  unity. 

Combining  eqs.  (6-7)  and  (6-12),  one  obtains 

lllj 

e = — (6-14) 

1 - u 

which  shows  how  e will  evolve  with  changing  u,  for  a given  lll.l. 
equation  (6-14)  also  shows  that  if  u is  equal  to  unity,  i.e.,  the  sliding 
speed  is  the  critical  defined  by  eq . (2-84),  e will  be  infinite  for 
finite  initial  waviness.  For  u < 1,  € will  be  positive.  For  u ^ 1, 
e will  be  negative.  This  means  that,  as  u increases,  the  surface 
waviness  e increases  wtiile  remaining  in  phase  with  IC(X),  until  u 
becomes  unity.  For  u > 1,  H^(^)  suddenly  goes  out  of  phase  with  11^  (X) 
and  then  t gets  smaller  as  u gets  larger.  If  £ does,  indeed,  become 
infinite  when  u reaches  unity,  it  might  mean  that  the  two  surfaces  will 


touch  and  the  analyses  of  Chapter  IV  will  break  down.  However,  such  a 
conclusion  cannot  be  drawn  from  eq.  (6-14)  because  all  the  variables  in 
it  are  normalized  with  respect  to  h,  and  the  behavior  of  a face  seal  in 
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practice  will  depend  upon  changes  in  h.  To  determine  whether  or  not  the 
system  can  actually  attain  unit  u,  and  how  it  behaves  If  it  does,  eq , 
(h-14)  will  have  to  be  rendered  in  terms  of  dimensional  variables.  Then 
Keynold's  equation  will  have  to  be  employed  to  determine  pressure  vari- 
ations and  consi'quent  variations  in  the  mean  f i Im- th  i ckne  s s . 

Written  in  terms  of  dimensional  quantities,  eq.  (b-l't)  becomes 


It'-l 


I - 


(b-15  ) 


where  it  should  be  remembered  that  the  modulus  symbcjl  represents  an 
amplitude  but  does  not  preclude  a sign,  positive  or  negative. 

As  shown  in  section  (5.2),  the  load  f is  given  by 


df  = p h dx 

‘ f\ 


which,  integrated  over  the  full-film  domain,  becomes  (fig.  47) 


,"l  ^ 2 


f = 


p L dx 


where  \ is  the  wavelength.  Then,  recalling  p^  from  eq.  (4-14),  one  has 


f = - 


\ 

3 /l  2 , d(h  + h ) 

P-lih  I o_ 

2 . /r  j I \ ^ 

X,  (h  + h ) 

i o 


dx 


li  there  is  only  one  wave  around  the  seal  circumference,  one  can  write 
(see  ilg.  47) 


h - h 


ML  ' 
2 


o'  d(h  + h ) 
o_ 

h + Ih  1 


b 


KiKure  47.  Definition  Sketch  for  Calculation  of  Load 

a)  Operating  Surface  Profile 

b)  Corresponding  Pressure  Wave 
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which  rcducfs  tn 
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i atul  (fi-16)  can  be  combint'd  to  ^ive,  after  st)me 

rea rran^in^' , 


') 


i*'- 


I 2- 


■> 

U">jn 


1 - 


M 


-2  2 


/'  1 h ,•  I = 


(h-  Ih) 


(h-  17  ) 


Kor  giv»'n  initial  waviness  Ih^j  and  load  f,  eq.  (6-17)  ttives  h at  any 
sliding  speed  U . ilien  the  correspt)nding  surface  waviness  can  be 

ditemiined  from  eq . (6-16).  Kquation  (6-17),  of  course,  Is  hardly 

tractable  by  analytical  means  and  is  numerically  solved  by  usin^  trial 
values  of  h and  determining;  tfiose  for  which  the  left  hand  side  of 
eq.  (6-17)  reduces  to  zero.  Four  real  solutions  are  obtained,  of  which 
two  are  untenable  [>ecause,  when  substituted  into  eq.  (6-lh),  they  give 
|h^|  • h.  Thus,  at  least  in  theory,  there  are  two  possible  combinations 

of  h and  Ih^l  for  most  sliding  speeds. 

Figure  (48)  shows  how,  for  a given  initial  waviness,  the  mean  film- 
thickness  varies  with  increasing  sliding  speed.  It  is  seen  that  film- 
thickness  increases  rapidly  with  U and  then  asymptotically  approaches 
the  straight  line  that  represents  eq.  (2-84).  So  long  as  the  film- 
thickness  evolves  along  the  lower  curve,  the  system  will  not  attain  a 
neutrally  stable  state  described  in  Chapter  V.  However,  since  the 
low«  r curve  hugs  the  straight  line,  instability,  if  and  wfien  it  occurs, 
must  always  fall  on  the  line  representing  family  IV  in  fig.  (15). 

This  might  explain  the  closeness  between  the  predictions  of  Chapter  II 
and  experimental  results,  as  shown  in  fig.  (23). 


In  tii;.  (nH),  the  uppiT  .T^NTnptntic  cnrvf  rcprc  si-nt  s the  second  set 
ol  solutions  to  eq.  (t)-17).  On  tliis  curve  lies  the  neutrally  stable 
point  described  in  Cliapter  V.  Instal)ility  prevails  to  the  left  of 
tills  point.  fhis  proves  that  tlie  lunitrally  stable  ojierating  conditions 
described  in  tdiapti^r  IV  are  mathematically  olit  a i nai>  le  . it  is  unlikely, 

howiver,  that  thev  will  be  attained  in  practice.  The  system  will  evolve 
alonn  the  lower  curve.  II  "jolted”,  it  mij;ht  jump  over  to  the  upper 
curve,  but  it  will  still  be  stalile  and  will  pass  through  a succession 
of  states  vi‘ry  close  to  those  represented  by  eq.  (2-H4;. 

Figure  (4H ) shows  the  curves  for  three  different  values  of  initial 
waviness,  but  the  same  load  (10  lbs.).  Load  is  important  only  for  the 
early  part  of  each  curve,  before  tlie  curve  approaches  the  straight  line. 
Ihe  four  experimental  data  points  from  fig.  (23)  are  also  shown;  during 
the  experiments,  however,  initial  waviness  could  not  be  measured. 

The  evolution  of  tiie  surface  waviness  |li^|  is  siiown  in  fig.  (49). 
Ihese  curves  are  obtained  from  eq.  (6-15)  by  using  the  values  of  h corre- 
sponding to  the  curves  of  fig.  (48).  The  neutrally  stable  points  pre- 
dicted in  chapter  V (figs.  41  and  43)  for  a face  load  of  10  pounds  are 
•also  shown,  as  are  the  points  obtained  in  the  experiments  of  Chapter  111, 
immediately  preceding  the  onset  of  instability,  wtiere  the  load  was 

6.25  pounds.  The  latter  suggest  tliat  the  polished  aluminum  face  in  the 

-4  -5 

experiments  had  an  initial  wave  amplitude  of  the  order  of  10  to  10  in. 
Figure  149)  also  shows  how,  with  increasing  sliding  speed,  the  surface 
waviness  will  increase  along  the  rightward  curves  and  never  reach  the 
neutrally  stable  operating  conditions  envisaged  in  Chapter  V. 
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riu-  experimental  points  in  fig.  (49)  lend  added  credibility  to  the 
analysis.  The  metal  face  was  lie  Id  by  hand  and  polished  on  a wheel.  The 
inferred  magnitudes  of  the  initial  waviness  seem  very  realistic. 


CflAPTCH  VII 


CONCLUSIONS 


Wlion  a thermally  conducting  body  slides  on  a thermal  insulator, 
thermoelastic  instability  is  found  to  occur  under  certain  conditions. 
Small  disturbances  on  the  conducting  surface  grow  and  the  surface  under- 
goes macroscopic  deformation. 

If  both  surfaces  are  perfectly  flat  and  separated  by  a film  of 
viscous  liquid,  the  regime  of  instability  is  determined  by  the  mean 
film-thickness  and  the  widtii  of  tlie  sliding  interface.  Wlu-n  tlie  width 
is  large  enough  for  cross-flow  to  be  neglected,  the  instability  will 
not  occur  at  any  sliding  speed  of  practical  interest.  The  sliding 
speed  at  which  neutral  stability  prevails  is  called  the  critical  and 
is  given  by 

The  instability  occurs  for  only  a range  of  h,  and  at  sliding  speeds 
greater  than  U 

cr  It 

Kxperiments  conducted  with  a polished  aluminum  face  sliding  on  a 

poll  shed, f lexib ly  mounted,  glass  surface,  where  the  surfaces  can  move 

apart  from  each  other,  show  that  the  aluminum  face  develops  a large 

waviness  wlilch  grows  very  dramatically  at  a certain  sliding  speed.  i'his 

speed  is  found  to  be  very  close  to  11  . . The  sudden,  dramatic,  growth 

j crit 

of  the  surface  undulation  into  spikes  or  asperities  is  interpreted  as 
the  occurrence  of  thermoelastic  instability.  These  spikes  break  through 


1J9 
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thf  liihruaiit  I i Im  .-iiKi  rub  directly  against  the  gla^;'..  On  slowing  down 
t roin  the  unslabli’  state,  the  metal  face  is  seen  to  rt'gain  its  original 
pro  tile. 

It  the  conducting  lace  is  wavy  at  the  very  outsi't,  Ll>eoretical 

eonsidt  rat  ions  show  that  thermoelastic  instal)  i 1 i t y , defined  as  the  growth 

ot  a small  disturbance,  can  occur  only  at  sliding  spei'ds  consideraiily 

higlu'r  than  r . The  initial  surface  profile  changes  wi  tli  sliding 
c r 1 1 

'•peed  and  grows  continuously  through  successive  operating  conditions. 

If  both  faces  are  rigidly  mounted,  a larger  operating  waviness  corre- 
sponds to  a lowi'r  critical  speed.  If  the  faces  can  move  normal  to  each 
otlier,  the  growtli  of  tlie  surface  waviness  causes  such  motion  and  a larger 
operating  waviness  gives  a higher  critical  speed.  The  critical  speed 
in  this  case  depends  on  the  face  load,  the  face  width,  and  the  amplitude 
of  the  initial  waviness.  Viscous  heating  is  calculated  with  the  film 
considered  full  throughout  and  also  with  the  film  assumed  to  be  non- 
existent in  the  cavitating  regions.  The  general  thermoelastic  heliavior 
is  seen  to  be  the  same  in  the  two  cases,  the  difference  being  in  the 
magnitudes  of  the  critical  sliding  speed. 

An  analysis  of  the  continuous  thermal  distortion  of  the  conducting 

face  shows  tliat,  as  sliding  speed  increases,  the  amplitvide  of  the  faci- 

profile  also  increases.  The  mean  film-thickness  increases  too,  and 

the  sliding  speed  is  then  found  to  be  very  close  to  U . , liut  n^ver 

cr  It 

greater.  The  surface  waviness  and  mean  film-thickness  thu-.  grow  with 
the  sliding  speed  and  asymptotically  with  the  threshold  defined  by 
For  the  given  initial  waviness  and  load,  the  neutrally  stable  operating 
conditions  mentioned  in  the  preceding  paragraph  lie  across  this  threshold 
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•ind  an-  novt-r  attained  in  praitiie,  altlioiinh  they  are  matiiemat  ica  I ly 
admissible  qvias  i - s tat  i c solutions.  Tliis  means  that  if  and  when  insta- 
bility does  occur,  it  occurs  at  a sliding  speed  very  close  to  IJ 
As  to  tiu'  precise  mechanism  that  triggers  the  instability,  several 
possibilities  i-xist.  These  possibilities  are  not  discussed  because, 
without  supportive  analyses,  they  are  speculative  in  the  main,  and 
such  analyses  have  not  lieen  conducted. 

rhe  last  statement  notes  one  of  a few  questions  that  have  been 
raised  but  not  answiTed.  It  is  worthy  of  note,  however,  that  tlie 
different  facets  of  the  present  work  are  mutually  supportive.  The  flat 
surface  and  initially  wavy  surface  analyses  and  the  experlmantal  data 
fit  togetlier  vc-ry  well  to  form  a comprehensive  picture. 


APPKNDIX  A 


IHP;  VAI.IDITY  OF  A PAKAliOMC  TKMPKRATURi;  U I STR  1 RUT  I ON 


In  section  (2.(>),  it  is  stated  that  the  tempi-rature  profile  in  the 
oil  film  is  approximated  hy  a polynomial  of  the  second  degree.  ft  is 
worthwhile  to  assess  the  departure  ol  this  approximation  from  the 
correct  temperature  distribution. 

I lu'  energy  equation  for  an  incompressible  flow  in  a thin  film  is 
as  written  in  eq.  (2-26)  where  T is  the  temperature  perturbation.  In 
writing  eq.  (2-26i  account  was  taken  of  the  fact  that,  for  a narrow- 
taci'  seal  ot  the  kind  studied  here  and  for  a small  perturbation, 

2 2 2 

i-L  111 

2 ''  2 ’ 2 
3y  5x  Sz 

u = u^  + 11^ 


where 


« 


Ie  it 

Yx  ^z 


w 


/ 


2 


^ i ~ 0(2)  in  the  perturbation  , 

where  U2  refers  to  the  Couette  component  and  u to  the  pressure-generated 
component  of  the  composite  flow.  It  should  be  noted  in  particular  that, 
for  the  narrow- lip  configuration  and  a small  perturbation,  the  contri- 
bution  of  u^  (i.e.,  ^)  to  the  frictional  heating  term  is  so  small  that 
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it  tan  bi-  discarded  from  the  energy  equation.  Equation  (2-2t)}  may  be 
revn  itten  in  d i f f rcn t ia  1 form: 


dT  3T'' 

\¥r  ^ ^2 


^ ^ hi 

«“ay'  h2  -- 


J 


(A-1) 


Kecogni:'ing  that  in  this  problem 


dT  dT 

At  Ax 


this  term  in  eq.  (A-1)  will  be  small  whenever  — is  small.  Then  eq.  (A-1) 

u 


mav  be  rewritten 


.r  He  - M1  fhl 

^^^2  Ax  ■ ^oil  . 2 “ -2  - J 

0 y h h 


(A- 2) 


In  most  lubrication  studies,  tlie  convection  term  on  the  left  hand 
side  is  left  out  because  it  is  usually  very  small  compared  to  tlie  con- 
duction term.  In  that  case  the  temperature  profile  will  be  precisely 
a second-degree  polynomial,  as  assumed. 

If  convection  is  non-zero  but  still  small  relative  to  conduction, 
and  one  wishes  to  allow  for  it,  one  may  write 

T = T + T' 

where  T is  the  parabolic  temperature  profile  that  corresponds  precisely 
to  the  case  of  zero  convection  and  T ' is  a correction  on  it  to  account 
for  the  small  convection.  The  energy  equation  (A-2)  may  then  be  written 


V -1-  ("?■  + T'  I - K fT  4-  T'  ) = ) 


- ><011-2  ('>■+■>''> 


>1 

Ay' 


,2  / 

h h 


(A-3) 


wtu‘re,  f)y  definition. 
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li  h 


Kqiiatioii  (A-U,  written  in  tlie  int<‘i;ral  f orm,  becumes 


Pi'l'  *■  ■"'T  31’  . 

^ — V ^ i I civ  - K . 

— . -•  X 3 X . oil.  _ 

h V,  y, 


2 ^2  . 

^ n + T' ) 


= 


or  , 


h 


5T' 

dx 


dv 


ni  I 


■ 

^y  J 

V 


= 


S i nee 


3T'  ^ 

^x 


(A-4) 


(A-5  ) 


and,  as  long  as  convection  is  much  smaller  than  conduction, 


f«  yi 


^3T’  ' 
oil  I ^ 


Then,  using  the  boundary  condition  given  by  eq.  (2-23),  one  can  write 


J2 

m , 


ST 


ST 


Sx  ''g  *^oil  Sy 


(A-b) 


y=y. 


It  slujuld  be  noted  that  in  e 
perature  distribution;  so  also  is 


q.  (A-6) 


q is  independent  of  tlie  tem- 

g 

which  is  provided  by  an 


independent  boundary  condition  (eq.  2-22).  Thus  the  only  term  dependent 
on  the  tempc'rature  profile  is  the  convection  term  and,  according  to  eq. 
(A-b),  the  parabolic  profile  of  T provides  the  appropriate  evaluation  of 
convection  when  T'  ■s<  T.  Stated  differently,  given  q^,  a parabolic 
temperature  distribution  in  eq.  (A-6)  will  provide  the  correct  tempera- 


ture gradient  at  the  conducting  surface. 
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I 111-  ahuvc  analysis  cannot  l>c  cxLenclocl  directly  to  tlu>  ease  where 
convection  is  comparable  to  conduction.  In  tliat  case,  the  correct  tem- 
perature profile  will  contain  terms  in  the  third  and  hiRher  powers  of  y. 
One  m.iv  still  replace  the  correct  profile  with  a paral>olic  one  that 
mi'et  s the  -anie  boundary  conditions  and  also  satisfies  eq.  (2-20).  The 
saciitice  in  tfte  shape  ol  the  temperature  profile  will  induce  a small 
errot  in  the  iiuiRn  i tude  of  the  temperature  and,  in  the  inteRrated  enerRy 


e(]uai  ion , 


and  will  remain  uncliaiiRed  while  the  convection  term. 


intCRrated  across  thi'  film,  will  undergo  a small  change.  One  set  of 
calculations  with  sizable  convection  indicated  that  the  error  induced 
by  the  paiabolic  approximation  is  less  tlian  307„  of  T'. 

An  order  of  magnitude  analysis  shows  that,  for  the  range  of  family 
IV  in  fig.  (15),  the  ratio  of  convection  to  conduction  is  of  the  order 
of  10  for  h = 1.20  X 10  ^ in.  and  10  ^ for  h = 3.4  x 10  ^ in.  At  the 
right  hand  I'xtremity  of  this  range,  convection  becomes  important,  as 
evidenced  by  eqs.  (2-72),  (2-73).  From  the  foregoing  arguments  it  would 
appear  that  for  almost  the  entire  range  of  the  solution  represented 
by  f.imily  IV,  which  is  tlie  threshold  of  principal  interest,  the  parabolic 
approximation  is  quite  appropriate. 

Georgopoulos  L34J  has  rigorously  solved  the  energy  equation  and 
shown  that  convection  has  little  effect  on  the  temperature  at  the  metal 


surface . 


ai>pi;ni)[X  li 


KlX’.K  I)1SPIAc;KM?;NT  of  a SKMI- ink  INITK  PIJVTK  IINDKR  NORMAL  STiu;ss 

Lot  the  x-coord iiiato  lie  along  the  edge  of  a semi  - inf  ini te  plate 
witl  the  v-coordinate  lying  in  the  plane  of  the  plate  and  pointing 
into  it.  Snppose,  also,  tliat  a normal  pressiire  p,  defined  as  a general 
harmonic  liinction,  acts  on  this  edge  and  that 
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wlionct',  lor  tin-  non-trivial  sohitioii, 

. 2,2  . 0 

7sy 

Tho  solution  ior  I'(y)  is  tlu'n 

f -/v  -^y  ’■'‘v  '^y\ 

F(y)  = ' I’o  Qyo  •^+  Mo  Nyo  ^ ^ 

so  tliaf,  now, 

^ -/v  -vfv  '■'v  , 

■r  = i Pe  •+  gyo  Mo  Nye  (C  sin  i^x  + D cos  ■'k) 

Kvulontly,  in  order  that  houndary  condition  (b-3)  may  ho 

M = N = 0 

Thus  , 

'j  = o"'^^(P  + Qy)(C  sin  yx  + I)  cos  i'"x) 


xy  ^y 

= /(c  cos  yx  - I)  sin  ‘'•x)!* -Kc  '^^(P  + Qy)  + Qc  i 

i -» 

In  order  to  satisfy  boundary  condition  (li-4)  , 

''■(C  cos  yx  - 1)  sin  '■'x)  (-yp  + Q)  = 0 
either  C cos  yx  - D sin  '‘‘x  = 0 
or  Q - yp  = 0 

If  the  first  condition  is  true,  boundary  condition  (B-2) 
satisfied  only  in  the  special  case  of 
A sin  yx  = -B  cos  Kx 

and  not  in  all  cases.  This  suggests  that 


sat  is  f led , 


will  be 


C cos  yx  - I)  sin  yx  ■/  0 
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blit  that 


or, 


g . KP  = 0 


Q = i^P 


Tin'll 


+ vy)(P(:  sin  /x  + PI)  cos  /x) 


or , 


^y=  “t 

^ ?1X 


’'^^(1  + '^y)v'^(PC  sin  + PD  cos  i^x) 


2 


rf  1 = - / (PC  sin  /X  + PD  cos  v^x) 

yiy=0 


Boundary  condition  (B-2)  then  requires 

-(PC)K^  = A 
2 

.(PD)K-  = B 


Then  becomes; 


f,,  = _ + Ky)  ^ (A  sin  Kx  + B cos  vx) 


It  follows  that 


- i^y)  (A  sin  vex  + B cos  Kx) 

^ ^y 
2 

= e"'^^(l  + i'yXA  sin  Kx  + B cos  i^x) 

The  deflection  at  the  edge  of  the  plate  is  then  given  by 

^ ) 

;^y  L ''  y x/ 
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AFPI-.NIJIX  C 


PHYSICAI.  PROPKRTIKS  i\Nl)  PARAMin'KRS 


Tho  fdl lowing  dimensions  and  material  properties  have  been  used 
■vervwhere,  except  in  Chapter  11. 


Mi’tal  taco:  aluminum, 


5 in.  diameter 


.(Jbb  in.  width 


Oi 1 : SAb  U)  at  150  F, 


= 28.74J  in- Ihs/sec-in-  F 

k.  =0. 133  in^/sec 
M 

K = 1 X 10^  Ib/in^ 

a = 1.24  X 1()‘^  “f'*^ 

M 

= 1.8  X 10  ^ Ib-sec/in^ 

k = 1.01  X 10  ^ in^/sec 
oi  1 

K = 0.016  in- Ibs/sec-in-'^F 


Parameters : 


W = 1.613  X 10 


f = 6.25  lbs 


<=0.67  in  for  one  complete  wave 
< = 1.33  in  ^ for  two  complete  waves 


In  Chapter  11,  the  diameter  is  2 in.  and 


< = 1 in 


n ~i  o o o o 


Al’l'KNIHX  » 


I'Kik;RAM  KOK  DKTKKMINING  critical  Ol’KRATING  CONDITIONS 
FOR  THK  INITIALLY  WAVY  StIKFAGF 


f ms  pooL'iMM  Allows  ful  unimiom  chaijoks  in  hhaw 

C hhaT  is  TMf  .^OnDIMENSIONaL  SmaLL  ChAMOE  IN  HHAH 

C hhaTuIM  is  TmF  dimensional  small  change  in  hhah 

r 
c 

<IM‘'NSIOO  MSI),  MO  f S)  ) , r (SI  ) ,PP  (SI  ) ,P1  (^1  ) ,HT  (SI  ) 

OIMfNSlON  H(Sl;,  A(SI,S1),  H(Sl,n,  XbiSI,!),  G(S1,S1),  IK(Sl) 
dimension  PF (Si ) ,PL (SI ) (51 ) ,BH1  (SI ) ,Al  (SI ) ,BB (51 , 1 ) ,aX (SI , 1 ) 

dimension  YY(SI),R0(S1M>»XL(S1,1) 
hhaT 1=0.0 
hH A ^ =1.0f.  - OS 
r'x  = ?.o^3.  h,isdtssa/so.o 
D(j  IJ  1 = 1,  SI 
DO  13  J=1 ,S1 
13  A(I,J)=0.0 

00  IS  I=^,so 
A ( I , I ) =-.'>.0 
A(I,IM)  = 1.0 

A ( I , I - I ) - 1 . 0 
IS  CONTINOI 

A ( 1 , 1 ) = I . 0 
A (SI ,S1 ) =-l  .0 
A (SI  , 1 ) = 1 .0 
DO  01  1=1 ,S1 

Ai  ( n = *Da‘‘  ( i-n 

x(l)=  + Dx'*(I-l) 

01  continue 

J is  The  SLiOlMD  bHl  Eu*  NUNO  I mLNS  I ONAl  I I T H RESPECT  To 
The  ThEOkUICAl  CRITKAL  SpFIO  FOR  INITIALLY  FLAT  SURFACES 

0=1.0 

IN()=  1 

Y0=H«AR'->r  .01 
YPR: ( YO/ 100.0) /HBAR 
SSSS  CO'jTlNUr 

H0(I)  IS  THk  non-dimensional  surface  profile.  H(I)  The  FILM- 
ThICkNESS  and  <(I)  a small  perturbation  on  HU(I) 

DO  0?  I = I .SI 

HO  ( I ) = -Y  )*SIN  ( A ( I ) ) 

H ( I ) = (HO ( I ) .Hb^R) /HHAP 

0^  CONTINUE 

DO  04  1=1. SI 

Y ( I ) r-YPRosiN(A ( I ) ) 

HT ( I ) =Y  ( I ) 

OY,  COMTINlIf 
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n n r,  n 
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MlM  = 1 

cofUiNur 
PWIST  IIUNUM’P 

1 1 1 r Or  MAT  ( 1 H(i , I n A < oMuM=o » I ? t «o=* » F 5.? ) 

x N(i  '-1=  I 

nn  cofjriNur 

c 

r :5l(I)  RtFTKj  To  THF  2tK0  Qf^OF  W PRE  SbURF.  Irj  THt  FIlM.  WHEpEVEp 

C 

r 1>1  I<-  fjLoATlVEt  IHE  FIlm  Is  ((JMSIOEPEU  BWOk  F r;  ArjU  Pi  SET  TO  /EPQ. 

C P.^(I)  IS  THE  .•'lufJ-OIMtNsIFirNlAL  PpFSSUpF  EEk  T UPPA  T I OM  AMO  IS  SE  T lO 

C /Ff^O  ‘H[  FltB  IS  hRUKFM. 

C 

ol<n=-3.0»U'->(rt(r')-H(Sn))/(UXOM(l)*o3.o) 

IF  (Pi  (1)  .ET.0.>.)S01  »5r)? 

SO?  niMTlMUf 

SI  ( I » = ( OP-  i.r^UO  ( Y ( 1 ) /H ( 1 ) ) > « <H ( 0 ) -H (SO ) ) 

so  M sc 

Pr  ( I ) =U« ( Y ( ? ) - I (bO  ) ) 

{ 1 ) = - 1.  n«  (PH  1 ) ‘PF  ( 1 ) ) / (fMOH  ( 1 ) *03,0  ) 
pp  ( I ) = ( 1 .0/K{  1 ) *<»A.O  ) * (H  (?)  -H  (SO  ) )/  (?.0*UX) 
bOl  PI ( 1 ) ^0.0 
P?( 1 ) =0.0 
PP ( 1 ) =0.0 
50  Cnr'TlN'JF 

PI (51)=P1  (1  ) 

P? (51 ) =P?  ( 1 ) 
pp ( 51 ) =PP  { 1 ) 
on  30  i=?»so 

(M  ( I ) = -3.0«U«  (>'(  I ♦ 1 ) -H(  I - 1 ) ) / (nx<^H  ( I ) *'i3.0> 

IF(P1 (I) .LT.0.")503»50^ 

50A  continue 

PF ( I ) =U* ( Y ( I ♦ 1 ) -Y ( I-l ) ) 

HI  ( I ) = (OP-  ( Y ( I ) /H  ( I ) ) ) '-»  (M  ( I ♦ 1 ) -H  ( I-l  ) ) 

P?  ( I ) = -3  .(l«  (pF  t T ) ♦PI  ( I ) ) / (UX«U  ( I ) <»*3.0  ) 

PP  ( I ) = ( 1 .0/H(  I)  »«A.O)  o (H  ( I ♦ 1 ) -H  ( I-l  ) ) / (?.0*l)X  ) 

(1O  TO  3? 

503  CONTINUE 
Pi  ( I ) =0.0 
P? ( I ) =0.0 
PP ( I ) =0.0 
3?  CONTINUE 
30  CONTINUE 
4?n  CONTINUE 
PAV1=0.0 
00  101  I=?,50 
PAvl =PAv I ‘P? ( 1 > 

101  CONTlNUF 

PAY  IS  THE  ufilFoPM  (AvFrAC.E)  COMPONENT  OF  P?  THAT  CAUSES  THE 
VERTICAL  LIFT  "HAT  wITH  A FLOATING  MOUNTING. 


15  J 


|JAV  = ( 0 .s*  (Pr’  ( 1 > ♦iv  ('51  ) ) ♦►"AV  1 ) /SO  . 0 

PPI'JT  9lr.«F-'AV 

91R  FOPHaT  ( lO^.^A'^L^AGi:  PPt  SSUPL  = " ♦ F.  1 . S ) 

1 F ( \MliM  ,(.T  . ! ) 1 9 , AOO 

^00  /OtjTIrAH 

'-‘Avi=c.S’^  (pr(  1 > ♦i’fm'31  ) ) 

HAV?=0.0 
OG  401  I^r’tSO 
HAV?  = riAV^*PP  { I } 

^01  coriTinut 

hAv=hAV 1 ♦HAv? 

PPI9T  AO^.HAV 

*♦07  FOP9AT  ( I'lX  ♦'•*HA»=*,F15.b) 

'■iHAr  = - (bo  .0"PAv  ) / (IH.n^U^HAV) 

1)0  414  1 = 1. bl 

*♦1^,  i.)p(j)=pp(j)oib  . 0*U*HHAT 
M 17  COllIlNUF 

PP ( I ) -P?  ( I ) .PP  I I ) 

4 14  CONTI  f'^UF 

»<Nij'^  = KNuM*  1 

GO  TO  ^30 
419  COllTUiuf 
C 

C HEPF  P ’ PFOuMfS  THl  /lP.)  AVilPAfjF.  ULFoPmATION  CAUSING.  CumPONFNT 

C OF  P?. 

C 

GO  40P  1=1. bl 
P? ( I ) =PP ( I ) -pAv 
r ( I ) = r ( I ) *HHAT 
*♦0  7 CONTI  NUF 

HhATOlM  = HHAT'»UPAP 

PPINT  40J,  (HHAf  ♦HHATl)IM,MrtAR) 

4 01  FORMAT  ( 1 Ok  .*HMMT  = ‘»»Elb,b»bx»*HPATOlM  = *.Elb.b»bX»*HBAR  = *,Flb.b) 
hma  T 1 =HhAT  1 ♦HH/,TD  I M 
lOl  lONTlNGt; 

40b  CONTINUL 

Pl=0.b*(Y(l)/(,i  (!)*«?. 0)»t(bl)/(H(bl)««2.0)) 

H?=0.0 

no  OG  1=7, bo 

= ( y ( I ) /H  ( 1 ) »«7.0) 

OG  CONTiNUf 

00  10  I = 1 , b 1 

M ( I ) = ( (U»U*nX)  / (?.0»3. 14lS9?Gb4  ) * ( B 1 *0^^ ) -U*U*  Y ( I > /H  < 1 ) * *2 . 0 ) 
ri  ( 1 ) =H  ( I ) *OX«'»t:  .0 

10  continue 
H ( 1 ) =0 .0 
H (bi ) =0.0 
C 

C xsd)  are  TrIE  TUERNOElASTIC  DISPLACEMENTS  CAUSED  BY  THE  SmAll 

C PFRTUf^BAT  ION  Td).  Bd)  ARE  THE  RIGHT  HAND  SIDES  OF  EQ.(4-37)  OF 

C The  text.  HEPi  . heating  IS  CALCULATED  ON  THE  BASIS  OF  A FULL, 

c coherei'JT  fiim  ihroughout.  xsamp  and  yamp  are  amplitudes. 

C 


154 


call  MVS  ( A.Sl  ♦‘-.1  ,B»  1 »XS*IjFT  » IUP  T «G»  IRt  If  ^') 
son  format  ( iox,»if:K  = »,  i^) 

IF  ( I Ft’. I IJ.  1 )(,o  TO 

<SMAX=AMAX 1 (xS( I ) ♦aS<^)  tXS(3)  »XS(4)  *XS(S)  tXS (0»  »XS(7)  .X5(H)  » 

1  xs(^)»XS(10),Xs(ll),xS(lf-')»XS(M),xS(l‘*)»XS(lS),XS(lF>),XS(W)» 

r xS(lH),XS(lQ)♦AS(?0)^XS(^l),XS(2^),XS(^'3)^xb(^"4),XS(?b),Xb(^IS), 
3 xc,(?7)*xr,  (2Q)*AS(29),XS(30),xS(31)*XS(Jx’)»XS(3T)»XS(J4),Xb(J">)t 
^ *S  ( Th  ) . xs  ( <7  ) » AS  ( 3H  ) , XS  ( J9  ) , XS  (40  ) » XS  (4  1 ) , XS  (4?)  ^ XS  (43)  . XS  (4>4  ) . 
s X S ( 4-5  ) ♦ xs  (4X,  ) , AS  (4  7 ) , XS  (4M  ) , XS  (49)  , XS  (SO  ) , XS  (M  ) ) 

'SMIM=AMirjl  (XS(1  ) tXSI/)  »A:-,(3)  »XS(4)  »XS(5)  »XS(f))  »XS(7)  ♦XS*'".)  ' 

1 xS(9).XS(10).A^(n),X5(l<")«XS(13)*XS(l4)»xS<l‘j)«XS(lh>.XS(W), 

2 xs(  IS)  «XS(  19)  *as(20)  ,XS(21  ) ♦xS(?2)  tXSlx’J)  tXS*?''*)  *XS(2S)  ,Xb(2h)  * 

3 xs(?7)»XS(2B)*aS(?9),XS(30),XS(31)iXS(32),XS(33)»XS(34),XS(3S), 

4 *s  ( 30)  , X',  ( 3 7 ) » AS  ( 3B  ) ,XS  ( J^)  ♦ XS  (40  ) , XS  (4  1 ) , AS  (42  ) , XS  (4  ))  , Xs  (4-4  ) , 

5 xs(4S)tXS(4(O*AS(47),XS(4M),xs(44),xS(‘j0),XS(‘3l)) 

XSAMP  = AHS  ( X5MAa-xsM  irj) 

yMAx  = AMAA  l(Y(l)»Y(2).T(3),r(4),Y(S),Y(t>)»Y(7),r(H),Y(9),Y(10), 

1 Y(ll)»Y(l?),Y(13),Y(14),t(lS),Y(10)»Y(l7),y(18)»Y(19),Y(20)» 

2 Y (21  ) . Y (22)  , Y (X.-3)  , Y (24)  , f (?S)  , Y (26)  , Y (?7)  , Y (28)  ,r  (?9)  , Y (30)  * 

3 Y ( 31  ). Y ( 32)  , Y ( j3)  ,Y  ( 34 ), T ( 3S) , Y ( 36)  ,Y (37) , Y ( 38) » Y ( 39)  , Y (40  ) « 

4 Y(4l).Y(42)*Y(H3)»7(44),r(4S),Y(46)»Y(47),Y(4H),Y(49),Y(50), 

C Y ( S 1 ) ) 

YMIM=AMINl(Y(l)»Y(2).Y(3)tY(4),Y(5)»Y(6)fY(7)»Y(8),Y(9),Y(10). 

1 Y(n)vY(12),Y(i3),Y(14),Y(lS)»Y(16),Y(17),Y(lH),Y(19),Y(20)» 

2 Y(2n.Y(22),Y(23),Y(24),Y(25),Y(26),Y(27),Y(28)»Y(29),Y(30). 

3 Y(3n*Y(  <2).y(  l3).f(34),Y(3S),Y(36)»Y(37),Y(38)«Y(39),Y(40)» 

4 y(41).Y(‘+2),Y(4  3)»Y(44),y(4S),Y(46),Y(47),y(48)»Y(49),Y(B0)« 

5 Y (SI ) ) 

YAMD=ARS( YMAX-rMIN) 

PRINT  1 12,xSAMo,yAMP 

112  format ( ihO. 1 OX ."XSAMPr®  tF lS.6»SXY*YAMP  = *Yt 15.6) 

1 F ( XSAMR  .1,  T . Y AMP  ) (-,0  TO  3333 
f.n  TO  9990 
33  33  corjTiour 
OOM=NUM* 1 
no  03  1=1,51 
Y ( I ) =Y ( I ) -HHAT 
0 3 '‘OflTlNOF 
11=1 ( ♦ n , 2 

IP (O.OT . 1 .8 ) 999V.6666 
HSMH  CONTINUf 
9999  COrjTlNiJF 
0=1.0 
I NO = I NO ♦ 1 
fO=Yn*KHAl *0.0 3 
IF  ( I.NU.Gl  . J)22p2,S5S5 
222?  CONTINUf 
STOP 
fno 


Al’PKNDIX  K 


riU.  KU-XIIU.Y  MOINTKI),  INITIALLY  WAVY,  SLAL  LACK 


In  practice,  one  of  the  faces  is  usually  flexibly  mounted  to 
accommodate  axial  misalignment.  If  both  faces  are  nominally  flat  but 
axial  misalignment  prevails,  the  fluid  pressure  creates  a couple  Lliat 
tilts  the  flexibly  mounted  face  and  corrects  the  misalignment  (fig.  16;. 
If,  however,  the  faces  are  axially  aligned,  tilting  will  still  occur 
if  one  face  is  wavy  and  the  wavelength  Is  equal  to  the  circumference  of 
tlie  seal.  If  there  are  two  identical  waves  around  the  circumference, 
there  will  be  no  tilt  because  the  circumferential  pressure  distribution 
will  be  symmetric.  If,  however,  the  face  deforms  during  operation  so 
as  to  destroy  this  symmetry,  tilting  moments  will  arise  to  counteract 
the  asymmetry  in  the  pressure  distribution. 

The  following  analysis  appertains  to  the  response  of  a flexibly 
mounted  system  to  small  perturbations.  It  is  assumed  that  a steady, 
non-tilting,  configuration  prevails  and  is  then  disturbed  slightly. 

Ihe  subsequent  pressure  fluctuation  and  tilting  moment  are,  therefore, 
of  the  order  of  the  perturbation.  The  tilting  moment  may  be  resolved 
into  its  components  about  the  A and  15  axes  (fig.  L-1)  as  follows; 

dT  = p'(L-dx)r  sin  b 
A * t 

(K-1) 

dt^  = p^(L'dx)r  cos  b 

where  p^  is  the  appropriate  pressure  perturbation.  In  terms  of  non- 
dinR'nsional  variables,  the  ccjmponents  of  the  net  tilting  moment  become 
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” sin  X dX 
L 


(K-2) 


2n 

P'cns  X dX 
I 


(K-'3) 


Thcsi'  components  tilt  the  face  so  that  the  film-thickness  changes 
'.hape.  Such  changes  are  also  of  the  order  of  the  perturbation  and  arc- 
denoted  by 


H,  = X 

A A 


(K-4) 


0:-5) 


The  faces  can  also  have  relative  axial  motion  (fi)  so  that  the  non- 
dimensional  film  thickness  is 


lUX) 


H(X)  -H  H + Y(X)  -b  H^(X)  + lIj^(X) 


(K-6) 


The  non-dimensional  pressure  perturbation  is 


T - 1-  - f - ^ 'h  <"a  " - 


dHN 

H dX/ 


(F.-7) 


where  1’"  and  P are  as  defined  by  eqs.  (4-24)  and  (4-26). 

The  tilting  moment  is  neutralized  by  the  tilt  defined  by  and 
Thus,  using  eqs.  (K-2)  and  (K-3),  one  has  the  conditions 

2n 


and 


= P'sin  X dX  = 0 
A t 

o 


2n 

/% 

. •=  P'cos  X dX  = 0 

B . t 
o 


(E-8) 


(E-9) 


1 
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It  lias  luHMi  i-xplaiiu'd  in  soction  (4.5)  tliaL  tlio  pressure  per tiirbaLion 
p'  must  havi'  a zi'ro  average.  1 ii  the  present  case,  hv  analu>;y,  I’^  must 
have  a zero  average.  Thus  a tiiird  condition  to  he  satisfied  is 

2’’ 

l’[  dX  = 0 (K-IO) 

o 

Tquations  (K-8),  (K-9)  and  (K-10)  contain  tlic  three  unknowns  II,  and 
ihey  cannot  he  solved  simultaneously  hy  analytical  means.  To 
solve  them  numerical Iv,  they  are  first  reduced,  hy  a straightforward 
but  cumbc'rsome  process,  to 


wilt  re 


II.  = 


A 6u5 


t 'X 

" ■ *:  "a  - dr 

H 4 


H..  = T-  (iji.  - t.,H  - 


B ■ V 6 A 


_9 

6uJ 


1 ~ *9*20  *10  ' *5*19 


i i 

5 18 


) =4$  - i +44 

2 6 20  lb  3 19 


4 4 
3 18 


4.. 


* ^4 

14  4^  *15 


I - 4 

12  4^  15 


11  4^  15 


(K-ll) 


(E-12) 

(E-13) 
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’ I ' V • V ^ ‘I"  • , -v 

— icos  X sin  X - 77-j^  sin  X / dX 

» , , J 11  cl  A ' 

11  It 


' 1 ^ ^ 3 dll  , V 

— 2 , sin  X + — ~ sin  X cos  X 1 dX 

o H 


, ' 1 clli  . „ 

o H 


P”sin  X dX 


10 


P"  dX 


11 


P"cos  X dX 


12 


'*  1 dll  „ 

3 —7  *77  cos  X dX 

I 4 dX 
o H 


13 


2’T 

J_  ^ HX 
o 11 


2n 

''  1 ^ 2^  3 dH  , V „ N 

“7  ' cos  X - -7  -77  sin  X cos  X 1 dX 
14  ,,3  H dX 

o H 

TT" 

" 1 . V V ^ d»  2,  ' 

f,_=  — :^^sinXcosX  + 77—  cosX.dX 

Ij  ^ II  uX  / 

o H 


16  ^ 3 1 

o 11 


i 1 y 
11  ^ dX/ 


■ cos  X - f:  sin  X dX 
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[sin  -'^  + 11  f''!- 


V dll  ' ,y 

c i>s  X — ; dX 
dX/ 


8 . 
'4 


20  1. 


KquiUions  (l.-lli,  (K-12)  and  (K-13)  are  solved  numerically  Cor  11^,  II 


and  H . 
H 


lo  analyze  the  flexibly  mounted  seal  for  thermoelastic  instability, 


eq.  (4-21  I is  solved  as  described  in  chapter  IV,  where  Y includes  H 

A 

and  H, . 
n 

In  describing  the  experiments  in  Chapter  111,  it  has  been  observed 
that,  as  a con  .equence  of  the  onset  of  instability,  the  two  faces  come 
into  direct  contact  (fig.  26).  To  test  whether  this  can  be  obtained  in 
theory,  an  experimental  profile  has  been  measured  at  an  operating  con- 
dition just  preceding  the  onset  of  instability.  The  non-dimensional 
critical  speed  has  been  taken  as  unity.  A small  sinusoidal  perturbation 
is  imposed  on  the  mi'tal  surface  and  eq.  (4-21)  is  solved  for  It 

is  lound  that  |y1.  The  following  definitions  are  then  introduced; 

Y (X)  = o-Ia  , Isin  X 
new  ' th' 


where 


O'  .05 


H (X)  = H (X)  + A , (X)  - Y (X) 
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llu'  solution  is  ropeated,  with  tlio  new  overall  film  thickness  H and 

new 

the  nt‘w  perturbation  Y . This  cycle  is  repeated  over  and  over  aeain 

new  ' ^ 

and  it  is  found  that  the  surfaces  do,  indeed,  touch.  Such  a solution 
has  been  obtained  for  the  case  of  the  floating  hut  non-tilting  mounting 
as  well  as  the  floating  and  tilting  mounting.  The  results  are  shown  in 
fig.  ( 1.  - 2 ) . 

It  has  to  be  pointed  out  tliat  this  procedure  does  not  take  full 
account  of  the  dynamics  of  operation  and  therefore  must  not  be  regarded 
as  a simulation  of  the  long-term  growth  of  a disturbance.  It  does, 
however,  simulate  the  short-term  growth  of  a surface  perturbation  and 
proves  that  if  one  surface  comes  very  close  to  the  other,  it  can  then 
also  touclj. 


o o o r>  ^ o 
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PROr>f?AM  aTLAST  ( INPUT  »OUTPUT  ) 

C 

r 

r This  program  ai  lovis  foi^  a nf  t lI^t  (hhad  as  <i»Fll  as  a tilT  (Composed 

r or  ( hAhai ) siNx* (hhhat ) COsA ) , 

^Id)  REFEF^S  to  The  total  -'CWO  orued  pressure,  that  determines  the 

REGIONS  OF  CA'^  • TaT  lOrj.  IN  THF  CAVITaTING  REuIONS*  NONE  OF  THE  EQOaTIUNS 

apPuT.  p?(i)  is  the  component  of  thf  first  order  perturbation  arising 
FROM  Th(  SUREAi  E PeRTUKraT  lOfJ  Y(I)  alone*  and  RS(I)  is  The  component  from 
MlT  and  net  LiFT.  wheRLVeR  PI  IS  ZfRO*  P2  Anu  PS  ARf  SET  TO  ZERO. 
FIt!ALLT»  P?  ANR  F’S  ARE  AUDED  To  PRODUCE  A NEiV  R^»  WHICH  Is  THE  TOTAL 
perturbation  R. ASSURE  ANu  HAS  A ZERO  AVERAGE  AND  PRODUCES  ZERO  TILT.  THE 
C (O^OITION  or  ZlRO-AvERAGE  pz,  together  with  The  zero  tilt  condition, 

c PROVIDE  Three  louations  for  hhat.  hahai*  hbhat, 

c Thf  new  PP  is  used  OIREcTLV  for  determining  elastic  deformation. 

c 

c 

DIMENSinn  X(Sl),  ho (SI ) »r (SI ) »PP1 (51 ) ,RPZ(S1 ) ,PRJ(51 ) ,PI (SI ) 
dimension  h(SI).  A(S1,5I)»  H(5l»l),  )(S(51,1),  G<S1»51)»  IR(51) 

UImEnSION  PF (Si ) ,Pl (51 ) »PP (51 ) ,BB1 (51 ) ,XI (51 ) ,HB(51 » 1 ) ,XX (51 . 1 ) 
u I KENS  I On  YY (SI ) ,H7 (51 ) , A3(51 ) ,H3(51 ) ,FZ(51 ) ,GZ (51 ) 
dimension  P 3 (PI ) »RA (51) ,P5(51 ) ,M  (51) ,A? (51 ) ,B1 (51 ) ,B? (51 ) 
fd  MENS  I ON  CZ (Si ) »DZ (SI ) ♦HA WAVE (51 ) ♦HhwAVE (51 ) ♦HMWAVE (51 ) *Cl (51 ) 

DIMENSION  TA(Sl)fTH(51) 

UP=0.0 

Hf3AR=l  .A8609E-VS 
Ox=P.O*3. lAlSR. GSA/SO.O 
HS=3. 1Z79PE-0A 
DO  13  1-1,51 
00  13  J=l,51 
13  A(1,J)=0.0 

DO  IS  I=Z,50 
A( 1 , I ) =-P.O 
A(I.IM)=1.0 
A( I , 1-1 ) =1 .0 

16  continue 

A ( 1 , 1 ) =1 .0 
A(51.51)=-1.0 
A (SI , 1 ) =1 .0 
00  0 1 1=1,51 
X 1 ( I ) =*DX*  ( i-n 
X ( I ) =*UX« ( I-l ) 

P5( I) =0.0 
01  CONTINUE 
5555  CONTINUE 
U = -l  .0 
NUM=  1 

YO=HBAR/10.0 


n o 


1 f)''* 


^ef'D  IN  l nTAi  I I Mf'ASInn  i)  DlMl.NSlUNAL  SUWFACr  PWOFIlt 

r p(.'t  croirjG  thf  unstahlf  mati 

c 

«LA0  ‘♦OP,  (H0(  n , 1 = 1 ,51  ) 

‘♦OP  FOR'^AT 

DO  0?  1=1 ,51 

Hf)  ( 1 ) =HO  ( I )♦(-,.  1 lr’5F-0t) 

M ( I ) = (HO ( I ) ♦HO, a-;)  /hhar 
00  COrjTiNUF 

rpR=H ( 1 J) /?5.0 

1)0  OA  1 = 1,51 

V ( 1)  =-yF^R*SIM  ( A ( I ) ) 

OA  CONTINUI 

ppi Nil  u 1,  (M ( n , I = 1 ,si ) 

U 3 format  (/,51  (/,lnx,tl5.5) ) 

CONTINUF 

PRl-gl  111,  MUM 

111  format  ( iHl  , 10A,*NI)M  = o,  I?) 
knjum=  1 

MNjijMr  1 

7 7 77  COfjTlNUF 

PI  (1  )=-3.0“0<»<h(P)-H(50)  )/(DX*H(l  )**J.O) 

IF  (PI  ( 1 ) .LT,0.>->501  ,50^ 

50?  CONTINUF 

PF  ( 1 ) =U» ( Y ( ? ) - I ( 50 ) ) 

PL  (1  ) = (OP- J.L’"u»  ( Y ( 1 ) / h(  1 ) ) ) <»  (H(?)-H(50)  ) 
p;?(l)=-3.0<*(pLil)*pF{l))/(i)XoH(l)o«3.0) 

PP3 ( 1 ) = ( 1 .0/H ( i ) o (H (?) -H (50) ) / (?.0*UX) 

A3 ( 1 ) =COS ( X ( 1 ) t 
h3( 1 ) =SIN(X ( 1 ) ) 

F^(1)  = (A3(1)-(j.O/h(1))*h3(1)»((h(/?)-h(50))  / (?.0*()X)  ) ) /H  ( 1 ) **3.0 
0?(l)=(H3(l)»(J. 0/m(1))oA1(1)*((H(?)-h(5U ))/(?. 0*DX)))/H(1)**3.0 
50  TO  50 
501  COrjTlNUf 
PI ( 1 ) =0.0 
P2( 1 ) =0.0 
PP3 ( 1 ) =0.0 
F? ( I ) =0.0 
( 1 ) =0.0 
50  CONTINUF 

PI  (51 )=P1  (1  ) 

P?(51 )=P^(  1 ) 

A3 ( d1 ) =A  3 ( 1 ) 

H3( 51 ) =H3  ( 1 ) 

PP3 (51 ) =rP3 ( 1 ) 

F?(5I )=F?(  1 ) 

( 51 ) =5P  ( 1 ) 

DO  30  I=^,50 

Pl(I)=-3.0*l)«(fi(I*l)-H(I-l))/(f)X*H(I)*o3.0) 

If  (Pi  ( I ) .1  T .G,-)5n3,50A 
50A  COMTINOF 


A 
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PF ( I ) =U* ( Y ( I ♦ 1 > -t ( I-l ) ) 

PL  ( I ) = (UP-3,3*U*  (V"  ( n/H(  I ) ) ) * ) -H(  I-l  ) » 

P?< I )=-3.0* (PF ( I ) ♦PL ( I ) ) / (UX*M( I ) **3.0) 
a3  ( I ) =cor.  (X  (1)1 

P3( I ) =SlN  (X  ( I ) j 

PP3(I)  = (1.0/H(i)»«»44.())o(h(I^l)_H(I-l))/(i:',0*L)X) 

F?(I)=(A3(I)-3.0/H(I)*H3(I)»(H(I^l)-H(I-l)) / (2.0*UX) ) /M ( I ) •* J.O 
0^( I ) = (H3 ( I ) ♦ J.O/H( I )*A3(1)*(H(I^1)-H(I-1) )/ (a.0*UX) )/H( I ) **3.0 

GO  TO  3? 

503  COnTlNUi 
Pi ( 1 ) =0.0 
P^( I ) =0.0 
A 3 ( I ) = 0 . 0 
P3 ( 1) =0.0 

PP3 ( I ) =o.n 

F?{ I )=0.0 
G? ( I ) =0.0 
3P  COIJTIGUE 
30  CONTINUF 
C 

C BEGIN  calculating  TIlTING  MOMEnT*  HHAT ♦ AnO  TIlT 

C 

Do  701  1=1.51 

PPl  ( I ) =t  f’3  ( I ) *o3(  1 ) 

%.  pp?(I)=PP3(1)*A3(1) 

P3 ( I ) =P2 ( I ) *HJ ' I ) 

P4 ( I ) =Pa ( I ) *A3  1 1 ) 

A2 ( I ) =FP ( I ) *53  ( I ) 

H?  ( I ) =G2  ( I ) «->3  I I ) 

( 1 ) =^^ ( 1 ) *A3» I ) 

D? ( I ) =GP ( 1 ) *A3  U ) 

701  CONTINUF 

hAV1=0.5* (PPl ( 1 ) ♦PPl (51 > ) 

HA\/?  = 0 . 0 

HA V 3=0.50 (PP^( 1 ) ♦pp?(Hi ) ) 

HAV4  = U . 0 

hAv7=0.5» (PP3( i ) ♦PP3(bl ) ) 
hA\/5  = 0 . 0 

P31NTl=0.5*(PJtl)^P3(51)) 

P31ST^=0.0 

P^r-IT  1=0.5*  (PA  (1  ) .06(51  ) ) 

P'.riT^  = 0 .0 

HlrjTAl  =0.5*  (At  i 1 ) .APlbl  ) ) 

hinta<?=o.o 

HlNTHl=0.5«(Ra(l ) ♦H?(51 ) ) 

HINT5i^  = 0.0 

Hiwrci=0.5* (C^( 1 ) .CPlSi ) ) 

HiNTry=o  .0 

hinTD1=0.g* (O^t 1 ) .uPlbl ) ) 

HINTOt:  = 0.0 

HINTF1=0.5*(F^^ 1 ) .F^lbl  ) ) 
mINTF^=0.0 

hINTG1=0.5»{GP ( 1 ) ♦G?(51 ) ) 

MlNTr-t?  = 0.0 


1 •)() 


00  '>01  I =<-'<.50 
hAV?  = hAV/‘>PP1  ( 1 ) 
mA\/>  = '-*A\/‘4  ‘PfV  ( i ) 
hAvH  = i-iA\/H*PP  1 ( i ) 

‘^JIMr^  = PJINT^>»0(  I ) 

P4nT^  = P^INT^♦^  ■4  ( I ) 

HiKjr  a^’=hinta?*-j^’  < 1 ) 

HI^'TR<^  = HI^Tb^♦c^  < 1 ) 

M ] tJ  I f P = H J uTCc’  * ( I > 

•<It;T[V  = HlMTI)c'>uP  ( I ) 

►iINTfc>  = HlMTf  I ) 

HirjTG<’  = HINTr, ( I ) 

401  CONTINUL 

HAVE)=(HA^1♦HAV^)C3.0^^I>X 
hAV6=<HAv  1*HA  ) * 3.0«('X 
MA\/4=(hAv?>HAVo)«  l.noDX 
3 I '4  T = < P 3 I NT  1 ♦ t'  3 1 N T ? ) <»  OX 
P<4r4T  = (P4  IMT  1 INT?)  «0X 

HINT)  = (HINT  A1  ‘MINT  A?)  ’*0X 
HlfjT?=  (HU.  TP  1 ♦NINTH?)  <»OX 
H If 'T  3 = (HINTC  1 ♦MifvlTC?)  '‘OX 
hjmT4=(hipT01  ♦nIMTu?)  '•‘OX 
NlNT5=(HIMTTl*fiINTF?)'‘0X 
hINTH  = (HlNTC)l  ♦hINTO?)  '‘(JX 
4?0  COfJTlNUL 
PAV1=0.0 
on  101  I=?.S0 
PAV I =PAv 1 ♦P?  ( I ) 

101  CUNIINUI 

pAV=(0.So(p?(l ) ♦p?(01 ) ) ♦pAvl )/S0.0 
PSlM  = PA y/oPO  . 0'*uX 
PRINT  '^IdfPAv 

HIR  PORHaT ( I OX , »AVuPAGt  PPL SGORf = » , F j ^ . 5 ) 

I ( { HNOM.cJ  . 1 ) uu  TO  Oil 

400  continul 

SlGHl  =P4I‘.T-  (h.NTO*P3INT)  / hint? 

SlGH?  = HlNT3-(HiMTA'‘HINTl  )/HINT? 

SIGH3=HAv6- (HA^5*HINT4/HINT?) 

SIGH4= (HINTG^nrVS) / (hlNT^^SlNHT) 

SIOHS  = HAVP/SIGti3 
SIGHG^HUfTb/hliiT? 

SIGH7  = SU,H?»51>.HS-HlMTS*nINT  1 «‘S  I GH6-S  I GM4*S  I GH? 

SlOHH  = SlGH!  oSI'->HS-pS0M*PiIMT«SlGH6-Sl&H4*SlGHl 
hAHAT  = 5IGHH/  (6.0'‘U»‘SIGH7> 

hhaT=  (HAhaT*S1^H?)  /GIGh3-SI(iH1  / (b.0*U*SlGH3) 

MhHAT=( (HAHAT*"!''! I I-PJINT/CH.O^hI-IHHAT^HAVS) )/mINT? 

DO  oON  1=1.51 

HAwAvL ( I ) =HAHA I*SlM<x ( I ) ) 

HP»f  AVI.  ( I t =HRhA  1 *COS  ( X ( 1 ) ) 

HMWAVr ( I ) rHAWA^r  ( I ) ♦HHwAVr  ( I ) 


n n o 
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bOq  CONTINUF. 

HHAT  AMD  TILT  HAVE  NOW  BEEN  OBTAINED 
DO  511  1=1.51 

P5(  I ) =-6.0*lJ<»  (liAHAT«F  2 < I ) -HHHAT*62<  I ) -HHAT*3.0*PP3  ( I ) ) 

511  CONTINUE 

DO  512  1=1.51 
P2 ( I ) =P2 ( I ) .Pb ( I ) 

51?  CONTINUE 

DO  516  1=1.51 
TA  ( I ) =P?  ( I ) <>B3  ( I ) 

TH(I)=P2(I)*A3(I) 

516  CONTINUE 

TAIN)T  = 0.b*^  (TA  1 i ) ♦TA(5l  ) ) 
lBI''iT  = 0.S«(Tfa<l).TH<51)  ) 

Do  51^  1=2.50 

TAI^JT  = TAINT♦TA^  I ) 

TBINT=TUINT*THlI) 

516  continue 

T0R<A  = 0)(<»TA  INT 
rOR<B  = DX<>TBINT 
mnU'1  = MNUM.  1 
GO  TO  A20 
513  continue 

)0  '-.02  1 = 1.51  ■ 

>’(  I ) =M  ( j ) ohhAR-HBAR 
Y ( I ) =r  < n *HHAR 
HMwAvE  ( I ) =HMwAvE.  ( I ) oHBAR 
602  CONTINUE 

hS  = HS*HBAR»«‘?.i' 

HHATDIM=hhaT»HoAR 

riBAR=HBAR.HHATuiM 

hS=H5/HBAr«o2." 

Print  ^oj. (hhai.hahat.hbhat.hbar.torka.torkb) 

60  3 FOR'^AT  ( 10;(  .*HH6T  = *,E12.6.2x»«HAHAT  = *,E12.6.2x.*HHHAT  = *,E1  2.6.2X. 
1 «HBAR  = *.Cl2.6.^rX.*T0RKA  = "»E12.6.2X»*T0RKB  = «.El2.6) 
r.iO  6 06  1-1.51 
HMWAVE ( I ) =HMWA¥E ( 1 ) /hBAP 
H ( I ) = (H  ( I ) ♦H^1A(^)  / HttAR*hMwAvE  ( I ) 

HT (I ) =H ( I ) *HBAK-HHAR 
Y ( I ) =Y ( I ) /HBAR 
6 06  conriNUF. 
lOT  CONTINUI 

DO  201  1=1.51 
-BPrO.O 
(<0  ?02  J=I.5I 
16 ( J.FO. 1 ) 203.206 
•>H  1 ( j)  =-R?  ( J)  .0/  Ox 
, TO  20  > 


1 


n o 
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HMi  ( j)  ( j)  / (mhs  ( X { I ) - X 1 ( j» ) ) 

PO^,  COfJIlNuf. 

HHP=BHr»eH 1 { J) 

POP  CONTIUUL 

hB^  = BB^-u.'5^  (ttol  ( 1 ) ♦Bbl  (^3l  ) ) 

hh  ( n =HHr'*'HS'»LiA*n>3.0 

?01  COrjTlfNlUf 

btt ( 1 ) =0,0 

bB(51)=0.0 

computation  of  flastic  oft lfct Ion 

CAl^  MVS  ( a I .-^1  «HH,  1 «XX  ,I)F.T  » lOET  *G»  IP»  ILp) 

PPIviT  GOl,  IL'P 
601  FOfiMAT(lox,"lF^'  = ^,I?) 

IF  ( IFP.E'J.  1 )G0  TO  999S 
'*05  CONTINUE 

-'i  = o.s'*(Y(n/(ii(i)»oB.O)*Y(5n/(H{5i)*"?,0)) 

H^-0.0 

00  06  I=^,50 
Ra  = B?*  ( Y ( I ) /H  ( 1 ) ) 

06  continue 

00  10  1=1*51 

B(l)  = ((U‘>U'>0x)/(^,0oJ,lAlS9?()SA)o(Hl*H/^)-IJ<nj*Y(I)/H(I)oo/^.0) 
H ( I ) =E  ( I ) ooxo'^c  .0 

10  continue 

H ( 1 I =0.0 
H<si )=0.0 
c 

c COMPUIATION  Of  THEI-'mAl  OEfLECTION 

C 

CALl  MVS(A*51*‘->l*H»l*Xb»  UE  T » I OE  T » G » I R » I CW  ) 

PPIMT  500*IEP 
500  roRMAT(10x,*lEM=«,I?) 

IF ( IEP. EQ. 1 )GO  TO  9999 

00  ?06  1=1,51 
?06  YY ( I ) =XS  ( I ) ♦XA  ( 1 ) 

OPINT  y00,(P^‘i),XX(I),XS(l),YY(I),Y(I)»H(I),HT(I),I=l,51) 
HOO  rORMAT(/,51 (/,10x,7C15.5) > 
mh=  Y (13)  ’*HHAP 
hH0= Y Y ( I J ) *HBA„ 

1 ^ ( ABS  (OH  ) .r,T  ...MS  (MHO)  > 464^  t 1 J 13 
COtlTlNUf 

PRINT  91/ 

917  FOPMAT(/,iOX,*KERTUPHATIOr*  IS  BEING  SUPPRESSED*) 

(lO  To  P.PPP 


z>  n n n 
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3333  CONTINUC 

Hm1NJ  = AmI.N1  (H(l)»H(2),H(3)»H(4),H(5),H(b),H(7),H(b),H(9),H(10)» 

1 H(ll),H(l?),H(13),h(l^),H(15)»H(16),H(17),H(18)»H{19),H(20)» 

2 H(ai ) ,H(2?) ,h(c3) tH(29) ,H(25) »H(2b) ♦H{27) ,H{2H) ,H(29) ,H(30»  » 

3 H ( 31 ) »r (32) ,H (j3) ,H(3A) ,H( 35) »H (3b) ,H( 37) ♦HC 38) ,H (39) ,H<40) ♦ 

A H(41  ) »H(42)  ,H(^3)  ,H(45)  »H(46)  »H(47)  tH(A8)  ♦H(A9)  ,H(50)  » 

5 H < 5 1 ) ) 

PRINT  702»HMIN 

70?  format (10X»«hMiN=«»e15. 5) 

YPR=HMIn/25.0 
DO  03  I=l»5l 
Y(I)=-YPR*SIN(a(I) ) 

YY(I)=YY(I)-Y(i) 

H(I)=H(1)*YY(1) 

03  continue 

YTE5T=HMIN 

IF ( YTE5T.lt. O.M  )GO  TO  8888 

NUM=nUM* 1 

IF  (NuN»GT»100)'Y999»b6bb 

8888  continue 

PRINT  31b 

31b  FORMaT(/.iOX»<»THE  surfaces  have  TOUCHED") 

9999  continue 
2222  COfJTlNUE 

print  314. (H(I> fY(I)f I=I »5]) 

314  format ( IHl ,51 (/ » 10X»2E15.5) ) 

STOD 

ENO 

OOOOOOOOOOOOOOOOOOOOOy  END  OF  RECORD 

experimentally  measured  surface  profile 

♦1 .b350lE-0b*l .3834 /r-Ob*2. T6694E-0 7-2.7 lb63E-0b-2. 76694E-0b 
-4,9050 3E-06-7. 54b02E -Ob- 1 . 04389F-a5-l .289] 4E-Ob-l .37089E-05-1 .b3501E-05 
-1 .73059E-05-1 .83501 E-nS-l .42372E-05-1 .22877E-05-1 .22870E-05-1 .2?8bOE-05 
-1 .23000E-05-1 .24000E-05-1 .2b770E-05-l . 34b99E-0b- 1 . 38347E-05- 1 .5092hE-05 
-1 .530b2E-05-l .3029HE-05-1 . 1533 IE -05-1 .0061 bE-0b-8.b5298E-06-7.20662E-0b 
-7,20602E-06-5.96148E-Ob-4.03719E-06-?.9baoOE-06-l .03129E-06*! ,60494E-1 I 
-1  , 2577 OE-0 7-1 .38  34  7E-0b*8 . 80390F-0  7 ♦H . 30082E-0b*9 .b8429E-0b« 1 . 18224E-05 
♦ 1 .14451E-05*  1 .02503E-05*9. I 4 34bF-0b*8 .4b4  32E-06^ 7 . 70970E-0b*6 , 565 1 9E-0b 
♦b,34523E-0b»3.b2I5bE-O6»O.OOOOOE*OO 
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APPKNDIX  K 


PK(K',R.\M  KOK  COMPUTING  FAGl-:  hOAD  /\ND  INITIAL  KACK  PROKILK 
CORRKSPONDING  TO  GIVKN  OPKRATING  CONDITIONS 


C-CG.--5M  SOPO-^TI  INPUT  , r.GTPLT  I 
C 

c 

C C :‘'‘'’UTiT:CN  OF  SOPPCFTF"'  FPP  VflPICDI  amplitudes  op  ncl- 

r IoICnal  suppace  wave  sl*"  a -ang^  pp  slioirg  speeos.  loap 

0 I'l  GIVEN  TV  HE  NET  CPESSUPE  IN  THE  plUIP  FILM.  pQS  GIVEN  OUASI- 

r PTAIIC  P CLOT  I ON*^  HFPPATINP  CONOITIONp*,  THEIMTIAL  SURFACE 

C P = CFILES  A=E  ALSO  SAlC'JLATFP.  GINGLE-LCGED  SURFACE  PROFILES  APP 

0 CCMir.EPEO,  WHERE  GOTH  SUPFACEp  APE  pESTcaINEP  frcm  TILTING 

r 

c 

CH.N^ICN  T(5i)  ,H0{F1)  ,P1  (PI)  hIEI)  ,y  (51) 

"'EM  NS  I CN  A (5  1,  PI  > ,G  (El  , P 1)  , IR  ( 51  > , Hu  (5  1,  1 ) , XL  f El  , 1) 

-iPAr=  1.  ;;r 

C^-  .G'^.l^lSTEPEu/EO.'] 


A ( 1 , 1 > = 1 . : 

A (5 1, PI ) =- 1 . a 
A ( 5 : , 1 ) = 1 . c 


''  ' .1  I = 1 , E 1 
f (:)=»nx»(:-i) 

: 1 c p T I N p E 

U IS  TM-:  eliding  SPFEC,  RCNOIMENSIONAlISED  kith  respect  TC  THc 
Tp'ECPETICAl  SJITICAL  SPPEO  FCR  initially  flat  SLRFACES. 

u p : . 1 . 

T 1 p - 4 

Y''=i-PAu»;;,  ir 

H AMFr  VO /Mg  AP 
535E  C''NTtnuE 

H(l)  IS  Thp  .NONniMENSIONAL  FILM  THICKNESS 

GO  :E  1=1, El 

hh  I ) =-  vn*  s I N ( X ( I ) ) 

H(I)=(hC(!»yh'>ap)/hpar 
CO  CONTINUE 

Ni;w-i 

5 h F f CONTINUE 

PPINT  1 1 1 , N'U*',U  , HAMP 

111  F CR^'AT  ( IMS  , lOX,  ‘NUMs*  ,I  2,  EX  , *U=*,F5.  ? ,5X,*HAMP=  *,  FG  .3) 

7777  CONTINUE 


PlC*  •’■.FE^'  TO  THE  PPPPPUFE  IN  the  FILM.  WHPREVER  PKD  IS 
N GATIVE,  the  film  IS  CCNSIP(RE0  iTRCKEN  AnO  PKD  IS  SET  TO  7ERO. 
r PI  IS  NCNDIMENS lORAL . TO  OpTAIN  niMENSIONAL  PRESSURE,  IT  HAS  TO 

C PP  riviCE)  HV  Hgiju/w»5,r3PT  (ALPHA/(MU*K)  . 
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I 


► 


I 


f 


c 

c 


1C  1 


r 

C 

c 

c 

c 

r, 

C 


:’Ci 
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IP(Pl(l).LT.o.0)5ri,5C 

p 1 ( 1 ) = o . 1 
CCK-^INUd 

CO  '0  : = c,5c 

IP  (PKI)  .Lr.3.>,  )Pn^,  12 

r cm:  NOE 

p 1 ( I)  = . . j 

c :nt:  N'JF 

CrNTIMIE 

s:;;::., ..  BEST'AVAIWBIE  COPY 

ravirofly/ifoi  (’’I 

C CM  INCE 

DflV-(:.5*(Pl(l)  ♦Pl(51>)  ♦PflWD/fO.O 


PPL'f'  -EP--  3'NTS  the  NChCTP'ENCIONaL  LOAD.  TC  OBTAIN  LCAO  IN 
PCl-'NPS,  PSO-  ^AS  TO  pp  OIVIOFD  BY  ( HO  Ap  • < /IPP  A ) / ( ><»L  ) * 

SOOT (ALPHA/ (P  J*K)  > 


PSU---  = PA  J*5C  . l‘Dx 

PLC')=ALCOi:  (pOUH) 

PPIM  31'A,PA\/,PSUM,PLCG 

cCP^ATdl*  ,*A7EpAG-  PPE?SUAr=*,El?.p,PY,*LOAC  SLPPOPT_n  = *,rl2.5, 

BX ,*L J& (LOAO) ,Pfe. 3) 

b:i=:.o*(i.c/h(i)»i.i/h(5i)) 

B :^  = c . 0 

f-xrN  7 4 T — > 

'-U  'cl  4.“C»  o 

p:2=p:2x(i.:/h(I)» 

continue 


WmE'EVE?  the  riLH  CAVITATps,  MEATINO.  IS  ASSDMEO  TO  3E  ZERO 


no 

I'^  ('  1(  I ) . £0 . 1.  0 t zc  1 » 70- 
ZCp  CCMIMjP 

p':(:>  = (<;*'n/H(i>-(u*u*O7W(2,a*3.i4iBP2fi54)*<B0i«-P02) 
b:(:»=‘U(:)*cy**2.3 
or  TO  7G7 
7'M  h:(!)=C..1 

T?7  CONTI  MJ- 

7C2  continue 

30(11=3.3 
B C ( , 1 ) = o . 3 

call  rV3(A,51.51,S3,l,xL,0ET,inPT,r,  ,IC,IEP) 


XLC)  HEPPESENTS  The  thep^-ally  GENEOATED  part  OE  H(I),  ANO  Y(I» 
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END 
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